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1.  INTRODUCTION 


The  performance  of  earth-orbiting  satellites  is  profoundly  affected  by  the  high-energy 
particle  environment.  At  high  altitudes,  spacecraft  charging  by  keV  particles  can  produce 
electrical  discharges  on  surfaces  and  insulating  materials  that  couple  to  satellite  electron¬ 
ics.  Impact  of  high-energy  particles  can  degrade  electronic  components,  cause  single  event 
upsets  and  produce  internal  charging.  The  high-energy  particle  environment  has  been 
studied  for  many  years  and  static  models  have  been  constructed  that  specify  the  time- 
averaged  conditions  in  the  radiation  belts.  These  models  have  aided  spacecraft  designers 
in  minimizing  the  deleterious  effects  of  high-energy  charged  particles.  However,  the  radia¬ 
tion  belts  are  very  dynamic,  undergoing  extreme  variations  in  response  to  magnet ospheric 
and  solar  activity.  With  the  coming  of  Air  Force  programs  involving  manned-presence  in 
space,  with  the  placement  of  spacecrzift  in  dynamic  regions  of  the  outer  radiation  belt,  as 
well  as  with  the  development  of  advanced  systems  that  may  be  sensitive  to  the  radiation 
environment,  the  understanding  and  possible  early  warning  of  this  dynamic  activity  is 
critical  to  successful  missions.  Identifying  the  sources  and  sinks  of  high-energy  electrons, 
and  understanding  their  relationship  to  the  observed  variations,  are  the  ultimate  goals  for 
developing  a  dynamic  model  of  the  radiation  belts. 

This  final  report  covers  a  45-month  research  effort,  undertaken  to  investigate  the  dy¬ 
namics  of  radiation-belt  particle  populations  using  an  existing  data  base  provided  by  the 
Lockheed  High  Energy  Particle  Spectrometer  (SC3)  on  the  Space  Charging  AT  High  Al¬ 
titude  (SCATHA)  satellite.  This  instrument  has  been  measuring  the  high-energy  electron 
population  from  47  keV  to  5  MeV  in  the  near-geosynchronous  region.  Because  of  the  high 
angular  resolution  (3°  field  of  view)  of  the  instrument  and  the  radial  movement  of  the 
spacecraft  over  the  L-shell  range  of  5.3  to  8.7,  the  measurements  constitute  a  unique  data 
base  that  idlows  the  interpretation  of  pitch-angle  diffusion  effects  in  the  presence  of  radial 
transport. 

The  main  objective  of  this  data  analysis/interpretation  effort  is  to  set  up  a  prototype 
of  the  data  effort  that  will  be  required  in  order  to  interpret  and  expand  the  single-point 
radiation-belt  measurements  to  be  acquired  by  the  Chemical  Release  and  Radiation  Ef¬ 
fects  Satellite  (CRRES)  into  dynamic  models  of  the  radiation  belts  as  a  whole.  The  long 
term  goal  of  the  CRRES  program  is  to  develop  useful  models  of  the  entire  radiation-belt 
environment.  While  a  static  phenomenological  model  can  be  constructed  without  much 
attention  being  peiid  to  characterizing  the  physical  basis  of  the  model,  such  phenomenolog¬ 
ical  models  suffer  from  the  major  drawback  of  not  being  able  to  characterize  the  dynamic 
changes  of  the  radiation-belt  environment  that  result  from  natural  geomagnetic  activity 
or  from  high-altitude  nuclear  explosions.  Therefore,  if  the  CRRES  data  are  to  be  used 
for  predictive  purposes,  there  is  no  ciltemative  but  to  interpret  the  data  first  in  terms  of 
physical  phenomena  and  subsequently  to  use  them  as  the  basis  for  a  predictive  dynamical 
model  of  the  radiation-belt  environment. 
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Such  a  program  for  CRRES  is  simple  and  effective  in  concept,  but  its  implementation 
is  a  different  matter  since  the  completeness  and  precision  of  the  expected  CRRES  data 
have  not  been  matched  by  any  other  previous  measurements  except  SCATHA.  Thus,  the 
SCATHA  energetic  electron  data  set  is  a  natural  testing  ground  for  the  CRRES  data 
analysis/interpretation  efforts.  Further,  since  the  development  of  such  efforts  requires 
some  lead  time,  the  present  SCATHA  data  analysis  is  the  perfect  choice  as  a  prototype 
for  CRRES  data  analysis  in  1990  and  beyond. 

The  basic  physiciil  requirement  for  implementing  such  a  prototypical  program  is  a 
fundamental  understanding  of  how  radiation-belt  particles  are  transported  in  space  and 
time,  how  their  motion  is  constrained  by  the  Earth’s  magnetic  field,  and  how  their  energy 
is  partitioned  in  the  complex  magnetospheric  plasma  envirionment.  In  short,  we  need 
to  “keep  book”  on  the  radiation-belt  particles  with  as  much  detEul  as  we  can  afford,  in 
order  to  understand  their  dynamical  interactions.  The  formulation  of  such  “bookkeeping” 
in  terms  of  the  behavior  of  statistical  aggregates  has  been  developed  in  the  late  1960’s 
as  radiation-belt  diffusion  theory  [e.g..  Reference  1].  In  the  statistical  picture,  energetic 
particles  are  transported  into  or  out  of  any  imaginary  volume  of  space  within  the  radiation 
belts  by  the  cumulative  effects  of  interactions  with  fluctuating  electromagnetic  fields.  Such 
transport  can  be  from  one  L-shell  volume  into  another  and/or  along  the  magnetic  field  into 
the  atmosphere.  The  diffusion  formalism  seeks  to  describe  these  transport  phenomena  by 
parametrizing  the  rates  of  transport,  in  its  various  modes,  in  terms  of  diffusion  coefficients. 
Thus,  in  the  diffusion  theoretic  formulation  of  radiation-belt  dynamics,  the  primary  task  is 
to  determine  the  diffusion  coefficients  from  data,  rather  than  construct  a  global  predictive 
model  based  on  first  physical  principles.  A  finite  set  of  diffusion  coefficients  comprise 
the  key  parameters  in  the  diffusion  equation  which  specifies  the  dynamical  state  of  the 
radiation  belt.  The  chief  goal  of  our  data  analysis/interpretation  effort  is,  therefore,  the 
characterization  of  the  diffusion  coefficients  for  the  SCATHA  data  set  on  outer-radiation- 
belt  electrons.  This  has  been  accomplished  and  discussed  in  detail  in  our  papers,  included 
in  the  appendices  of  this  report,  resulting  from  this  study. 

The  adoption  of  a  diffusion  point  of  view  implies  that  the  particle  distribution  be 
specified  with  respect  to  a  set  of  phase-space  coordinates  or  parameters:  any  of  which 
may  vary  during  the  physical  interactions  which  lead  to  diffusion.  The  ideal  coordinate 
frame  for  the  description  of  charged-particle  motion  and  diffusion  in  the  radiation  belts  is 
the  space  of  three  adiabatic  invariants.  Each  invariant  corresponds  to  one  form  of  cyclic 
motion  in  a  static  electromagnetic  field.  The  deviations  from  purely  cyclic  motion  in  the 
presence  of  time  varying  fields  lead  to  particle  diffusion  with  respect  to  one  or  more  of 
the  adiabatic  invariants.  In  the  inner  radiation  belts  the  adiabatic-invariant  formalism  is 
especially  valuable  because  there  is  a  simple  mapping  of  the  adiabatic-invariant  space  onto 
the  configuration  space  defined  by  the  L  shells:  the  mix  of  particles  on  any  given  L-shell 
is  invariant  with  respect  to  longitude  under  static  conditions. 

For  the  outer  radiation  belt,  and  for  geomagneticsilly  active  times,  however,  the  mat- 


ter  is  not  so  simple.  The  primeiry  problem  is  that  the  definition  of  an  “L-shell”  is  not  very 
meaningful  under  such  circumstances  because  the  mapping  of  adiabatic-invariant  space 
onto  the  configuration  space  (i.e.,  radial  distance,  longitude,  and  latitude)  is  not  invari¬ 
ant  vi'ith  respect  to  longitude.  The  mapping  also  changes  whenever  the  geomagnetic  field 
varies  w’ith  time  on  a  scale  that  cannot  be  described  as  “fluctuation”.  Consequently,  L  is 
not  a  good  coordinate  with  which  to  characterize  radial  diffusion  “across  L-shells”.  The 
concept  of  adiabatic  invariants,  however,  is  still  well-defined.  The  breakdown  of  the  map 
ping  between  invariant  ajid  configuration  spaces  can  lead  to  coupling  between  the  modes 
of  diffusion;  a  phenemenon  which  can  be  treated  correctly  only  with  adiabatic-invariant 
coordinates.  Therefore,  the  generalized  description  of  diffusion  in  adiabatic-invariant  space 
is  the  proper  choice  to  address  radiation-belt  diffusion  in  the  outer  regions  at  geomagnet- 
ically  active  times.  A  prerequisite  to  our  goal  of  data  interpretation  is  the  capability  of 
transforming  between  the  two  sets  of  coordinates  for  realistic  geomagnetic  fields  relevant 
to  various  active  conditions.  The  difficulty  of  generating  an  efficient  transformation  from 
configuration  space  to  a  fully  adiabatic-invariant  coordinate  frame  is  the  primary  reason 
why  there  has  been,  until  this  time,  no  method  for  routinely  characterizing  diffusion  in  the 
analysis  of  outer-radiation-belt  data  for  geomagnetically  active  times. 

The  diffusion  analysis  has  been  attacked  on  several  fronts.  Our  efforts  for  the  first 
nine  months  of  this  study  were  mostly  directed  at  a  diffusion  description  in  the  static 
geomagnetic  field.  the  same  time  wc  proceeded  to  develop  basic  tools  to  cope  with 
changes  in  the  magnetic  field  during  active  times.  A  first  step  toward  a  fully  dynamic 
diffusion  model  was  the  acquisition  for  this  study  of  a  divergence-less  realistic  geomagnetic 
field  valid  for  high  Kp  values.  That  requirement  seems  to  be  fully  met  in  the  form  of 
the  Tsyganenko- Usmanov  model  [References  2  and  3]  as  modified  by  Dr.  David  Stern 
of  NASA  Goddard  Space  Flight  Center  [private  communication].  Ihe  second  step  was 
to  acquire  the  capability  to  compute  adiabatic  invariants  in  this  model.  Only  efficient 
computation  of  the  third  invariant  was  needed  since  the  computation  for  the  first  and 
second  invariants  is  routine.  The  third  adiabatic  invariant,  which  is  proportional  to  the 
area  vector  integral  of  the  magnetic  flux  enclosed  by  the  particle  drift  orbit  around  the 
Earth,  is  rather  cumbersome  and  time-consuming  to  compute  on  a  routine  basis.  We  are 
not  aware  of  previous  attempts  to  do  so  on  a  routine  basis.  An  important  development 
of  the  present  effort  was  that  we  have  been  able  to  develop  an  accurate  transformation  to 
reduce  the  two-dimensionzd  integration  for  the  third  invariant  into  a  simple  one-dimensionad 
integration,  thus  opening  the  way  to  analysis  of  outer-belt  diffusion  at  active  times,  as  well 
as  for  other  applications,  such  as  for  the  proton  belts. 

Detauled  descriptions  of  these  achievements  are  described  in  our  interim  report  for 
the  first  nine  months  [Reference  4].  These  include  a  brief  survey  and  description  of  the 
SCATH.4  SC3  electron  data  along  with  characterization  of  some  active  periods  of  interest; 
modehng  of  the  temporal  variations  by  mapping  the  phase-space  distribution  functions 
in  terms  of  observable  parameters  at  fixed  first  and  second  adiabatic  invariants;  a  com- 
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putationftUy  efficient  procedure  for  calculating  the  third  adiabatic  invariant  and  the  drift 
motion  of  the  particles;  and  the  first  phase  of  fitting  a  period  of  SC3  electron  data  with 
an  eigenfunction  formulation,  combining  both  radial  and  pitch-angle  diffusion. 

Following  this  initial  phase  of  our  study,  a  general  solution  of  the  bimodal  (radial  and 
pitch-angle)  diffusion  equation  for  the  radiation  belts  was  developed  with  special  regard  to 
the  requirements  of  satellite  radiation  belt  data  analysis.  Utilizing  a  period  of  SCATHA 
SC3  electron  data,  the  radial  and  pitch-angle  diffusion  coeffients  were  determined  and 
compared  to  previous  results  as  a  first  test  of  our  solution  This  was  presented  in  our  first 
paper  [Reference  5],  which  can  be  found  in  Appendix  A  and  will  be  hereafter  referred  to 
as  Paper  I.  These  results  were  expanded  on  in  our  second  paper  [Reference  6i,  which  also 
included  a  discussion  of  issues  concerning  dynamic  modeling.  This  work  will  be  referred 
to  as  Paper  II  and  can  be  found  in  Appendix  B.  The  development  of  a  nonlinear  response 
model  of  radiation  belt  particles  to  wave  growth  and  particle  precipitation  [Reference  7 
was  partially  sponsored  by  this  study,  can  be  found  in  Appendix  C  and  will  be  referred  to  as 
Paper  III.  As  we  have  progressed  in  this  study  towards  a  computationally  efficient  dynamic 
model  of  the  outer  electron  belt,  we  have  refined  our  earlier  work  on  the  solution  of  the 
bimodal  diffusion  equation  and  its  uniqueness  into  our  fourth  paper  (Paper  IV)  [Reference 
8].  This  work  can  be  found  in  Appendix  D.  In  Section  II  of  this  final  report,  the  bimodal 
diffusion  equation  and  our  general  solution  is  presented,  as  developed  in  Papers  I,  II  and  IV. 
This  general  solution  can  serve  as  the  basis  of  a  new  approach  to  dynamic  modeling  of  the 
radiation  belts,  having  the  capability  to  include  temporal  variations,  diffusion  and  trapping. 
The  next  four  sections  summarize  the  results  of  each  of  our  papers  from  this  study.  Our 
results  and  conclusions  from  this  Space  Radiation  Test  Model  Study  are  presented  in 
Section  VII.  Recommendations  for  implementation  of  dynamic  modeling  in  the  CRRES 
time  frame  are  given  in  Section  VIII. 
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II.  A  GENERAL  SOLUTION  OF  BIMODAL  DIFFUSION 


In  order  to  organize  the  d'-^cussion  in  this  report,  we  briefly  summarize  the  salient 
features  of  simultaneous  bimodal  diffusion  theory  [References  1,  9  and  10'  and  its  specific 
form  developed  for  efficient  data  representation  in  Paper  I. 

By  introducing  an  approximate  invariant  valid  over  a  sufficiently  thin  slice  of  L-space, 
Walt  Reference  9  showed  that  the  bimodal  diffusion  equation  can  be  formulated  in  con¬ 
figuration  and  pitch-angle  coordinates  where  x  =  cos  a.  Such  a  formulation  is 

particularly  convenient  for  data  applications  because  the  coordinates  are  simple,  although 
it  needs  to  be  emphasized  that  its  validity  is  restricted  to  slices  in  L-space,  over  which 
energy  gained  or  lost  from  radial  diffusion  does  not  impact  Walt’s  invariant.  As  discussed 
in  Paper  I  the  simultaneous  multi-mode  diffusion  equation  for  the  electron  distribution 
function  /  can  be  written 


dt 


dx^ 


X  dx 


A 

^  dL 


n  +  7 


dL 


(1) 


where  p  and  y  are  considered  to  be  constant  parameters.  In  Paper  1,  we  showed  that  if  the 
pitch-angle  and  radial  diffusion  coefficients,  and  Dpt,  respectively,  take  on  certain 
generic  forms 

=  D^JE) 

zz\  )  ^2) 

Dll  =  i{E)L^ 

then  the  general  solution  to  the  bimodad  diffusion  equation  can  be  expanded  in  terms  of 
eigenfunctions  of  various  Bessel  forms.  This  convenient  property  forms  our  basic  principle 
to  characterize  the  state  of  the  outer  electron  belt  in  a  piecewise-continuous  faishion  in 
terms  of  a  data  base  measured  along  a  satellite  trajectory  L{t). 

The  general  solution  consists  of  the  linear  sum  of  three  modal  expansions:  (1)  single¬ 
mode  pitch-angle  diffusion.  (2)  single-mode  radial  diffusion  and  (3)  simultaneous  bi-mod.al 
diffusion.  Using  the  device  of  the  Kronecker  delta  (fno.  the  three  modal  expansions  can  be 
written  in  terms  of  a  single  eigenfunction  expansion 


£  M  Joih,,x/Xc] 

CnO  -  (  t  -  CriO  I -  •  - r-Tt — ^ 

Xc  Jl(kn) 


(c„(l 


where  is  the  zero  of  the  Bessel  function  Jo  and  Xc  = 
functions  are  defined  as  follows; 


(3) 

cosQc.  The  parameters  and 


e  =  3(1  -  |)/(p  -  2) 
p  1  3 

^  =  _  _  _  _  ^  -ly 


(•i) 
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Y„(L)  =  L’  IL^)]^-[h^,[\l3„  \  IL^)]^  (5) 

Z„(i)  =  (6) 

J  is  the  regular  Bessel  function,  I  and  K  are  the  modified  Bessel  functions,  and  N  is  the 
Neumann  function.  The  usual  notation  [p]>  is  defined  in  this  case  by 


[gh  =9  if  >  0 


[5]>  =  1  if  <n  <  0 


(7) 


and  analogously  for  [5j<.  The  separation  constants  and  parauneters  of  the  solution 

are  related  by  the  algebraic  relations 


in'  =  (1  -  ^no)T-'  -  (8) 

=  -Dzzkilx\  (9) 

The  parameter  vaJue  A  =  4  (i.e.,  p  =  10)  is  chosen  to  conform  with  the  consensual  L- 
dependence  of  the  radial  diffusion  coefficient  [Reference  11].  From  Paper  I  the  parameter  u 
is  set  to  2  (i.e.,  e  =  0)  for  ease  of  computing  the  Bessel  functions.  Other  values  of  v  are  easily 
accommodated  at  the  expense  of  computational  speed.  Our  work  so  far  has  demonstrated 
the  economy  of  using  u  =  2,  but  further  refinements  can  be  readily  accommodated  in  our 
procedure.  Equation  (8)  will  be  referred  to  as  the  separation  equation  because  the  entire 
separable  diffusion  differential  equation  is  reduced  by  the  eigenfunction  expa  asion  into  this 
constraint  on  the  diffusion  parameters.  It  is  seen  from  (9)  that  the  single  pitch-angle  mode 
time  constants  are  negative  as  expected,  since  single-mode  pitch-angle  diffusion  always 
leads  to  decay.  For  radial  diffusion,  however,  the  time  constants  can  be  either  positive 
or  negative  over  a  finite  time  interval.  This  is  due  to  the  fact  that  the  phase  space  density 
in  the  finite  L  domain  can  rise  or  fall  in  time  in  response  to  particle  transport  into  and 
out  of  the  adjacent  volumes.  With  the  sign  of  each  not  predetermined,  the  separation 
constant  /?*  can  be  positive  or  negative.  This  is  the  origin  of  the  (J, /)  and  [N.K)  options 
of  (5)  and  (6).  This  general  solution  to  the  simultaneous  bimodaJ  diffusion  equation  and 
the  data  representation  were  introduced  and  utilized  in  three  of  the  four  papers  to  be 
discussed  below  as  part  of  our  test  model  study. 
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III.  FIRST  RESULTS  FROM  SIMULTANEOUS  BIMODEL  DIFFUSION 


The  generaJ  solution  of  the  simultaneous  bimodal  (radial  and  pitch-angle)  diffusion 
equation  for  the  outer  electron  radiation  belts  developed  in  the  previous  section  was  first 
put  forth  in  Paper  I  [Reference  5]  and  is  included  in  this  report  as  Appendix  A.  The  so¬ 
lution  was  developed  with  special  regard  for  the  requirements  of  satellite  radiation  belt 
data  analysis  and  is  used  to  test  the  bimodal  theory  of  outer  electron  belt  diffusion  by 
confronting  it  with  satellite  data.  Satellite  observations,  usually  over  finite  volumes  of 
{L.  t)  space,  are  seldom  sufficient  in  space-time  duration  to  cover  the  relaxation  to  equi¬ 
librium  of  the  entire  radiation  belt.  Since  time  scales  of  continuous  data  coverage  are 
often  comparable  to  that  of  radiation  belt  disturbances,  it  is,  therefore,  inappropriate  to 
apply  impulsive  semi-infinite  time  response  solutions  of  diffusion  theory  to  interpret  data 
from  a  finite  .vindow  of  (L,  t)  space.  Observational  limitations  indicate  that  appropriate 
solutions  for  the  interpretation  of  satellite  data  are  general  solutions  for  a  finite-volume 
boundary  value  problem  in  bi-modal  diffusion.  In  the  paper  we  tested  such  a  solution  as 
the  prime  candidate  for  comprehensive  radiation-belt  dynamic  modeling  by  applying  the 
solution  and  developing  a  method  of  analysis  to  radiation  belt  electron  data  obtzdned  by 
the  SC.4TH.4  satellite  at  moderate  geomagnetic  activity.  The  success  of  our  first  results 
and  the  generality  of  our  solution  indicate  its  promise  as  a  new  approach  to  dynamic  mod¬ 
eling  of  the  radiation  belts.  Such  an  approach  should  be  based  on  our  bimodal  diffusion 
solution  and  on  satellite  data  taken  on  a  routine  global  basis. 

Our  test  data  base  from  the  SCATHA  SC3  experiment  was  discussed  in  Paper  I, 
followed  by  aji  in-depth  introduction  to  bimodal  diffusion  in  the  outer-belt  region.  When 
a  period  of  data  from  day  165,  June  13,  1980,  was  fit  to  the  general  solution  expansion 
(Equations  3  through  9),  exclusive  of  the  separation  equation  (8),  the  electron  distribution 
function  was  characterized  by  the  set  of  parameters  (o^,  bn,  c„.tn,  Tn,  3n)-  The  phase  space 
volume  for  this  fit  covered  the  region  (5.3  <  L  <  6.7;  A<  =  5  hrs.;  200  keV  <  E  <  2000 
ke\':  Oc  <  Q  <  90°).  where  q  is  the  pitch-angle,  Oc  is  the  edge  of  the  loss  cor.t  and  E  is 
the  energy.  It  wa.s  found  that  five  orders  of  eigenfunctions  (n  =0-4)  per  energy  channel, 
i.e..  30  parameters  per  energy  channel,  were  sufficient  to  characterize  the  full  data  base  of 
5040  data  points.  Note  that  t„  and  3„  were  independently  derived  from  the  fit  and  (8) 
was  not  imposed.  No  iterations  were  made  and  a  value  of  ^  was  determined  for  each  order 
from  (8).  i.e.  a  set  of  From  Equation  (9)  the  pitch-angle  coefficient  £)„  can  be  derived 
from  the  data  set  for  each  energy  channel.  Figure  1  shows  the  derived  Du  together  with 
the  fit  errors  as  functions  of  ei.  rgy.  In  the  solution,  Pn  w'as  assumed  independent  of  L  in 
order  to  achieve  separability.  From  the  figure,  it  is  seen  that  Z)„  is  constant  with  energy, 
except  for  the  446  ke\’  channel.  No  iterations  were  made  in  this  early  work.  In  Paper  1 
these  values  of  were  found  to  agree  with  other  determinations  giving  a  pitch-angle 
diffusion  lifetime  of  about  5  days. 

The  derivation  of  the  radial  diffusion  coefficient  serves  two  roles  of  providing  for  phys- 


ical  data  in  the  outer  belt,  which  can  be  compared  with  other  determinations  from  direct 
particle  observations,  and  of  testing  independently  the  adequacy  of  our  representation  to 
the  data.  Figure  2  shows  a  scatterplot  of  Du  versus  energy  for  L  —  5.3.  As  mentioned 
above  the  scatter  values  of  Du  result  from  our  deliberately  not  enforcing  (8)  in  an  itera¬ 
tive  manner.  A  worst  case  variance  of  the  fit  is  obtained  from  the  set  of  at  each  order. 
It  was  found  that  these  order-dependent  values  of  ^  clustered  closely  around  those  deter¬ 
mined  by  other  workers,  an  indication  that  the  representation  is  reaisonable.  The  results 
of  Paper  I  were  encouraging  in  that  solutions  of  simultaneous  bimodal  (pitch-angle  and 
radial)  diffusion  theory  were  shown  to  have  the  potential  to  provide  an  efficient  description 
of  a  SCATHA  outer-belt  electron  data  set.  The  description  was  achieved  in  essentially 
the  initial  pass  without  a  systematic  fitting  procedure.  These  initial  results  established 
plausibity  for  the  approach. 


9 


Radiol  Diffusion  Cocff.  (day**) 


10 


Energy  (keV) 


Figure  2.  Radial  diffusion  coefficient  as  function  of  energy  at  Z,  =  5.30.  The  spread  of 
values  for  each  energy  is  a  measure  of  the  worse-case  variance  in  its  determi¬ 
nation. 


IV.  HIGH-ENERGY  OUTER  RADIATION  BELT  DYNAMIC  MODELING 


Paper  II  [Reference  6],  which  can  be  found  in  Appendix  B,  originated  from  two  papers 
we  presented  at  the  Conference  on  High  Energy  Radiation  Background  in  Space  in  antici¬ 
pation  of  the  dynamic  modeling  objective  of  the  upcoming  CRRES  program.  We  described 
the  good  results  of  our  initial  attempt  at  constructing  the  dynamics  of  the  outer  electron 
radiation  belt  by  combining  diffusion  theory  with  data  as  found  in  Paper  I  and  summa¬ 
rized  in  the  previous  section,  as  well  as  discussed  issues  concerning  dynamic  modeling  in 
general.  The  latter  was  pertinent  to  the  conference  because  dynamic  modeling  should  pro¬ 
vide  a  means  to  map  the  particle  phase-space  distributions  away  from  the  satellite  track. 
However,  the  uniqueness  of  such  an  off-trajectory  projection  needs  to  be  determined  in 
light  of  the  fitting  errors  involved  in  the  data  representation.  These  issues  were  brought 
up  here  and  analyzed  in  more  detail  in  Paper  IV,  summarized  in  Section  VI  below. 

The  formulation  of  the  simultaneous  bimodal  diffusion  problem  and  our  approach  is 
introduced  in  a  different  way  from  that  of  Paper  I.  This  may  provide  a  better,  or  at  least 
different,  perspective  of  our  solution  representation.  Further,  we  introduced  the  fact  that 
the  radial  portion  of  the  solution  is  explicitly  a  function  of  t  by  expressing  L  as  a  function 
of  i.  We  may  rewrite  the  radial  portion  of  Equation  (3)  as 

h{L{i),  t)  =  {Lit))  +  bnZn  (!(<))]  (10) 

Utilizing  L{t)  minimizes  the  mixing  of  space-time  in  our  fitting  For  the  neau’-geosynchronous 
SCATHA  satellite  covering  only  a  small  range  of  L,  the  effect  of  space-time  mixing  is  mini¬ 
mal.  This  allowed  us  to  complete  the  first  test  case  shown  in  Paper  I.  However,  for  the  more 
comprehensive  situation  of  highly-elliptic  orbits  such  as  that  for  CRRES,  which  attempts 
to  specify  the  radiation  belts  over  a  wide  range  of  X,  the  trajectory  L{i)  is  a  nonlinear 
function  of  i,  requiring  a  nonlinear  fit  for  the  determination  of  the  diffusion  parameters. 
This  is  illustrated  in  Figure  3,  which  shows  plots  of  h{L)  for  three  situations:  (a)  X  and  t 
independent,  (b)  SCATHA  orbit  on  which  X  is  nearly  independent  of  t,  and  (c)  CRRES- 
like  orbit  on  which  X  is  strongly  dependent  on  t.  The  strong  distortions  on  panel  (c)  for 
hiL{i).t)  relative  to  the  same  function  h{L,t)  on  panel  (a)  signifies  that  dynamic  model¬ 
ing  in  conjunction  with  a  satellite  data  base  involves  much  more  than  obtaining  solutions 
to  diffusion  theory.  Using  L{i)  results  in  a  better  fitting  scheme  and  supports  unique  off- 
trajectory  mappings.  Paper  II  also  briefly  touched  on  the  problem  of  magnetic  and  electric 
X-shell-splitting,  which  needs  to  be  studied  further  so  that  it  can  be  taken  into  account  in 
the  dynamic  models. 

The  discussion  in  Paper  II  showsed  that  dynamic  modeling  differs  fundamentally 
from  radiation  belt  analysis  in  that  a  theoretical  representation  of  the  on-trajectory  data 
set  is  required  to  predict  fluxes  outside  of  the  satellite  orbit.  The  fjiithfulness  of  the 
projection  depends  on  whether  the  representation  is  a  good  description  of  the  physical 
diffusive  processes  seen  in  the  on-trajectory  data.  Our  test  case  study  yielded  results 
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Figure  3.  Illustration  of  the  nonlinear  nature  of  fitting  the  radial  representation  to  data. 

In  all  panels,  the  function  h{L,t)  is  plotted  with  the  same  set  of  arbitrarily 
chosen  constant  parameters.  The  solid  lines  are  contours  of  h{L,t)  and  the 
dashed  lines  are  L{t).  Panel  (a)  is  for  an  independent  {L.i)  grid.  Panel(b) 
is  for  L{i)  specified  by  SCATHA  orbit  and  panel  (c)  is  for  the  CRRES  orbit. 
Note  the  drastic  change  of  time  scale  required  in  order  for  the  CRRES  case  to 
have  similar  contours  for  the  h{L,t)  function  as  the  nearly  identical  panels  (a) 
and  (b). 
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consistent  with  basically  all  other  measurements  in  the  moderately  active  outer  electron 
radiation  belt.  The  resolution  of  a  number  of  important  issues  brought  by  examination  of 
the  test  case  awaits  further  work  before  operational  dynamic  models  can  be  constructed. 
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V.  A  NONLINEAR  MODEL  OF  WAVE-PARTICLE  INTERACTIONS 
IN  THE  TRAPPED  RADIATION  BELTS 


The  interactions  of  waves  and  trapped  electrons  are  described  in  Paper  III  [Reference 
7],  which  is  attached  in  Appendix  C,  by  a  closed  system  of  three  nonlinear  coupled  equa¬ 
tions.  The  model  has  applications  to  many  aspects  of  wave-particle  interactions  in  the 
magnetosphere.  Nonlinear  numerical  solutions  pertinent  to  auroral  pulsations  have  been 
computed  for  realistic  ranges  of  physical  parameters.  The  results  confirm  that  precipitat¬ 
ing  pulsations  can  be  initiated  by  either  an  injection  of  energetic  trapped  particles  or  an 
increase  in  the  plasma  ionization  desity.  Relaxation  of  the  system  results  in  cyclic  phase 
trajectories  about  a  (quasi  equilibrium)  point  that  is  not  located  at  the  initial  equilibrium 
point.  That  observation  may  explain  the  tendency  of  auroral  pulsations  to  appear  super¬ 
imposed  upon  an  increase  in  the  general  rate  of  electron  precipitation.  The  results  also 
imply  that  precipitation  pulsations  are  most  readily  initiated  by  an  encounter  of  drifting 
energetic  electrons  with  a  region  of  enhanced  cold  plasma  density. 

The  interaction  between  waves  and  trapped  particles  in  the  magnetosphere  is  a  key 
problem  in  magnetospheric  physics  [Reference  12],  whose  basic  elements  are  considered  to 
be  well  understood.  However,  the  complexity  of  a  full  treatment,  based  on  first  principles, 
of  the  diffusion  of  particles  into  the  loss  cone  and  the  response  of  waves  to  particle  diffusion 
has  restricted  quantitative  studies  to  a  handful.  Davidson  et  al.  [Reference  13]  utilize  the 
SCATHA  SC3  data  base  to  test  the  theoretical  trapping  limit  proposed  by  Schulz  and 
Davidson  [Reference  14).  The  results  of  the  study,  under  a  different  sponsorship  thajj  that 
of  this  report,  were  found  to  be  consistent  with  the  proposed  model  and  the  numerical 
resiilts  agreed  with  the  earlier  predictions  of  Reference  8.  In  this  study  the  pitch-angle 
distribution  was  fit  to  determine  the  parameters  of  the  shape.  The  representation  of  our 
bimodal  diffusion  solution  also  determines  shape  parameters  of  the  pitch-angle  distribution. 
Thus  our  representation  can  monitor  the  state  of  the  radiation  belt  and  indicate  when  the 
particle  populations  are  approaching  isotropy,  i.e.,  approaching  the  trapping  limit. 

The  modeling  applications  for  which  the  solutions  of  Paper  III  can  be  applied  are 
many.  In  fact,  the  physical  principles  and  the  code  used  in  the  paper  are  applicable  with 
virtually  no  modifications  to  radiation  belt  particle  response  to  wave-particle  interactions 
and  flux  limiting.  The  nature  of  the  detailed  wave  modes  enters  only  in  the  determination 
of  the  physical  parameters.  This  wave-particle  interaction  model  may  play  a  role  in  our 
dynamic  modeling  in  the  future. 
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VI.  TOWARDS  DYNAMIC  MODELING  OF  THE  OUTER 
ELECTRON  RADIATION  BELT 

We  show  in  Paper  I\’  [Reference  8j,  which  is  in  Appendix  D,  that  the  general  solution 
of  a  simple  form  of  the  simultaneous  bimodad  (radial  and  pitch-angle)  diffusion  equation 
can  represent  the  SCATHA  outer  belt  electron  distribution  to  a  high  degree  of  accuracy  for 
periods  up  to  10  hours.  The  data  representation  requires  diffusion  parameters  that  are  in 
agreement  with  all  previous  experimental  and  theoretical  determinatiuns  of  such  diffusion 
coefficients.  Such  a  representation  satisfies  the  basic  requirements  towards  the  dynamic 
modeling  of  the  outer  electron  radiation  belts  for  X  <  7  at  quiet  and  moderate  geomagnetic 
conditions.  For  L  ~~  1  the  representation  has  difficulties  with  “butterfly”  distributions  that 
may  signify  energetic  electron  encounters  with  the  magnetopause.  It  is  demonstrated  with 
a  typical  example  of  computer  simulation  of  electron  encounters  with  the  magnetopause 
that  the  effects  of  the  encounter  may  be  highly  anisotropic  in  configuration  or  in  invariant 
space. 

This  paper  addresses  many  of  the  issues  still  remaining  in  our  test  model  study  for 
developing  the  tools  needed  for  dynamic  modeling  of  the  radiation  belts.  Because  the  res¬ 
olution  of  several  of  these  issues  axe  new  in  this  study,  we  therefore  include  major  sections 
of  the  paper  where  they  are  relevant  to  our  results  and  conclusions.  At  least  three  major 
issues  require  attention  if  the  bimodal  diffusion  approach  to  dynamic  modeling  is  to  be 
pursued.  First,  we  must  establish  a  systematic  iterative  procedure  to  achieve  a  fit  of  the 
data  consistent  with  all  the  constraints  of  the  solution.  This  requirement  involves  not  only 
specification  of  the  accuracy  of  the  fit  in  four- dimensional  phase  space  but  also  satisfaction 
of  diffusion  parameter  constraints  imposed  by  bimodal  theory.  Second,  the  form  of  the 
bimodal  diffusion  theory  used  was  computationally  efficient  but  its  validity  was  recognized 
to  be  restricted  to  a  sufficiently  thin  slice  of  X-space  [Reference  9].  Such  a  property  is  not 
incompatible  with  dynamic  modeling  with  a  satellite  data  base  since  satellite  data  are  nec¬ 
essarily  restricted  in  (X,t)  space.  Nevertheless,  it  is  crucial  to  determine  the  observational 
extent  of  X-space  over  which  the  version  of  bimodal  diffusion  theory  is  Vcdid.  Third  and 
most  important,  the  degree  of  uniqueness  of  the  data  representation  must  be  established 
since  dynamic  modeling  implies  the  projection  of  electron  fluxes  from  measurements  on 
a  satellite  trajectory  to  other  regions  of  phase  space  where  direct  measurements  are  not 
available.  Note  that  the  solution  representation  of  the  theory  is  necessarily  unique  every¬ 
where  if  the  on-trajectory  specification  is  unique.  However,  the  data  fit  to  the  solution  is 
over  a  trajectory  on  which  X  and  i  are  mixed;  therefore,  the  on-trajectory  specification 
is  not  necessarily  unique  insofar  as  the  separation  of  X  and  t  effects  are  concerned.  In 
order  to  insure  that  our  iterative  procedure  for  untangling  X  and  t  effects  does  not  alias 
the  projection  to  other  regions  of  space-time,  we  must  at  least  demonstrate  that  our  data 
representation  yields  flux  projections  independent  of  the  projection  path. 

In  Paper  I\'.  we  have  addressed  the  above  three  issues.  We  have  demonstrated  that 
simultaneous  bimodal  diffusion  theory  is  an  efficient  representation  of  SCATHA  outer-belt 
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electron  data  in  the  configuration  space  range  (5  <  L  <  7;  A/  <  10  hrs.)  for  Day  165. 
1980  and  Day  110,  1979.  Iterative  procedures  have  been  devised  to  improve  the  fit  and  to 
satisfy  diffusion  parameter  constrciints  that  were  not  enforced  in  Paper  I.  The  result  is  a 
reduction  of  parameters  (by  ~  15%)  to  achieve  better  fits  over  a  much  larger  data  base  (by 
~  100%).  We  were  also  able  to  demonstrate  that  fit  errors  do  not  impact  the  uniqueness 
of  off-trajectory  flux  projections  based  on  the  solution  representation. 

In  Paper  I  it  was  found  that  five  orders  of  eigenfunctions  (n  =0-4)  per  energy  chan¬ 
nel,  i.e.,  30  parameters  per  energy  channel,  were  sufficient  to  characterize  a  data  base  of 
5040  data  points.  For  Paper  4,  we  included  the  effects  of  the  separation  equation  (8)  in 
constraining  the  diffusion  parameters  (t„,/3„,^).  Note  that  the  diffusion  parameter  ^  is 
independent  of  order  n.  whereas  the  others  are  order  dependent.  In  Paper  I,  Tn  and  /3„ 
were  independently  derived  from  the  fit  and  (8)  was  not  imposed.  No  iterations  were  made 
and  a  vzJue  of  (  was  determined  for  each  order  from  (8),  i.e.  a  set  of  It  was  found 
that  these  order-dependent  values  of  (  clustered  closely  around  those  determined  by  other 
workers,  an  indication  that  the  representation  is  reasonable.  Here  we  used  the  constraint 
(8)  as  the  focus  for  an  iterative  fitting  scheme  in  determining  the  best  order-independent 
value  of  Since  (8)  was  imposed,  (t„,4)  and  0n  were  no  longer  independent,  thus  the 
number  of  required  parameters  per  energy  channel  was  reduced  by  4  in  this  work.  The 
iteration  also  improved  the  fit  to  the  data  set.  The  process  of  fitting  a  data  set  to  the  rep¬ 
resentation  begins  with  decomposing  the  electron  distribution  functions  into  pitch-angle 
eigenfunctions  and  their  amplitudes.  We  have  demonstrated  in  Paper  I  that  this  proce¬ 
dure  is  highly  accurate.  We  shall  focus  our  discussion  below  on  tne  decomposition  of  the 
pitch-angle  amplitudes  into  eigenfunction  expansions  in  L. 

We  have  devised  a  nonlinear  fitting  procedure  keyed  to  the  iterative  determination  of 
the  order-independent  radial  diffusion  parameter  $  in  the  separation  equation  (8).  In  the 
first  pass,  the  procedure  described  in  Paper  I  is  applied  and  the  separation  equation  (8) 
is  solved  for  each  order,  resulting  in  a  set  of  order-dependent  parameters  Iterations 
begin  with  the  choice  of  an  interim  mean  value  ^  =  (^n)  where  the  averaging  is  weighed  by 
the  eigenfunction  amplitudes  (an,  in)-  The  interim  value  ^  is  enforced  in  (8)  to  determine 
new  values  of  the  parameters  {Tn,0n)  in  a  new  fit.  With  the  new  parameters,  a  new  set  of 
^n  is  determined  by  solving  (8)  exactly  and  another  iteration  cycle  begins.  The  iterative 
procedure  continues  until  the  residual  error  in  the  equality  between  the  left  and  right  hand 
side  of  (8)  can  no  longer  be  improved  upon.  The  final  value  of  ^  is  taken  to  be  the  fit  value 
of  the  diffusion  parameter.  This  is  a  nonlinear  iterative  procedure;  therefore,  there  is  no 
a  priori  guarantee  that  it  is  convergent.  In  practice,  it  converges  because  the  set  of  in 
the  first  pass  is  closely  clustered  to  begin  with.  Figure  4  shows  an  example  of  the  degree 
to  which  the  constraint  (8)  is  obeyed  by  the  initial  and  final  iterations  of  the  procedure. 

The  iterative  scheme  has  been  applied  to  ~  4  hours  of  SCATHA  electron  belt  data 
for  Day  165,  1980  and  to  9  hours  of  data  for  Day  110,  1979.  Figures  5  and  6  show  three- 
dimensional  plots  of  the  data  and  the  fit  reconstructions  for  these  two  days,  respectively.  It 
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should  be  noted  that  reconstruction  errors  for  both  days  are  primarily  concentrated  in  the 
high-X  and  q  ~  90“  region.  We  found  that  this  region  is  rich  in  “butterfly"  distributions 
[References  15  and  16].  These  will  be  discussed  in  more  detail  below.  Figures  7  and  8  show 
the  distributions  of  the  fit  errors  for  each  data  point  in  the  initial  and  final  passes  of  the 
iterations  for  the  two  days,  respectively.  The  data  points  associated  with  the  large-error 
tail  of  the  error  distributions,  Jess  than  1000  data  points  in  each  data  set,  are  agzun  related 
to  “butterfly”  distributions,  although  a  substantial  portion  of  the  large  fit  errors  are  also 
contributed  to  by  some  noisy  data  in  the  data  base.  Unsmoothed  raw  data  is  used  to  tally 
the  percentage  error  in  these  figures.  Figure  9  shows  the  comparison  of  radial  diffusion 
coefficients  determined  in  this  work  with  all  other  known  observational  and  theoretical 
determinations  [References  17  through  29].  The  mottled  band  of  Du  values  in  the  figure 
corresponds  to  this  work.  The  spread  of  the  band  is  due  to  energy  dependence  of  the 
diffusion  coefficient.  The  larger  cross-hatched  band  of  Du  values  was  determined  in 
Paper  I  and  is  due  not  only  to  the  spread  of  energy  channels  but  also  primarily  to  the 
dependence  of  ^  on  order.  It  is  seen  that  our  enforcement  of  (8)  allows  for  a  much  more 
precise  determination  of  the  radial  diffusion  coefficient.  Energy  dependence  is  not  always 
addressed  in  other  work. 

For  a  data  repesentation  scheme  to  qualify  as  a  dynamic  modeling  scheme,  it  must 
demonstrate  the  ability  to  uniquely  predict  (or  obtain  by  projection),  the  off-trajectory 
electron  distributions,  from  a  set  of  on-trajectory  electron  distributions.  The  use  of  the 
general  solution  to  bimodal  diffusion  theory  yields  a  naturally  unique  off-trajectory  projec¬ 
tion  within  the  limitations  of  the  assumed  physics.  Whether  the  assumed  physics  correctly 
describe  the  actual  processes  prevailing  in  the  particular  region  of  the  outer  electron  belt 
is  a  sepuate  issue  to  be  discussed  below.  Within  the  scope  of  diffusion  theory,  if  the 
on-trajectory  representation  is  exact,  the  projection  to  any  off-trajectory  point  is  also  in¬ 
dependent  of  path  because  the  solution  is  uniquely  specified  everywhere  that  it  is  valid. 
Since  the  on-trajectory  representation  is  not  given  a  priori  (but  given  as  a  fit  based  on 
a  data  set  that  mixeb  i  and  t),  an  important  question  arises:  Can  the  fit  errors  due  to 
adiasing  X  and  i  destroy  the  uniqueness  of  the  projection  by  making  it  path  dependent? 

We  have  investigated  this  problem  of  projection  consistency  (or  path  dependence) 
for  our  representation  by  projecting  to  an  off-trajectory  point  from  two  well-separated 
on-trajectory  points:  one  having  the  same  X  and  the  other  having  the  same  t  as  the  off- 
trajectory  point.  The  projection  paths  are  thus  along  the  t  and  the  L  axes,  respectively. 
Figure  10  shows  the  results  of  such  a  uniqueness  test  for  Day  110,  1979.  The  projection 
in  our  representation  is  indeed  path  independent.  This  test  clears  a  hurdle  for  using  the 
bimodal  diffusion  solution  for  outer-belt  dynamic  modeling,  since  it  demonstrates  that 
uniqueness  is  not  destroyed  by  the  fit  errors.  It  must  be  noted  that  this  is  only  an  internal 
consistency  test,  not  a  vzdidity  test  that  diffusion  theory  correctly  describes  the  prevailing 
physics  of  the  outer  belt.  A  validity  test  requires  data  sets  from  more  than  one  satellite. 
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Three-dimensional  presentation  of  the  phase  space  density  of  (a)  the  data  base 
and  (b)  its  reconstitutions  for  the  28S.5  ke\’  channel  of  Day  165,  1980  in  pitch 
angle  and  L.  Large  errors  are  found  primarily  in  the  corner  of  (b)  pertaining  to 
large  L  and  large  pitch  angles.  Noisy  data  near  L  —  5.7  also  resulted  in  some 
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Figure  6.  Three- dimensionad  presentation  of  the  phese  space  density  of  (a)  the  data  base 
and  (b)  its  reconstitutions  for  the  448  keV  channel  of  Day  110,  1979  in  pitch 
angle  and  L.  As  in  Figure  5.  the  large  errors  are  mainly  confined  to  the  corner 
of  (b)  pertaining  to  large  L  and  large  pitch  angles. 
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Figure  7.  The  distribution  of  fit  percentage  errors  in  the  reconstitution  of  all  the  data 
points  in  the  Day  165,  1980  data  base.  The  error  distributions  of  the  initial 
and  final  iterations  are  shown  as  dots  and  Eisterisks.  respectively. 
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Figure  8.  The  distribution  of  fit  percentage  errors  in  the  reconstitution  of  all  the  data 
points  in  the  Day  110,  1979  data  base.  The  error  distributions  of  the  initial 
and  final  iterations  are  shown  as  dots  and  asterisks,  respectively. 
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Figure  9.  Comparison  of  the  radial  diffusion  coefficients  determined  in  Paper  IV  for  the 
energy  range  0. 2-2.0  keV  (mottled  area)  with  Paper  I  (cross-hatched  area)  and 
with  the  compilation  of  West  et  al.  !l981j,  whose  individual  analysis  yields 
curve  1.  Curves  2-8  are  experimental:  2,Tomassian  et  al.  [1972];  3,  Farley 
1969;:  4.  Newkirk  and  W'alt  .'1968bi;  5,  Lanzerotti  et  al.  1970:;  6,  Newkirk  and 
Walt  ;1968bj,  Falthammar  [1970];  7.  Newkirk  and  W'alt  '1968a];  8,  Kavanaugh 
1968;.  Curves  9-14  are  theoretical  or  semi-theoretical:  9,  Nakada  and  Mead 
'1965::  10,  Tverskoy  [1965];  11,  Birmingham  1969]:  12,  13,  Cornwall  [1968];  14. 
Holzworth  and  Mozer  [1979];  [References  17  through  29,  respectively]. 
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Figure  10.  Demonstration  of  path  independence  in  projecting  the  electron  distributions 
from  on-trajectory  to  off-trajectory  points.  Panels  A  and  B  are  measured 
and  fitted  phase  space  distributions  at  the  trajectory  locations  shown  on  the 
SCATHA  trajectory  panel.  Using  the  eigenfunction  representation,  the  distri¬ 
butions  are  projected  along  the  i  and  L  axes  respectively  to  the  off-trajectory 
point  C.  The  comparison  of  the  two  projected  distributions  is  shown  on  Panel 
C. 
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A  valid  dynamic  model  should  be  able  to  “predict”  one  satellite  data  set  from  the  solution 
fit  to  the  other  satellite’s  data  set,  and  vice  versa. 

We  have  repeatedly  stated  from  the  outset  that  the  vzilidity  of  the  form  of  simultaneous 
bimodal  diffusion  theory  chosen  for  this  study  is  very  severely  limited,  particularly  in 
regions  of  the  outer  electron  belt  at  L  >  7.  The  first  limitation  is  that  Walt’s  invariant 
[Reference  9]  remains  so  only  within  a  limited  slice  of  L.  This  property  of  piecewise 
invariance  is  exploited  as  much  as  possible  in  our  application  in  the  region  5  <  X  <  7,  but 
there  must  be  constraints  to  the  exploitation.  In  the  stable  regions  of  the  outer  belt  at  X  < 
5.  we  do  not  see  any  reason  why  the  application  of  the  representation  will  be  different  from 
that  shown  here,  although  it  is  necessary  to  test  it  with  another  satellite  data  set  because 
SCATHA  does  not  cover  that  region.  For  X  >  6,  we  sometimes  encounter  difficulties  with 
“butterfly”  pitch-angle  distributions.  We  recognized  at  the  outset  that  an  important  factor 
limiting  the  validity  of  the  representation  is  the  effects  of  third  invariant  violations  that  may 
or  may  not  be  describable  in  terms  of  diffusion.  Such  effects  can  be  due  to  severe  magnetic 
or  electric  shell  splitting  or  to  encounters  with  the  magnetopause.  In  such  encounters, 
the  electron  distribution  on  a  given  X-shell  can  be  changed  catastrophically.  Since  such 
third-invariant  violating  effects  are  deemed  to  be  most  effective  at  the  equatorial  region,  it 
is  usually  assumed  that  their  signature  is  a  loss  to  the  drift  shell  of  equatorially  mirroring 
particles  (those  with  ~  90°  pitch  angle)  leading  to  the  so-called  “butterfly”  distributions 
[References  15  and  16’.  Such  an  expectation  seems  reasonable  as  far  as  the  pitch-angle 
distribution  is  concerned;  however,  the  transport  in  the  configuration  part  of  phase  space 
is  a  matter  that  has  not  been  investigated  beyond  the  expectation  that  the  third  invariant 
is  drastically  violated. 

Since  our  representation  is  not  expected  to  be  valid  for  describing  such  drastic  viola¬ 
tions  of  the  third  invariant,  the  inclusion  of  such  distribution  functions  in  our  data  base 
incurs  large  errors  in  the  representation.  These  constitute  a  part  of  the  large-error  tail  of 
the  error  distributions  shown  in  Figures  7  and  8.  The  successful  treatment  of  the  “butter¬ 
fly”  events  constitutes  the  remaining  hurdle  towards  dynamic  outer  electron  belt  modeling. 
The  solution  to  this  important  problem  is  beyond  the  scope  of  this  brief  report  but  we 
shall  address  some  salient  points  here.  An  initial  reaction  to  the  failure  of  the  simple  form 
of  radial  diffuion  used  is  to  assume  that  another  form  of  radi2il  diffusion,  perhaps  in  third 
invariant  coordinates,  may  be  a  better  description  of  the  physics  involved  in  “butterflies”. 
Such  may  indeed  be  the  case,  but  a  more  primitive  question  needs  to  be  considered  first: 
.4re  encounters  involving  drastic  violations  of  the  third  invariant  necessarily  diffusive  in 
the  configuration  part  of  phase  space? 

In  order  to  answer  this  and  related  questions,  we  have  initiated  computer  simulations 
of  10  MeV  electrons  encountering  a  realistic,  non-turbulent  paraboloidal  magnetopause 
[Reference  30  in  the  configuration  of  a  Tsyganenko  [Reference  3]  field  model.  The  mag¬ 
netic  field  magnitude  drops  by  30  nT  across  the  magnetopause  at  the  sub-solar  point. 
The  field  configuration  outside  merges  smoothly  into  a  uniform  southward  interplanetary 
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configuration.  The  important  feature  of  the  magnetopause  encounter  is  that  the  elec¬ 
tron  executes  trapped  bounce  motion  occasionally  on  inward  field  lines  and,  in  between 
bounces,  it  skips  along  the  equatorial  regions  sometimes  inward  and  sometimes  outward 
of  the  magnetopause  boundary  with  differing  amounts  of  residence  time,  dependent  on 
pitch  angle  and  gyration  phases;  i.e.,  there  is  little  evidence  of  scatter  perpendicular  to  the 
boundary.  Such  anisotropic  scatterings  in  encounters  with  the  magnetopause  can  hardly 
be  described  in  simple  terms  of  isotropic  diffusion,  in  invariant  coordinates  or  otherwise. 
The  encormter  effects  are  somewhat  more  isotropic  in  configuration  space  if  the  magne¬ 
topause  is  not  smooth,  but  the  anisotropy  will  remain  since  the  average  magnetic  field  is 
non-zero  both  inside  and  outside.  Clearly  the  equatorial  skippings  are  likely  to  lead  to 
enhanced  pitch-angle  diffusion;  however,  pitch-angle  diffusion  cannot  produce  a  “butter¬ 
fly”  population.  On  ti.e  other  hand,  skippings  along  the  equatorial  magnetopause,  which 
cuts  through  many  i-shells,  clearly  constitute  transport  in  L  and,  consequently,  viola¬ 
tion  of  the  third  invariant,  zis  expected.  The  loss  to  a  given  T-shell  of  electrons,  which 
skip  to  a  different  T-shell  along  the  equatorial  magnetopause,  will  clearly  produce  a  “but¬ 
terfly”  distribution.  The  surprising  feature  of  the  simulation  is  that  this  transport  in  L 
is  anisotropic  in  configuration  space  and  can  hardly  be  considered  diffusive.  Based  on 
computer  simulations  of  magnetopause  encounters,  the  “butterfly”  distribution  is  but  a 
superficial  signature  peculiar  to  the  process  in  which  violations  of  invariants  occur  with 
drastic  anisotropy  in  configuration  space.  Further,  the  channeling  of  electron  motion  pro¬ 
vided  by  the  sharp  magnetic  field  gradient  at  the  magnetopause  boundary  surface  indicates 
that  it  acts  macroscopically  as  a  scattering  boundary  to  outer-belt  electrons,  having  simul¬ 
taneously  the  properties  of  a  sink  to  electrons  approaching  the  boundary  on  given  initial 
X-shells  intercepting  the  magnetopause  and  of  a  highly  anisotropic  source  of  electrons  to 
X-shells  intercepting  the  magnetopause  inward  of  the  initial  X-shell.  In  the  non-turbulent 
magnetopause  case,  virtually  no  electrons  were  lost  to  interplanetary  space.  It  appears 
that  diffusion  in  invariant  space  is  not  the  panacea  that  cursory  considerations  make  it  out 
to  be. 

By  testing  against  a  SCATHA  electron  data  b^lse,  we  have  found  a  computationally 
practicable  representation  for  the  outer-belt  electron  distribution  from  the  general  solution 
of  a  simple  version  of  the  simultaneous  bimodal  diffusion  theory  that  satisfies  the  basic  re¬ 
quirements  for  outer-belt  dynamic  modeling  inwards  of  X  ~  7.  The  representation  fadls  to 
characterize  the  small  portion  of  outer-belt  electron  populations  known  as  “butterfly”  dis¬ 
tributions.  The  nature  of  outer-radiation-belt  electron  encounters  with  the  magnetopause 
needs  to  be  investigated  beyond  the  superficial  formation  of  “butterfly”  pitch-angle  dis¬ 
tributions.  In  particular,  the  distribution  of  “butterflies”  in  configuration  space  can  shed 
light  on  the  scattering  processes  in  encounters  with  the  magnetopause. 
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VII.  RESULTS  AND  CONCLUSIONS 


Early  in  this  Space  Radiation  Test  Model  Study  we  developed  the  necessary  tools  for 
handling  the  SCATHA  SC3  electron  data  in  both  configuration  and  invariant  coordinates. 
We  mapped  the  phase-space  distribution  functions  in  terms  of  observable  parameters  at 
fixed  first  and  second  adiabatic  invariants.  We  have  subsequently  developed  the  inverse 
algorithms  for  mapping  the  distribution  in  terms  of  the  invariants  at  constant  energy  and 
pitch  angle,  the  observable  parameters.  A  computationally  efficient  procedure  for  calculat¬ 
ing  the  third  adiabatic  invariant  and  the  drift  motion  of  particles  was  also  demonstrated. 
Using  these  tools  the  SCATHA  environmental  data  base  was  surveyed  for  periods  of  inter¬ 
est  to  this  study.  We  then  proceeded  in  an  iterative  fashion  to  solve  the  bimodal  diffusion 
equation  and  test  it  with  the  data  base.  The  following  results  and  conclusions  can  be 
drawn  from  this  study: 

1.  A  computationally  efficient  representation  was  obtained  based  on  a  solution  of  the 
simultaneous  bimodal  (radial  and  pitch-angle)  diffusion  equation  in  the  configuration 
and  pitch-angle  coordinates  {L,x,t). 

2.  The  solution  representation  was  tested  against  the  SCATHA  SC3  data  base  and  was 
found  to  be  good  for  L  <  7. 

3.  The  Bessel  function  expansion  in  the  pitch-angle  coordinate  provides  a  means  for 
identifying  the  pitch-angle  distribution,  e.g.,  for  a  normal,  “butterfly”  or  isotropic 
distribution,  and  thus  for  monitoring  the  state  of  the  radiation  belt  over  time. 

4.  Data  obtained  on-trajectory  as  a  function  of  L{t)  and  i  do  not  necessarily  lead  to  a 
unique  off- trajectory  projection.  However,  we  have  explicitly  tested  our  representation 
and  procedure,  finding  that  they  provide  a  unique  projection  away  from  the  spacecraft 
orbit. 

5.  Based  on  the  off-trajectory  uniqueness,  our  representation  wiD  provide  dynamic  mod¬ 
eling  for  L  <  7,  if  there  are  no  “butterfly”  pitch-angle  distributions  present  in  the 
data  base. 

6.  “Butterfly”  distributions  lead  to  large  errors  when  they  are  included  in  the  data  set 
to  which  the  simultaneous  bimodal  diffusion  solution  is  applied. 

7.  The  nature  of  “butterfly”  distributions  has  been  elucidated  with  a  particle  simula¬ 
tion  code.  They  are  characterized  by  third  invariant  violation,  as  well  as  anisotropic 
transport,  in  configuration  space.  This  effect  has  not  been  discovered  until  now. 
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VIII.  RECOMMENDATIONS 

From  the  results  and  conclusions  in  the  previous  section  of  the  report,  the  following 

recommendations  are  suggested  for  implementation  in  the  near  future  in  order  to  perform 

dynamic  modeling  of  the  radiation  belts  in  the  CRRES  time  frame. 

1 .  In  order  to  implement  dynamic  modeling  of  the  outer  radiation  belt,  our  representation 
of  the  electron  data  should  be  used  for  X  <  6. 

2.  The  “butterfly”  distributions  should  be  investigated  further  with  computer  simulation 
and  the  findings  should  be  incorporated  into  the  diffusion  theory. 

3.  Shell-splitting  effects  on  lower-energy  particles  due  to  steady  electric  fields  and  mag¬ 
netic  fields  should  be  investigated. 

4.  A  procedure  should  be  developed  to  couple  all  of  these  effects  into  the  diffusion  theory 
and  our  data  representation,  leading  to  the  capability  for  dynamic  modeling  of  the 
radiation  belt. 

5.  The  dynamic  model  should  be  tested  with  SCATHA  SC3  data  for  the  outer  belt 
region. 

6.  The  dynamic  model  should  be  tested  with  CRRES  high-energy  electron  data  over  the 
whole  radiation  belt. 
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Simultaneous  Radial  and  Pitch  Angle  Diffusion 
in  the  Outer  Electron  Radiation  Belt 

Y  T  CuR',  R.  W.  Nightingale,  and  M.  A.  Rinaldi 


Lfickheed  Pah  AUo  Re\earfh  Laharutorv.  Cahtorma 


A  ^('lulKin  of  Ihc  hrmodal  (radial  and  pilch  anglel  diffusion  equation  for  ihe  radiation  bells  is  devel¬ 
oped  vsith  special  regard  for  Ihe  requiremenis  of  saielliie  radiation  bell  dala  analysis.  In  this  paper,  we 
use  this  solution  to  test  the  bimodal  theory  of  outer  electron  bell  diffusion  by  confronting  ii  with  satellite 
data  Satellite  vvbservations.  usually  over  finite  volumes  of  (/,.  f|  space,  arc  seldom  sufTicienl  in  space-time 
duration  to  cover  Ihe  relaxation  lo  equilibrium  of  the  enure  radiation  bell  Since  lime  scales  of  continu¬ 
ous  data  coverage  are  oflen  comparable  lo  that  of  radiation  bell  disturbances,  it  is  therefore  inappro¬ 
priate  lo  apply  impulsive  semi-mtinile  lime  response  si'luiions  of  dilTusion  theory  lo  interpret  dala  from  a 
hniie  window  of  IL.  r)  space  Observational  limitations  indicate  that  appropriate  solutions  for  the 
iiiicrprctaiion  of  saicIliie  dala  are  general  solutions  for  a  finite-volume  boundary  value  problem  in 
bimodal  dilliision  Here  we  lest  such  a  solution  as  the  prime  candidate  for  comprehensive  radiation  bell 
dynamic  modeling  by  applying  Ihe  solution  and  developing  a  method  of  analysis  lo  radiation  bell 
electron  dala  obtained  by  the  S(  A  I  H  A  satellite  at  moderate  geomagnetic  activity  The  results  and  Ihc 
generality  of  our  solution  indicate  its  promise  as  a  neix  approach  lo  dynamic  modeling  of  the  radiation 
bells 


I  IsrROUUfTloN 

In  the  three  decades  since  the  discovery  of  the  radiation 
bells,  the  physics  of  geomagnetically  trapped  radiation  has 
evolved  from  the  position  of  being  the  motivating  force  of 
sp.ice  plasma  physics  into  the  current  state  of  relative  quies¬ 
cence  This  IS  no!  necessarily  due  to  a  matured  understanding 
of  radiation  belt  dynamics:  rather,  it  is  due  to  a  relative  lack 
of  opportunities  lo  place  into  suitable  orbil  advanced  instru¬ 
ments  that  can  cover  ihe  large  dynamic  range  of  radiation  belt 
particle  variations  in  space,  time,  energy  and  pitch  angle 
(phase  .(luvei  New  discoveries  in  auroral  phenomena  in  the 
late  1970s  and  early  l9S0s  h.ive  also  allracled  research  interest 
away  from  the  radiation  belts  toward  the  low-energy  high- 
latilude  plasma  Indeed,  aside  from  geo-synchronous-orbit  in¬ 
struments  which  are  especially  suitable  for  studies  of  sources 
feg  .  Raker  er  a/..  I97f<],  (here  exist  few  electron  radiation  belt 
dala  seis  since  in  19bh  [eg.  IV  e.st  er  al .  1981,  and 

references  therein  |  that  can  provide  (he  appropriate  dynamic 
range  in  phase  space  coverage  suitable  for  specifying  the  dy¬ 
namics  of  the  outer  electron  radiation  belt  Among  these,  (he 
data  from  ihe  Lockheed  SL'-3  instrument  on  board  the 
SCATHA  lP78-2i  spacecraft  launched  in  1979  is  unique  in 
ih.ii  It  has  excellent  pitch-angle  resolution  covering  the  very 
n.iriow  losv  cone  .il  llic  qiiasi-cqii.itori.il  orbit  of  the  space¬ 
craft  [Rrai/aa  <■(  al .  198IJ  f  urthcr.  since  the  spacecraft  is  not 
exactly  geosynchronous,  the  Z.  shell  coverage  of  5  3  lo  8.7  in 
the  data  provides  a  testing  ground  for  further  specifications  of 
Ihe  dynamics  of  the  outer  radiation  bell 

Relative  inactivity  in  Ihe  scientific  aspects  of  radiation  bell 
physics  in  the  last  decade  does  not  imply,  however,  a  decrease 
of  applicalional  demand  for  radiation  environment  specifi¬ 
cation  Indeed,  as  a  result  ol  increasingly  complex  and  sensi¬ 
tive  space  insirument.ition  and  .iciivities.  demands  for  im¬ 
proved  radiation  environment  specification  have  become  Ihc 
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driving  force  for  a  revisit  and  a  further  cflcirt  al  scvlvmg  the 
important  problems  of  radiation  belt  dynamics  left  unsolved 
in  the  early  1970s  The  upcoming  Air  Force  NAS/s  Combined 
Release  and  Radiation  EITects  Satellite  (CRRES)  crystalirtcs 
these  research  trends,  for  the  goal  of  its  mam  complement  of 
instruments  is  to  develop  a  dala  base  for  comprehensive  speci- 
lieaiion  of  the  particle  populations  in  Ihe  space  radiation  cnvi- 
tonmcnl. 

The  collection  of  data  by  a  single  satellite  mission,  however, 
comprehensive,  is  not  sufficient  to  specify  the  worldwide  radi¬ 
ation  enviropmrn’  in  space  and  time  This  can  only  be 
achieved  by  Ihe  Iheorelical  explication  of  the  dy  namics  with  a 
priori  physical  laws  so  that  piedictivc  models,  based  upon  and 
verilied  by  dala.  can  be  constructed  for  applicalional  needs 
Such  an  cbjcclive  is  the  eventual  goal  for  which  the  present 
work  IS  intended  as  a  first  step 

In  this  paper,  we  focus  on  deriving  a  practicable  solution  of 
Ihc  diffusion  physics  ft'r  outer  radiation  bell  electrons  The 
solution  must  be  appropriate  for  the  finite  space-lime  cover¬ 
age  of  Ihc  dala  so  as  to  limit  theoretical  assumptions  that  arc 
not  directly  relevant  lo  Ihe  data  We  shall  dcmtvnsiralc  that 
these  goals  can  be  achieved,  al  least  for  the  region  of  space 
covered  by  the  SC  ATHA  spacecraft  in  a  moderately  active 
outer  radiation  bell,  by  .ipplic.ition  of  a  general  eigeiifuiiction 
expansion  solution  to  the  simultaneous  radial  and  pitch  angle 
diffusion  equation  Although  our  solutions  are  admittedly  ap 
proximations  lo  the  complex  bimodal  dd'asion  physics,  our 
work  here  represents,  insofar  ,is  we  can  determine,  a  first  com¬ 
prehensive  allempi  to  confront  simultaneous  bimodal  ditfu- 
sion  theory  with  data  As  such  ,ii  the  curreni  primitive  si.ige 
of  Simultaneous  bimodal  diffusion  theory,  we  lound  il  neccs- 
sarv  to  m.ike  a  number  of  simplify  ing  assumptions  m  order  lo 
achieve  a  complete  solulii'n,  which,  insofar  as  we  c.in  detci 
mine,  has  not  been  wrillcn  down  before  We  impose  siikI 
discipline  in  filling  Ihc  dala  lo  Ihe  solution  however  Ihe 
result  ts  a  til  lo  the  entire  d.ila  in  eiieigv  s|vecliiiiii  /  sluTI 
pilch  angle  anti  lime  1  he  dillusum  parameters  viclded  bv  the 
III  compaie  well  with  other  dctermirialions  from  single mode 
diffusion  theory  II  is  impoit.mi  lo  iioie,  however,  ih.it  the 
icsults  of  the  present  woik  .ire  not  me, ml  lo  be  applic.d'le  lo 
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Fit:  I  Ten-minutc  averugcd  eleciron  fluxes  inlegraied  over  all 
pitch  angles  for  seven  energy  channels  measured  by  the  Lockheed 
SC'l  insirumeni  on  board  SCATHA  during  day  165,  1978  This  is  the 
data  set  to  be  analyzed  in  Ihis  paper 


outer  bell  electrons  under  very  active  geomagnetic  conditions 
when  shell  splitting  and  niiigiielopause  transit  efiects  become 
dominant. 

In  section  2,  the  unique  characteristics  of  the  Lockheed 
SCATHA  data  base  are  discussed  in  terms  of  their  applica¬ 
bility  to  explicating  dynamics  of  the  outer  radiation  bell.  In 
section  3.  a  brief  review  of  the  theoretical  development  in 
radiation  belt  diffusion  physics  is  followed  by  the  derivation  of 
a  general  btmodal  simultaneous  diffusion  solution  suitable  for 
situations  encountered  in  satellite  sampling  of  radiation  belt 
dynamics.  The  most  important  point  is  that  satellite  observa¬ 
tions  in  a  linite  segment  of  space-time  are  not  amenable  to 
intcrprctaiions  in  terms  of  a  relaxation  to  equilibrium  The 
application  of  theory  upon  the  chosen  data  set  is  presented  in 
section  4  In  section  5.  liic  physical  implication  of  the  data 
representation  and  its  applications  to  comprehensive  dynamic 
modeling  of  the  outer  radiation  belt  will  be  discussed 

2  Daia  Hasi  Chakacii  RISIIf  s 

The  IJSAf'NASA  l’78-2  spacecraft  studying  satellite 
charging  at  high  altitudes  (SCATHAl  was  launched  on  Janu¬ 
ary  30.  1979.  In  its  final  orbit  the  satellite  nominally  spins  at 
~  I  rpm  and  traverses  an  altitude  range  between  ~  27,500  km 
at  perigee  and  -45.2(X)  km  at  apogee  twice  a  day  within  a 
latitude  range  of  +  15  .  This  coverage  corresponds  to  an  L 
shell  range  of  about  5  3  to  8  7  in  the  near-geosynchronous 
region  Keganlless  of  the  orientation  of  the  orbit  with  respect 
to  the  F.irth-Suii  line.  SCATHA  spends  ~83''.,  of  the  time 
beyond  I.  -  b 

The  Lockheed  high-energy  panicle  spectrometer.  SC3.  is  a 
solid-state  particle  telescope,  employing  a  scries  of  surface  bar¬ 
rier  silicon  detectors,  that  primarily  measures  energetic  elec¬ 
tron  fluxes  as  a  function  of  energy  and  pitch  angle  The  instru¬ 
ment  normally  operates  in  two  energy  modes  with  the  lower 
mode  covering  the  electron  energy  range  between  47  to  299 
keV  in  12  channels,  each  about  21  keV  wule  The  higher 
cnergv  mmlc  extends  over  the  eleciron  energy  range  from  256 
lo  4'l7i)  keV  in  12  channels  each  393  keV  wide  SC3  has 
excellent  pilch  angle  resolution  and  pointing  alignment,  suf- 
hcicnl  to  look  into  the  loss  cone,  thereby  .illowing  precise 
lining  of  the  pilch  angle  disirihuiioiis  I  he  lield  of  view  of  the 
deieclor  is  3  .  providing  an  elledive  angul.ii  width  for  the 


detector  response  of  ~4  ,  which  includes  the  3  sweep  of  the 
delector  for  each  0.5-s  accumulation  period  This  detector  re¬ 
sponse  should  be  compared  with  the  expected  loss-cone  width 
of  4'  to  6'  in  the  near-synchronous  region  The  spectrometer 
is  aligned  perpendicular  to  the  spin  axis  of  the  satellite  sweep¬ 
ing  through  almost  all  pitch  angles  (from  0  to  9U  I  four  limes 
per  minute.  The  SC3  instrument  has  been  acquiring  data  and 
operating  normally  since  launch  A  more  complete  description 
of  the  instrument  can  be  found  in  Reagan  et  at  [1981] 

The  Lockheed  SCATHA  SC3  data  base  consists  of  ~9.50 
days  of  digitized  data  on  magnetic  tape  from  launch  in  1979 
to  the  present  in  several  intervals.  The  pitch  angle  for  each 
accumulation  period  is  determined  using  data  from  the  o.n- 
board  vector  magnetometer  (SCIl).  From  the  data  tapes  a 
condensed  data  base  of  10-min  averaged,  pitch  angle  binned 
electron  fluxes  for  the  24  energy  channels  over  the  full  energy 
range  of  47  to  4970  keV,  has  been  generated  for  on-line  stor¬ 
age  on  a  VAX  11/780.  For  each  10-min  lime  interval  the  data 
are  grouped  by  pitch  angle  values  into  3°  bins  from  0  to  90  . 
making  30  pitch  angle  bins  for  each  energy  channel  in  the 
interval.  Every  10-min  interval  also  includes  the  universal  time 
(UT),  local  time  (LT),  latitude,  longitude.  Kp.  Dtt.  magnetic 
latitude,  B/B^,  value,  and  the  L  shell  value  for  the  SCATH  A 
spacecraft.  The  quiet  time  Olson-Pfiizer  [O/sort  ei  al .  1979] 
magnetic  field  model  together  with  the  formulations  of  A/i  //• 
wain  f!961]  for  L  shell  were  used  lo  calculate  the  last  three 
ephemerides.  This  mulliparameier  data  base  facilitates  the 
generation  of  energetic  eleciron  distribution  functions  in 
adiabatic-invariant  coordinate  space 

The  SCATHA  orbit  has  a  daily  perigee  of  L  ~  5.3  and  an 
apogee  that  often  carries  it  beyond  L  =  8.6.  The  orbit  pre- 
cesses  eastward  at  ~5.3'  of  longitude  per  day;  so  the  line  of 
apsides  gradually  rotates  eastward  with  respect  to  the  noon- 
midnight  meridian  at  -4  3'  per  day  An  especially  important 
region  is  near  7  R^,  where  irappcd-parlicle  drift  orbits  can 
cross  the  magnetopause  boundary,  the  latter  depending  on  the 
activity  level  and  located  somewhat  beyond  L  =  7  on  the  day- 
side.  Electrons  with  large  pitch  angles  that  encounter  the  mag¬ 
netopause  boundary  are  presumed  lo  be  scattered  and  ef¬ 
fectively  removed  from  the  drift  shell  on  which  they  had  start¬ 
ed  As  a  result  of  losses  al  the  magnetopause,  the  electron 
distributions  on  the  evening  side  are  often  depleted  at  large 
pitch  angles;  this  depletion  is  very  pronounced  when  the  mag- 
entic  field  is  distorted  in  the  development  of  a  subsiorm  (or 
storm),  and  is  an  important  indicator  of  geomagnetic  activitv 
Siheck  et  al  [1987]  discussed  such  a  case  of  anomalous  shell 
splitting  but  our  distributions  in  the  selected  data  base  do  not 
show  such  features 

The  SCATHA  SC3  temporal  data  set  was  from  the  laiicr 
portion  of  day  165,  June  13,  1980,  a  moderately  active  period 
with  Kp  =  4-4  and  Dsi  =  -.50  Two  days  earlier  a  geomag¬ 
netic  storm  had  occurred,  but  the  magnetosphere  had  sirue 
relaxed  enough  that  the  observed  electron  fluxes  on  sl.iv  If'S 
were  smoothly  varying  with  no  pitch  angle  signatures  of  p.oii- 
cle  loss  in  the  magnetopause  Figure  I  shows  the  In  min 
averaged  data  integrated  over  all  pitch  angles  for  seven 
channels  of  electron  flux  IJuring  this  period  the  sp,Kc..rafi 
was  drifting  inward  from  /.  --  7  2  lo  /.  =  5  3  I  he  /  shell  i.inge 
included  in  the  analysis  has  been  restricted  lo  I.  --  t- 48  to 
eliminate  the  large  variability  of  the  eleciron  fluxes  bev'"d  the 
geosynchronous  region  fit  esi  ci  al .  1981  ]  and  the  low  ^ount 
rales  due  It'  the  approximately  exponential  fall  oil  wuh  in 
creasing  I.  shell  [Reuivuri  et  al .  |98|]  In  addition  we  have 
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Aelecled  for  the  formal  analNsis  only  ihose  data  channels 
above  25()  keV  so  as  to  eliminate  channels  which  are  most 
likely  to  be  aflecled  by  convective  forces  or  resulting  from 
consequences  of  the  major  storm  preceding  the  data  period 
Data  with  electron  energies  above  22US  keV  have  low  statis¬ 
tics  and  are  also  excluded  from  the  data  set  Data  from  energy 
channels  outside  of  these  limits  are  also  analyzed  to  show  the 
extraneous  factors  in  data  analysis  that  a  theoretical  diffusion 
picture  must  contend  with  Figure  2  illustrates  typical  pilch 
angle  distributions  during  this  lime  for  the  selected  energies. 
.Although  the  period  was  moderately  active,  the  magnetic  field 
was  smoothly  varying,  with  a  field  direction  that  remained 
perpendicular  to  the  spacecraft  spin  vector.  This  allowed  SC3 
to  look  into  the  loss  cone  for  the  whole  period,  providing  data 
over  the  full  pilch  angle  range  In  addition,  the  pilch  angle 
distributions  were  smoothly  varying  and  without  signatures 
that  indicate  loss  at  the  magnetopause 

I  in.illv.  It  nitisi  be  reiterated  that  the  above  test  data  set 
w,is  preselected  on  the  bases  of  moderate  geomagnetic  ac- 
tivitv.  nonexistence  ol  magnetopause  boundary  crossing  signa¬ 
tures,  avoidance  of  convective  etlects  and  statistical  sufTiciency 
.iltme  '  'nee  the  data  set  was  chosen,  no  other  data  sets  were 
analyzed  until  the  tests  presented  below  were  cotnpicted.  In 
the  course  I'f  the  work,  no  other  criterion  was  applied  to  the 
data  and  no  search  for  another  data  set  was  conducted. 

1  HivtiiiiAi,  Dirnsitiv  is  iiii  OciiiK  Bili 

The  framework  of  radiation  bell  difTusion  theory  has  been 
based  on  the  well-developed  Jacobi-Hamillon  formulation  for 
motion  of  magneticallv  trapped  particles  and  the  Boltzmann 
equation  lor  an  ensemble  of  such  particles  in  space  [cf.  Schul: 
(itiJ  /.iiii;i r-,uii.  I4'4)  In  the  ideal  case,  because  of  the  exis¬ 
tence  of  cyclic  co-('rdinales  ca'rrespondiag  to  the  gyration, 
bounce  and  drift  periodic  motions,  the  slate  of  the  radiation 
bell  IS  given  by  the  srilution  to  the  collisionless  Boltzmann’s 
equation  fi'r  the  distribution  function 

/(  \f,  ./.  d'.i  =  const  (I) 

wliere  ( '/  J.  ‘I’>  are  respectively,  the  (irsi,  second  and  third 
adiabatic  iiivarianls  (/|.  I /,!  to  ihe  motion  To  the  extent 
that  scallcrings  of  Ihe  trapped  particles  thy  plasma  waves  m 
the  magnetosphere  and  by  interactions  with  the  collisional 
.itrnosphere  at  the  ntirrorsi  violate  the  invariants,  the  distri¬ 
bution  /  undergoes  redistribution  in  (M.  J.  d'l  space,  or  equiv¬ 
alently.  in  phase  space  This  redistribution  of  /  is  most  rigor- 
ouslv  treated  in  kinetic  space  plasma  physics  as  solutions  to 
ftoli/iii.iiiii  s  equation  with  collision  integrals  However,  a  tra¬ 
dition. il  F  okkcr  i’ianck  development  of  microscopic  kinetic 
iheorv  Ic.ids  to  .i  humulation  of  the  familiar  diffusion  picture 
for  which  Ihe  elleds  of  the  above  scallcrings  m  random-phase 
.igcrcgale  are  assumed  to  be  described  by  Ihe  diffusion  equa¬ 
tion 


whore  we  h.ivc  .idopicd  the  convention  that  repeated  indices 
.ire  Slimmed  .wet  'he  canonical  iiuari.int  variables  In  (2),  the 
eiicigies  ol  i'k  imiIilIc  population  arc  assumed  sufTicienllv 
Inch  ih.il  ihi  o' vsics  of  the  scallcrings  is  contained  m  the 
eomponenis  oi  ilie  diiliision  matrix  0,,  The  scallcrings  that 
violate  Ihe  .uh.ib.iiic  in-.inanls  e.in.  in  cenct  it,  violate  anv 
eornbin.Uioii  .-I  them  .it  the  s.imc  lime  and  the  existence  of 


I  ig  2  Typical  pilch  distributions  for  Ihe  seven  energy  channels  of 
the  data  set  referred  to  in  Figure  I 


many  plasma  wave  modes  imply  that  the  nine  terms  on  Ihe 
right-hand  side  of  (2)  may  in  general  be  hopelessly  tangled  As 
a  result,  il  can  almost  be  said  that  radiation  bell  diffusion 
Ihcivry  consists  of  attempts  to  consider  approximations  in  iso¬ 
lating  Ihe  dominant  contributions  to  the  right-hand  side  of  (2| 
Included  in  such  efforts  are  alicmpis  to  discover  coordinates 
other  than  invariants  in  which  only  one  or  two  dominant 
diffusion  terms  and  diagonal  diflusion  matrices  can  be  ob¬ 
tained.  linfortunately.  only  one-dimensional  considerations,  in 
radial  or  in  pitch  angle  difTusion,  have  previously  been  at¬ 
tempted 

The  intractability  of  (2)  in  any  form  more  complex  than  that 
of  a  single  mode  has  given  rise  to  the  traditional  treatment  of 
single-mode  radiation  belt  diffusion  in  some  single  variable, 
say  y.  in  the  form 


(T 


1  0  «V  / 

-  -  KD,,  —  -  -  +  i' 
K  ly  <y  r 


t.f) 


where  K  is  a  Jacobian  sh.ppropriaie  for  the  transformation 
from  particular  canonical  terms  in  |2|,  /  r  is  a  loss  term  and  S 
IS  a  source  term  Both  loss  and  source  terms  arc  ad  hoc  repre- 
scnlalions  under  which  the  effects  of  other  untreated  ditiusion 
terms  are  swept  F.laborale  solutions  of  initial  value  problems 
for  impulse  responses  in  the  semi-inhnilc  lime  domain  for 
single-mode  radial  diffusion  in  Ihe  ad  hoc  form  ofl.fl.  m  whicli 
V  -  L.  have  been  evolved  [eg.  Schul:.  1986;  Schuh  unit 
\ewman.  1987].  Single-mode  pilch  angle  diflusion  has  also 
been  treated  in  detail  teg,  see  references  in  St  huh  tinil  I  iin 
•III  (  1974,  chap  2]|.  .Such  elaborate  solutions  to  an  ad  hoc 
onc-dimensional  approximation  are  clearly  of  value  in  deline¬ 
ating  the  theoretical  behavior  of  radiation  belt  dilfusion  in  the 
single  chosen  variable  i  in  the  semi-infinile  lime  domain,  but 
progress  in  developing  Ihe  multimode  area  of  radiation  bell 
ditlusion  Iheorv  is  also  called  for  Let  us  now  consider  Ihe  role 
ol  data  in  the  next  step  in  complexity,  i  e  ,  m  multimode  dillu 
Sion  in  finite  space-time  II  is  crucial  to  note  that  this  new  sup 
calls  for  two  new  elements  that  will  require  special  consider 
alion  III  nuillimode  dillusion.  and  (2|  finite  space  time  sp.m 
Since  our  general  solution  is  developed  to  address  llicse  two 
issues,  a  discussion  of  each  is  necessary  at  Ihe  outset 

Bimodal  diflusion  Isimultaneous  in  radial  and  pitch  anglel 
forimilalions  were  discussed  early  in  the  dcveli'pment  ol  ladi 
Jtion  bell  phvsics  [lim  ri  n.lcl.  I96K.  ']  hciuiiiriitis.  |96S.  I  rank 
cnihni  1 1  ill.  1U68.  /  ufl/ijnnnnr  iiiut  Wall.  1969,  Wall.  197(1] 
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lli'U'  llic  Lliicf  diiriciilly  kkas  in  developing  ii  scheme  lo  urt- 
condilion.illy  diagonuli/e  llie  dillusion  operator  so  that  solu¬ 
tions  hy  the  separation  of  independent  variables  could  be  used 
fv>r  the  dillusion  equation.  H'all  [1970]  introduced  a  variable 


C  =  mbjb. 


(4) 


which  IS  shown  to  be  approximately  conserved  by  both  radial 
and  pitch  angle  dilTusion.  In  (4).  the  magnetic  field  (S|  sub¬ 
scripts.  Ill  and  c.  refer  lo  the  mirror  and  the  equator,  respec¬ 
tively.  Under  the  approximate  invariance  of  C  the  variables  or 
and  L  act  as  approximate  canonical  coordinates  so  that  the 
dillusion  equation  can  approximately  be  written  in  diagonal 
form  [f  ullliaiiiiimr  and  H'all,  1969;  Roederer  and  Schulz,  1969] 


-V 

rr 


I  r 
vT'tvirv 


7(.x)D„ 


L. 


2  ♦  (3 


(51 


where  v  =  cos  y.  and  T(.v)  is  the  quarter-bounce  integral, 
which  can  be  conveniently  approximated  by  expansion  in  sev¬ 
eral  slightly  dilVercnl  forms  [e  g..  Schulz  and  Lanzeroui,  1974, 
p  19;  Oiiriilvoii.  1976]  Here  we  use  Davidson's  formula 


7(v)  -  I  .1802  -  0.6.197(1  -  *  (6) 


In  {.S).  a  source  term  Six.  L)  (or  a  loss  term)  is  not  needed  to 
represent  ellccls  of  other  modes  of  dillusion,  as  is  customary  in 
singic-modc  dillusion  treatments,  unless  a  real  source  (such  as 
neutron  decay)  exists  in  the  volume  under  consideration. 
Indeed,  if  pitch  angle  (x)  and  radial  (L  shell)  diffusions  repre¬ 
sent  all  the  particle  transport  in  the  solution  volume,  the 
boundary  conditions  for  /  specified  by  satellite  observations  in 
a  given  finite  space-time  volume  are  the  "sources"  of  (5).  This 
point  will  be  discussed  more  thoroughly  below  Because  {  in 
(4)  IS  not  an  exact  invariant  corresponding  to  which  the  x  and 
L  terms  of  (5)  are  exactly  separable,  the  index  v  in  (5)  has  a 
vestigial  dependence  on  x.  For  example,  v  =  J  at  x  =  0  for 
which  (.S)  is  exact  and  v  =  2  at  v  =  I  In  the  region 
II  <  \  <  (194,  1  i  1  IS  a  very  good  approximation.  In  the 
region  t)94  <  v  <  I.  there  are  cross-dilfusion,  i.e.,  /7,,,  terms 
in  (S)  proportional  to  (i  -  .1)-(l  -  v'')’.  These  are  small  and 
will  be  Ignored 

Aside  from  the  above  approximations  arising  from  repre¬ 
sentational  and  kinematical  considerations  m  the  reduction  of 
the  full  multimode  dillusion  equation,  the  bimodal  equation 
|5|.  similar  to  the  ad  hoc  single-mode  equation  (.1|.  is  also  beset 
with  Ihe  same  physical  uncertainties  in  the  proper  forms  for 
the  dillusion  coeiriciciils  Wilhoul  simiillaiicous  and  global 
measurements  of  wave  modes  and  their  amplitudes,  there  is 
virtually  no  information  on  the  nature  of  D,,  and  0,^  for  a 
given  radiation  bell  data  set  This  major  lack  of  supporting 
inlormalion  for  all  known  data  sets  implies  that  one  of  the 
"/croih  order"  tasks  in  radiation  belt  diffusion  theory  is  to 
dcicrniine  the  approximate  properties  of  the  coefficients.  In  ad 
hoc  single-mode  diffusion  theory,  it  is  often  assumed  that  a 
diffusion  parameter  can  be  derived  simply  by  spatial  quadra¬ 
ture  of  the  diffusion  equation  [eg.,  Schulz  and  Lanzcrniii. 
1974,  chap,  5).  provided  that  the  data  is  sulficicntly  compre¬ 
hensive  lo  allow  for  accurate  determination  of  derivatives  of 
disliibiilion  functions  and  that  only  a  single  mode  of  dillusion 
IS  I'pei.iiivc  on  the  data  set  However,  both  of  these  are  very 
rcsiriciive  assumplions  C/uadralure  techniques  arc  based  on 
the  assumption  Ih.il  salcllitc  observations  provide  HI..  I)  m 


such  a  way  that  derivatives  in  space  and  lime  can  be  per¬ 
formed  independently.  Except  by  ad  hoc  assumption  and  per¬ 
haps  for  geosynchronous  satellites  under  certain  geomagnetic- 
conditions,  no  satellite  data  can  satisfy  such  a  condition  As 
for  the  assumption  that  only  single-mode  diffusion  is  oper¬ 
ative,  it  is  easy  and  tempting  to  argue  that  the  time  scales  of 
radial  and  pitch  angle  diffusion  are  widely  separated  so  that, 
for  example,  magnetic  impulses  can  cause  violation  of  the 
third  invariant  without  disturbing  the  gyration  and  bounce  of 
the  particles.  Such  arguments  are  well  and  good  in  theory,  but 
the  fact  of  nature  is  that  electric  and  magnetic  fluctuations  of 
all  time  scales  have  a  tendency  to  occur  together.  A  lime  span 
in  which  magnetic  impulses  occur  is  also  a  likely  period  for 
increased  plasma  wave  activity;  therefore  it  is  perhaps  futile 
(or  somewhat  vexing)  to  force  what  occurs  naturally  (simulta¬ 
neous  multimode  diffusion)  into  too  simple  a  mechanism  even 
though  simplicity  allows  for  a  higher  degree  of  theoretical 
sophistication  in  the  method  of  analysis.  These  are  the  reasons 
that  we  choose  to  investigate  multi-mode  diffusion,  w  hich  is  of 
course  also  beset  by  its  own  problems. 

For  multimode  diffusion,  the  quadrature  technique  simply 
involves  too  many  partial  differentiations  of  data  to  be  of  any 
value.  Thus  to  obtain  dillusion  parameters  In  multimode  dilfu- 
sion,  the  only  available  technique  is  to  hi  the  data  to  solutions 
of  (5)  under  various  models  of  the  diffusion  coefficients  This, 
in  turn,  implies  that  a  measure  of  uncertainty  associated  with 
the  variations  of  diffusion  coefficients  of  (5)  is  necessarily 
mixed  with  small  residual  variations  associated  with  the  defi¬ 
nition  of  separable  variables  discussed  above.  These  variations 
and  uncertainties  can  only  be  handled  in  the  manner  of 
“learning  curves,"  i.e.,  by  sequential  improvements  built  upon 
a  “zeroth-order"  analysis  The  main  purpose  of  this  paper  is  to 
develop  this  "zeroth-order"  analysis  for  multimode  dilfusion 
As  such,  we  shall  make  simplifying  assumptions  on  (.‘'I  that  are 
approximate  when  considered  in  the  applied  mathematics 
viewpoint  but  will  clearly  demonstrate  the  physical  character 
of  the  solution.  Further,  a  comprehensive  solution  technique 
and  its  application  will  be  developed  to  test  the  reasonable¬ 
ness  of  the  multimode  dilfusion  paradigm 

Ihe  guiding  principle  of  the  following  development  is. 
therefore,  to  obtain  an  exact  solution  to  a  simplihed  version  of 
(5)  that  IS  easily  applied  to  analyze  data  without  changing  the 
physical  character  of  the  dillusion  equation  Noic  that  the 
emphasis  is  to  obtain  an  exact  solution  lo  an  approximate 
equation  rather  than  an  approximate  solution  lo  an  exact 
equation  because  we  have  demonstrated  above  that  (.‘'I.  as  in 
the  case  of  singlc-niodc  dilfusion  (.').  is  by  no  means  an  exact 
equation  7  hercforc.  lo  Ihc  exicnl  lhal  (.ll  is  approvim.ilc 
anyway,  further  simplifying  assumplions  ,ire  acccpi.ihle  as 
long  as  they  do  not  change  the  character  of  ihe  equation  Wc 
have  attempted  many  optional  simplilicalions  of  ifi.  bui  scl- 
tlcd  on  Ihe  following  for  boih  m, iihcm.ilic.il  ,ind  physic.il  rea¬ 
sons  (I)  Wc  assume  7(xl  is  const, ml  for  m.ilhemalic.il  con¬ 
venience  because  the  Ihei'ry  of  ;inalylic  conimualion  over  a 
nonsingular  domain  can  be  applied  lo  show  that  the  solutum 
character  is  not  changed  if  7ixl  is  in  the  neighborhov'd  of 
unity  and  nonconslani  No  singularity  of  Ihe  dillusion  equa¬ 
tion  IS  encc'unlcred  when  7(v|  is  varying  frv'm  a  consiani  in 
the  vicinity  of  unity  Ihe  impact  of  this  assumplion  on  the 
pilch  angle  eigenfunciions  is  very  small  as  .i  comp.insv'n  be¬ 
tween  Ihe  exact  cigenlunciion  of  ihe  pilch  angle  dilfusion  op- 
crali’t  inl'ilvsilh  variable  fivland  Bessel  functions  which  are 
Ihe  cigenliiiiclions  for  our  .ippiovim.ilion.  h.is  been  given  bv 


A-A 


Cnii  ET  al  ;  Simultaneous  Radial  and  Pitch  Anole  DirrusioN 


Sihul:  amt  Lmizcroiii  [  1974.  p  166]  Kurther,  the  Bessel  func¬ 
tions  form  a  conveniently  available  complete  set  which  can 
span  the  domain  1  >  v  >  0  aside  from  being  eigenfunctions  of 
the  approvimate  pitch  angle  diffusion  operator.  (2)  The  pitch 
angle  dilTusion  coefficient  D.,  is  assumed  to  be  independent  of 
1.  since  the  range  of  L  covered  in  the  data  set  is  fairly  small. 
The  energy  dependence  of  D,,  is  to  be  determined  by  the  data 
analysis  (3l  The  diffusion  equation  cannot  be  solved  unless  a 
form  IS  specified  for  the  radial  diffusion  coefficient.  This  is  true 
in  single-mode  radial  diffusion  also  The  radial  diffusion  coef¬ 
ficient  D,i  IS  assumed  to  be  of  the  general  form  [e  g.,  Schulz 
and  l.anzcraiu.  1974.  p.  89;  Theodoridis.  1968] 

0,1.  =  s‘0  (7> 

where  ;  and  are  parameters  to  be  determined  by  the  data 
analysis  In  general.  4  can  be  an  arbitrary  analytic  function  of 
V  and  energy  The  energy  dependence  is  determined  by  apply¬ 
ing  the  solution  to  each  energy  channel.  If  i  depends  on  x,  (5) 
IS  still  separable  but  a  set  of  parameters  corresponding  to 
matrix  elements  of  s  in  the  eigenfunctions  of  the  pitch  angle 
diffusion  operator  must  be  determined  in  addition;  a  special 
case  of  which  is  considered  by  H  all  [1970]  These  matrix 
elements  cross-couple  pitch  angle  modes  [Fdllhammar  and 
H  all.  1969]  Inclusion  of  cross-mode  coupling  is  a  compli¬ 
cation  but  IS  not  included  in  this  initial  analysis;  therefore,  for 
this  analysis,  we  assume  .;  to  be  independent  of  .x  to  simplify 
expressions  It  should  be  noted  that,  although  there  is  no 
cross-mode  coupling,  there  is  simultaneous  bimodal  diffusion. 
There  extsts  scarce  thcorettcal  guidance  on  the  magnitudes  of 
/).,  and  ;  in  the  outer  bell,  but  there  are  theoretical  reasons  to 
expect  that  /i  ~  6  -  Id  [e  g  .  Si  liul:  and  I.anzt’rolli,  1974.  chap. 
7.  H  I'si  et  111.  1981]  We  assume  p  =  10  in  this  work.  Under 
these  simplifying  assumptions,  the  diffusion  equation  can  be 
written 


Cl  1  ll 

r V *  X  r\ 


(  L 


(^n 

JU  *  2-  i.\  2»V 

eL] 


(8) 


where  L),,.  4.  /i  and  i  are  considered  to  be  constant  parame¬ 
ters 

This  simplilicd  equation  is  to  be  solved  in  a  small  finite 
domain  /.,  <  /.  <  /-„  expected  to  be  in  the  outer  bell  and  in 
the  pilch  angle  domain  0  <  x  <  x,.  where  x_  is  cosine  of  the 
loss  cone  angle  The  cllccl  of  pitch  angle  dilfusion  is  to  define 
pitch  angle  behavior  of  /  inside  and  outside  of  the  loss  cone 
boundary  x  =  x,{/.l  1  he  dependence  of  x,  on  1.  is  yet  again 
another  f.ictor  which  mixes  the  x  and  L  operators  in  (8|;  how¬ 
ever.  (or  the  outer  /one  region,  where  (/.„  -  /.,)  <  2.  is  very 
iie.irly  constant  and  unity  We  assume  that  x,  is  constant  in 
the  data  span  so  that  separable  solutions  can  be  obtained  for 
(8i  with  the  present  lest  case  This  is  an  important  point  when 
we  deal  with  data  from  a  large  domain  of  /,  (such  as  that  from 
CRRE.Sl  which  involves  large  variations  of  the  loss  cone  both 
from  the  large  variations  in  1.  and  in  distances  to  the  mirrors 
for  our  SC'ATHA  data  set  the  loss  cone  dependence  on  1.  is 
entirely  negligible 

In  general  many  different  kinds  of  solu'uins  lo  |8|  can  he 
developed  from  the  purely  theoretical  standpoint,  a  favorite 
solution  IS  to  consider  an  initial  value  problem  for  an  impul¬ 
sive  source  delincd  over  a  semi-inlinite  time  domain  and  valid 
fivr  the  sp.u.r  o(  iltc  enure  radiation  bell  Such  a  solution  is 
easily  wrillen  formally  in  terms  of  an  eigenfunction  expansion 
of  the  Green  s  function  for  the  equation  (e  g  .  C'liiu  and  Hiliun, 


Fig  3  Schematic  illustration  of  the  boundary  volume  in  |v.  i| 
space  in  which  the  general  bimodal  diffusion  solution  is  developed 
Saielliie  data  lo  be  matched  to  the  solution  span  lie  on  the  daia 
surface 


1977],  However,  we  can  think  of  several  factors  that  vitiate 
the  meaning  of  such  a  solution  and  its  relevance  to  data  appli¬ 
cation.  First,  radiation  belt  diffusion  processes  are  currently 
not  well  known  enough  to  permit  single  forms  and  values  of 
the  diffusion  coefficients  to  be  specified  for  the  enure  bell, 
after  all  the  purpose  of  data  application  is  to  determine  the 
ensemble  of  coefficients.  Second,  the  source  time  scale  is  niH 
well  defined  for  radiation  bell  particles  What  seems  to  be 
impulsive  in  one  phenomenon  need  not  be  so  when  beheld 
with  the  vantage  point  of  another.  An  important  factor  which 
determines  the  appropriate  solution  to  (8)  is  the  space-time 
span  of  the  data  set.  Here  two  features  must  be  recognized 
with  the  help  of  Figure  3.  This  figure  illustrates  the  satellite 
data  span  and  its  relation  lo  the  appropriate  solution  volume 
in  (x,  L.  I)  space.  The  fundamental  nature  of  nonsynchronous 
satellite  data  spans  is  that  they  lie  on  a  data  surface  {ABCD  on 
the  figure)  defined  by  the  satellite  trajectory  Ui)-  Further,  the 
lime  span  defines  the  (L.  x)  boundary  surfaces  of  the  solution 
volume  labelled  by  the  Roman  letters  A  lo  H  on  the  figure  In 
such  a  source-free  volume,  particles  are  transported  into  and 
out  of  the  surfaces  ABCH  and  CDEF  by  radial  dilfusion.  in 
other  words,  these  boundary  surfaces  are  also  source  surfaces 
Particles  are  lost  through  the  surface  BdC'F  which  sliouUI 
rigorously  lie  on  x,  as  discussed  above.  If  the  data  time  span  is 
less  than  or  comparable  to  the  injection  and  loss  times  in  the 
volume,  as  is  the  case  in  most  satellite  observations  such  as 
SCATHA.  the  appropriate  solution  to  be  developed  is  the 
general  source-free  complete-set  expansion  solution  in  a  finite 
space-lime  interval,  and  not  an  initial  value  solution  in  the 
semi-infinite  time  domain  corresponding  lo  response  lo  a  tran¬ 
sient  (i.c.,  della  function)  source.  In  the  finite  spacc-lmtc 
domain,  the  entire  .spectrum  of  lime  constants  (rising  as  well  as 
falling)  must  tic  included  in  the  general  complete-set  expansion 
because  unknown  sources  (and  sinks)  outside  of  the  volume 
are  diffusing  particles  into  (and  out  of)  the  finite  volume  caus¬ 
ing  the  phase  space  density  or  distribution  function  lo  increase 
(decrease)  with  time  scales  that  can  be  comparable  to  the  d.il.i 
time  span;  whereas,  in  the  semi-infinite  time  domain  for  the 
Green's  function,  only  the  falling  time  response  function  is 
appropriate  To  be  sure,  the  Green's  funcion  can  be  applied 
lo  obtain  the  finite  space-time  solution  by  applying  Gieen's 
theorem  onto  the  boundary  surfaces,  which  involve  inicgr.i 
tion  of  Green's  function  over  the  surface  bounilary  fluxes  or 
their  derivative  whenever  appropriate  I  his  crucial  inicgi.iiion 
imparls  the  source  lime  scales  onto  the  solution  I  be  impi't 
lanl  pxiinl  is  that  this  integration  convolves  the  rising  and 
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falling  lime  constants  of  the  surface  sources  with  the  natural 
falling  time  constants  of  dilfusion  to  result  in  a  spectrum  of 
both  rising  and  falling  lime  constants  for  the  spatial  eigenfunc¬ 
tions.  As  a  result  of  this  lime  convolution,  the  hnal  solution 
for  linite  volume  need  not  look  anything  like  the  Green's  func¬ 
tion  that  would  be  developed  for  the  initial  value  problem 
over  scmi-inlinite  time.  To  be  sure,  all  excited  stales  eventually 
decay  by  dilTusion  to  equilibrium,  but  the  satellite,  with  its 
finite  space-time  data  span,  does  not  know  whether  it  has 
caught  the  rising  or  falling  phase  of  a  modal  disturbance  of 
unknown  time  scale. 

These  esoteric  characteristics  of  solutions  to  partial  differ¬ 
ential  equations  are  familiar  enough  to  the  theorist;  however, 
it  may  be  relevant  here  to  relate  the  situation  to  a  more  intui¬ 
tive  and  mundane  example.  Consider  the  case  of  heat  conduc¬ 
tion  in  a  bar  of  metal  occupying  the  volume  ABCDEFCH 
above.  The  heat  diffusion  equation  is  similar  in  structure  to  (8) 
and  a  transient  response  to  a  heal  pulse  (Green's  function)  is 
composed  of  a  spectrum  of  sinusoids  all  exponentially  decreas¬ 
ing  in  lime  with  the  appropriate  rales.  However,  in  order  to 
obtain  the  solution  for  the  case  of  an  impressed  temperature 
held  on  the  surface  ABGH  that  is  rising  in  time,  say  exponen¬ 
tially.  the  spectrum  of  sinusoidal  modes  of  heat  conduction 
will  ri.se  or  fall  in  lime  each  according  to  their  ability  to  con¬ 
duct  the  heat  to  the  volume  and  to  the  boundary  surfaces. 
Intuitively,  it  is  easy  to  imagine  that  if  (he  impressed  temper¬ 
ature  is  sulTiciently  high  it  may  cause  the  temperature  field  of 
portions  of  the  volume  to  rise  exponentially  for  a  finite  time, 
although  rising  time  constants  are  definitely  not  in  any  mode 
of  the  Green's  function  response.  This  case  of  transient  heat 
conduction,  and  the  relation  between  the  diffusion  Green's 
function  and  boundary  conditions,  are  given  by  Morse  and 
teshhach  [I95.T  equation  12.1.5,  p.  1587].  The  reader  can 
easily  carry  out  the  Green's  function  convolution  specified 
there  to  obtain  (he  time  response  to  an  impressed  temperature 
at  the  boundary  surface  and  demonstrate  that  over  finite  time 
some  modes  may  be  rising  while  others  are  falling. 

The  consequence  of  the  esoteric  points  discussed  above  is 
that  It  would  be  wrong  to  force  satellite  data  obtained  over  a 
finite  space-time  volume  to  fit  an  impulsive  transient  response 
solution  such  as  the  Green's  function.  For  purposes  of  satellite 
data  analysis  in  a  finite  space-time  interval,  the  appropriate 
solution  to  (8)  is  the  solution  corresponding  to  a  general 
boundary  value  problem,  and  not  corresponding  to  a  global 
response  to  an  impulsive  source.  Since  a  general  solution  can 
be  obtained  by  eigenfunction  expansion  of  (8)  with  unrestric¬ 
ted  separation  constants,  which  will  be  constrained  by  data, 
this  convenient  properly  forms  the  basic  principle  of  data 
chatacicri/alion  analysts  in  this  paper.  It  will  be  shown  later, 
based  on  the  success  of  our  data  analysis,  that  this  principle 
can  be  used  to  characterize  the  stale  of  the  entire  radiation 
belt  in  a  piecewise-conlinuous  fashion. 

tinder  the  precautions  discussed  above,  a  general  solution 
to  18)  in  limie  space-time  can  be  developed  [cf  Morse  and 
it'shhoik.  195.11  m  terms  of  complete-set  eigenfunction  ex¬ 
pansions  In  general,  there  are  three  modal  expansions:  (I) 
single-mode  pitch  angle  diffusion.  (21  single-mode  radial  diffu¬ 
sion  and  l.f|  simultaneous  bimodal  diffusion  Using  the  device 
of  the  Kronecker  della  <he  three  modal  expansions  can  be 
written  in  terms  of  a  single  eigenfunction  expansion 
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where  k,  is  the  nth  zero  of  the  Bessel  function  The  parame¬ 
ters  and  new  functions  are  defined  as  follows: 

c  =  3(^l  -^j/(/i-2) 


y,(t)  =  f  ■  [J , , .(I  /?.  I  /L-')]  .  -  [/ ,  , .(I  /? J  / L'* )]  .  (II) 

Z,(L)  =  L*  [N , .  .(I  p,  I  /L^)]  .  •  [  A , .  .(I  i  //-•')]  (12) 


J  is  the  regular  Bessel  function,  /  and  K  are  the  modified 
Bessel  functions,  and  Af  is  the  Neumann  function.  In  the 
above,  the  separation  constants  are  r,  and  p,.  The  usual  nota¬ 
tion  [(/]>  is  defined  in  this  case  by 
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and  analogously  for  [(/].,.  The  separation  constants  and  pa¬ 
rameters  of  the  solution  are  related  by  the  algebraic  relations 
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which  will  be  referred  to  as  separation  equations.  It  is  seen 
that  the  single  pitch  angle  mode  time  constants  are  negative  as 
expected,  since  single-mode  pitch  angle  diffusion  always  leads 
to  decay.  For  bimodal  and  radial  diffusion,  however,  the  time 
constants  i,  can  be  either  positive  or  negative  over  a  linite 
time  interval.  This  is  due  to  the  fact  that  the  phase  space 
density  in  the  finite  L  domain  can  rise  or  fall  in  time  in  re¬ 
sponse  to  particle  transport  into  and  out  of  the  adjacent  vol¬ 
umes,  as  discussed  above.  With  the  sign  of  each  (,  not  prede¬ 
termined,  the  separation  constant  p,^  can  be  positive  or  nega¬ 
tive.  This  is  the  origin  of  the  (7.  I)  and  IN.  K)  options  of  (II) 
and  (12).  The  solution  (9)  is  not  the  most  general  solution,  for 
a  separate  sum  of  pure  radial  modes  can  be  added.  We  have 
not  done  so  for  two  reasons:  (I)  we  expect  that  episodes  of 
radial  diffusion  must  accompany  pitch  angle  diffusion  because 
of  the  association  of  magnetospheric  fluctuations  of  all  scales 
with  each  other  and  (2)  this  creates  more  parameters.  As  it 
presently  stands,  only  one  pure  radial  mode  is  included,  i  e.. 
for  n  «  0 

Finally,  we  return  to  the  second  feature  of  satellite  data 
illustrated  on  Figure  3.  In  a  data  analysis  effort,  the  (it  of 
satellite  data  is  on  the  IL.  l)  data  surface  ABCD.  which  sam¬ 
ples  the  boundary  surfaces  and  the  volume  in  a  wav  that 
mixes  space  and  time.  The  fit  to  the  general  solution  i9i  is 
supposed  to  yield  the  parameters  (a,.  A,.  i„.  f..  r,i.  which  would 
allow  lime  constants  and  diffusion  coefficients  il),,.  O,  ,  i  to  be 
determined.  Further,  the  solution  in  the  enure  volume  is  sup¬ 
posed  to  be  determined  via  the  expansion  coefTicienis  (u...  h,. 
c,)  Clearly,  the  solution  to  (9|  is  not  specihed  b>  the  data  m  a 
wav  that  is  usual  to  a  boundary  value  problem  Bcc.iusc  of  the 
mixing  between  L  and  r  on  (he  data  surface,  it  is  difficult  lo 
devise  a  completely  unambiguous  procedure  to  determine  all 
of  the  fit  parameters  An  approximate  procedure  that  works 
will  be  outlined  in  section  5  The  generality  of  (9|  must  not 
overshadow  the  restrictions  under  which  |9)  is  derived,  i  e  .  the 
assumption  of  separability  of  the  v  and  1.  operators  Thus,  the 
data  analysis  must  be  contemplated  in  tcriiis  of  a  (iiccewisc 
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procedure,  dealing  willi  an  (L.  (f  volume  in  which  the  quasi- 
invariance  of  .T  is  preserved  Our  attempt  here  is  to  approach 
ihe  specification  with  an  appropriate  piecewise  method,  es- 
i.ihhshing  parameter  sets  h,.  c„.  r„)  for  many  cases,  from 

which.  It  IS  hoped,  a  "learning  curve"  will  allow  us  to  deter¬ 
mine  approximate  specifications  of  the  radiation  belts  under 
various  classes  of  conditions. 

The  development  up  to  this  point  under  the  simplifying 
assumptions  is  strictly  theoretical  and  exact.  From  this  point 
onward,  we  shall  attempt  to  test  whether  bimodal  theory  (and 
Its  solution)  IS  a  good  description  of  the  data.  This  confron¬ 
tation  IS  the  main  objective  of  this  paper.  It  may  be  referred  to 
as  a  "fit.”  which  is  really  a  misnomer  since  our  primary  objec¬ 
tive  from  this  point  onward  is  to  see  if  fundamental  limitations 
of  the  data  in  mixed  space-time  can  be  overcome  in  the  face  of 
purely  theoretical  requirements  of  the  exact  solution.  No  over¬ 
all  sophisticated  "htting"  has  been  attempted  in  the  radial 
dilTusion  segment  of  the  analysis  Indeed,  we  deliberately 
avoided  "fitting"  in  this  segment  in  order  to  determine  the 
maximal  error  spread  of  the  exact  bi-modal  diffusion  descrip¬ 
tion  of  the  data  base.  This  maximal  data  spread  will  be  illus¬ 
trated  in  two  ways:  1 1 1  as  error  spreads  in  the  determination  of 
the  radial  diffusion  coefficient,  and  (2)  as  error  in  the  three- 
dimensional  reconstitution  of  the  entire  data  base. 

4  APFIICATION  to  Data  REPRFSENTAtlON 

Under  the  ass.  mption  that  bimodal  d  lfusion  is  operative, 
the  application  of  the  general  solution  (9)  to  satellite  measure¬ 
ments  involves  yet  other  approximations  and  artifices  that  do 
not  relate  to  the  basic  mathematical  and  physical  structures  of 
the  solutions  themselves  Rather,  these  arise  chiefly  from  the 
limitaiions  of  satellite  measurements 

Data  limitations  due  to  instrument  and  orbit  coverage  have 
been  addressed  in  section  2.  Limitations  on  space-time  inter¬ 
vals  of  coverage  have  been  addres.sed  in  section  3.  These  limi¬ 
tations  are  relatively  easier  to  deal  with  than  the  fundamental 
mixing  of  space  and  time  in  satellite  observations  due  to  the 
orbiting  motion  of  satellites.  For  radiation  bell  observations, 
it  IS  somewhat  unfortunate  that  the  timescale  of  the  order  of 
hours  IS  common  to  both  physical  phenomena  of  interest  in 
diffusion  theory  and  to  satellite  orbits  Thus  additional  ap¬ 
proximations  and  assumptions  have  to  be  made  to  bring 
theory  and  data  into  contact  with  each  other  The  use  of 
invariant  space  to  examine  the  corresponding  space-lime  be¬ 
havior  of  the  distribution  function  will  help  us  resolve  the 
mixing  of  space  and  time  in  the  data  This  technique  of  map¬ 
ping  between  (L.  x,  r|  and  invariant  space  is.  however,  the 
topic  of  another  paper  For  the  present  analysis,  we  shall 
assume  a  procedural  chain  which  will  be  described  below  We 
shall  depend  on  the  internal  consistency  of  our  results  and 
agreement  with  analyses  of  other  data  sets  to  justify  our  pro¬ 
cedure 

The  initial  steps  of  ihe  procedural  chain  for  data  analysis 
under  the  guidance  of  uillusion  theory  suggest  themselves  in 
the  eigenfunction  expansion  solution  (9)  After  apprevpriate  or- 
gau  /ation  and  prc'cessing,  the  set  of  radiation  belt  electron 
fluxes  such  as  that  shown  on  Figure  2  for  a  single  lO-min 
inicrval.  arc  converted  to  the  distribution  function  /  specified 
over  the  pitch  .ingle  and  energy  ranges,  and  over  the  intervals 
of  space-time  in  which  the  data  are  taken  The  pilch  angle 
diM.'ii’.,iions  arc  decomposed  into  amplitudes  y,  of  the  com¬ 
plete  set  of  UesscT  luiicUnns.  which,  under  our  staled  premises, 
are  also  the  complete  set  of  eigenfunctions  of  pitch  angle  diffu¬ 
sion  This  decomposition  of  Ihe  observed  distribution  function 
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yields  the  pilch  angle  amplitudes  y,,  which  are  nominally 
functions  of  L.  l  and  £.  The  complete-set  properly  of  the 
Bessel  functions  can  be  invoked  for  (16)  even  if  one  docs  not 
invoke  the  solution  (9). 

The  set  of  amplitudes  Q„.  separated  into  energy  channels  /;. 
must  now  be  resolved  into  separable  functions  of  lime  and  /. 
value,  as  indicated  in  the  second  major  bracket  in  (9)  This  is  a 
somewhat  ambiguous  task.  Several  simplifying  assumptions 
need  to  be  made,  even  just  to  ascertain  that  the  space  variable 
is  meaningfully  represented  by  £  in  a  sialic  magnetic  field 
model  over  Ihe  lime  span  in  question.  The  next  step  is  i.i 
assume  a  procedure  to  disentangle  the  mixing  of  L  and  /  by 
the  motion  of  the  satellite.  From  (9). 


y.(L.  n  =  c.(I  -  +  e"-[a,y,(£)  -l-  b,Z,(L)]  li", 

Because  of  the  mixing  due  to  satellite  motion,  spatial  elfcciv  ol 
scale  AL  such  that  (|t,|.  |(„|)  ^  AL/f,.  where  r,  is  satcTliic 
speed  perpendicular  to  the  L  shell,  cannot  be  distinguished 
from  temporal  effects.  However,  effects  for  which  Ihe  Iwn 
scales  are  distinctly  different  can  be  identified  in  the  (it.  To  til 
the  data,  an  ad  hoc  but  reasonable  choice  must  be  made  t.. 
identify  the  time  or  space  scales  to  select  or  to  sacrifice  in  the 
procedure.  For  the  SCATHA  data,  Ihe  coverage  is  typically 
AL  ~  I  —  2.  unlike  satellites  at  high-latitudc  elliptical  orbits 
which  cover  much  larger  AL.  Therefore  for  analysis  of 
SCATHA  data,  we  have  chosen  to  emphasize  the  hi  of  long 
lime  scales  and  short  L  scales,  i.e..  we  have  chosen  to  ascribe 
slow  variations  along  the  SCATHA  trajectory  to  temporal 
changes  and  Ihe  fast  variations  along  the  trajectory  to  spatial 
changes.  Our  procedure  to  effect  this  assumption  is  to  identily 
the  time  constants  to  each  order  by  first  doing  a  nonlinear  fit 
of  each  y,  to  a  simplified  version  of  (17)  in  which  the  data 
variations  over  the  entire  space-lime  span  arc  ascribed  to  time 
only.  This  enables  us  to  pick  out  the  time  scales  of  the  data, 
i.e.. 


Q,  =  <,ii  -  <5,0 »<■"■  + cy'"  ii><> 


where  r„.  r„,  <„  and  C„  arc  constant  parameters  for  this  step  If 
values  of  and  i„  arc  longer  than  or  of  order  of  Ihe  data  lime 
span,  which  happens  to  be  the  case  in  our  data  set,  we  can 
proceed  with  our  analysis  hypothesized  above.  If  not,  Ihe 
mixing  of  space  and  time  scales  must  be  disentangled  by 
invariant-space  mapping  techniques  and-or  the  use  of  proton 
data,  both  of  which  are  outside  the  .scope  of  this  paper.  Since 
this  step  in  the  fit  does  not  involve  the  L  dependence,  the 
amplitudes  C„  obtained  are  of  no  significance.  On  the  other 
hand,  since  e,  is  the  amplitude  for  single-mode  pilch  angle 
diffusion  it  is  not  expected  to  depend  on  L,  therefore  i,  is 
accepted  as  a  (it  parameter  The  time  constants  in  the  above  lit 
step  are  constrained  not  only  by  the  data  time  span  but  also 
by  (14)  and  (15).  The  constraint  (14)  on  is  somewhat  com¬ 
plex  and  will  be  discussed  below  fcqualion  (15).  on  the  other 
hand,  constrains  the  dependence  of  r,  on  n  according  to  the 
zeroes  of  the  Bessel  functions,  resulting  in  only  one  paramclei. 
D,,  This  constraint  is  easily  applied  on  the  nonlinear  lit  ( IS) 
Having  obtained  the  fit  parameters  (/7,,.  f,,  i,)  from  llic 
above  procedure,  we  next  obtain  a  lit  of  the  /.  depeiuleiRc 
according  to  Ihe  si'luluui  (<t)  lor  this  step,  we  consiriicl  .i 
residue  Hjl.)  by 


ii'h 


«.(£!  =  [y.lf-  'I  -  <.(l  "■ 


2f>2h 
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which,  according  lo  (17)  and  our  tit  hypothesis  above,  is  solely 
a  function  of  L. 

The  construction  of  this  residue  is  algebraically  straightfor¬ 
ward.  but  care  must  be  taken  not  to  amplify  the  natural  fluc¬ 
tuations  in  the  bracketed  dill'erence.  The  above  method  of 
derivation  forces  |r,|  and  |/,|  to  be  of  the  order  of  or  larger 
than  At.  i.e., 

(U.l.lfj)s  At  (20) 

This  constraint  also  assures  that  the  natural  fluctuations  in  the 
residue  RJL)  will  not  be  amplified  by  the  factor  exp  ( -  t/i,)  for 
I,  <  0.  The  criterion  (20)  allows  for  a  reasonable  unraveling  of 
the  space-time  scale  mixing  in  single  nonsynchronous  satellite 
observations. 

Hitherto  in  this  section,  we  have  treated  the  application  of 
the  general  solution  (9)  in  the  form  (16)  to  the  data  without 
making  any  reference  to  the  pro(>erty  that  the  expansion  is  a 
general  solution.  Rather,  we  have  been  exploiting  the  property 
that  it  is  a  complete-set  expansion  not  tied  to  a  specific  solu¬ 
tion.  In  other  words,  the  development  in  this  section  up  to  this 
point  can  be  effected  without  reference  to  the  bimodal  diffu¬ 
sion  equation,  because  (9)  can  be  viewed  as  a  general 
complete-set  expansion  when  all  the  (it  parameters  are  essen¬ 
tially  considered  free.  However,  our  aim  is  to  apply  the  solu¬ 
tion  (9),  which  is  of  similar  form  to  (16),  with  the  associated 
constraints  of  (9)  to  the  data. 

The  next  step  of  the  procedure  is  fitting  R,{L)  to  the  L 
eigenfunctions  of  the  diffusion  equation 

R.(L)  =  o,V.(L)  +  h,Z,(T)  (21) 

1  roni  ildl.  III).  (12)  and  (21).  we  note  that  the  parameters  of 
this  lit  arc  (r,  /(,.  <i,.  b,\  which,  together  with  t,.  arc  constrained 
by  (14).  I'tilizing  a  large  computer,  a  matrix  of  these  fit  pa¬ 
rameters  can  be  cycled  and  the  best  ht  chosen  for  each  energy 
channel  of  the  data.  Such  a  procedure  can  be  devised  as  a 
matrix  diagonalizing  routine.  Tor  our  test  case,  however,  the 
priority  is  not  to  attempt  such  large-scale  computations  at  the 
outset  but  to  examine  the  key  relationships  of  subsets  of  these 
parameters  lo  the  data  and.  foremost  among  the  various 
issues,  lo  demonstrate  that  a  physically  reasonable  fit  of  (he 
data  can  be  represented  by  the  general  bimodal  diffusion  solu¬ 
tion 

We  begin  the  procedure  by  reducing  the  number  of  parame¬ 
ters  lo  be  determined.  Previous  analyses  indicated  that 

2  <  <■  <  .7  [e  g..  WT'.s/  <7  ill,:  1981].  The  expected  value  of  v 
leads  lo  an  important  simplification  with  respect  lo  the  order 
of  the  Bessel  functions  in  (I  II  and  (12).  From  the  delinition  of 
I  III  the  first  of  the  triplet  ( 10).  we  obtain  —  3  16  <  /  <  0  which 
indicates  (hat  i  is  a  small  number  and  the  corresponding  order 
of  the  Bessel  functions  is  nearly  unity.  ITnm  (10).  the  order  is 


Since  the  Bessel  functions  are  analytic,  slowly  varying  func¬ 
tions  of  the  order,  we  choose  lo  constrain  the  fit  with  i;  =  0 
aiul  use  Bessel  functions  of  order  I  lo  approximate  those  of 
order  I  -E  i..  thus  resulting  in  a  large  reduction  of  compu¬ 
lation  With  these  simplifications  the  fit  parameters  are  now 
UK-  ",  f',)  which  arc  specified  by  energy  channel  T  and  order 
of  pilch  angle  cigcnmode  n  Once  the  til  for  the  enure  data  set 
IS  obi. lined,  the  parameters  are  used  lo  rcconslilule  the  data 
via  (9)  A  residue  distribution  is  then  cal,.ul.i'cd  lo  test  the 
goodness  of  til 


At  this  point,  the  analysis  can  be  viewed  on  two  related 
levels.  The  first  level  is  to  view  (9)  as  a  convenient  data  repre¬ 
sentation  tool  which  compresses  many  thous^.ids  of  data 
points  into  a  set  of  about  a  hundred  fit  parameters  amenable 
lo  .averaging  and  other  forms  of  processing  l  urlhcr.  the  tool 
isscpanilcd  but  not  divorced  from  the  solution  ul  the  bimodal 
diffusion  equation.  As  such,  this  data  representation  is  also 
convenient  for  interpretation,  and  for  construction  of  maps  of 
radiation  belts  under  various  conditions.  The  numerical  exam¬ 
ples  of  this  paper,  to  be  shown  below,  arc  at  this  level  of 
analysis.  In  this  approach,  the  final  product,  aside  from  de¬ 
monstrating  the  efficacy  of  bimodal  diffusion  in  specifying  the 
dynamics  of  the  outer  belt  electrons,  is  the  determination  of 
diffusion  coefficients  which  can  be  compared  to  theory.  The 
second  level,  a  natural  consequence  of  the  first,  is  to  view  (9) 
as  the  basis  of  modeling  the  dynamic  behavior  of  the  radiation 
belt  volume  appropriate  to  the  data.  The  data  fit  to  the  gener¬ 
al  solution  of  the  diffusion  equation  can  be  used  to  (ll  extract 
the  surface  sources  of  the  linite  space-time  volume  correspond¬ 
ing  to  the  solution  domain.  (2)  relate  the  ensemble  of  such 
extracted  boundary  conditions  to  the  set  of  diffusion  coef¬ 
ficients  from  a  set  of  satellite  data  spans,  and  (3)  integrate  the 
ensemble  of  such  verified  transport  parameters  for  the  ditlu- 
sion  equation  into  a  dynamic  model.  We  shall  limit  ourselves 
to  the  first,  less  ambitious,  level  of  application  in  this  paper 
The  data  set  consists  of  5040  radiation  belt  electron  Ilux 
values  derived  from  instrument  count  rates.  These  are  distrib¬ 
uted  among  the  seven  energy  channels,  the  24  ten-minute  flux 
averages  over  the  data  time  span  and  the  30  pitch  angle  bins 
of  3'  each.  i.e..  7  x  24  x  30  =  5040.  The  fit  procedure  is,  of 
course,  not  lo  deal  with  the  5040  Ilux  values  as  independent 
pieces  of  data,  for  these  are  organized  into  7  x  24  =  168  pilch 
angle  distributions  such  as  those  shown  on  Figure  2  These 
pitch  angle  distributions  are  subjected  to  the  Bessel  function 
decomposition  (16).  We  found  that  the  truncation  of  the  sum 
to  n  <  4  yielded  sufricicnl  lit  accuracy  without  causing  undue 
burden  on  a  VAX/780  computer  Next,  the  amplitudes  (1,  are 
fitted  nonlinear!)  according  lo  (18)  lo  extract  the  set  of  time 
constants  r,.  the  pilch  angle  dilfusion  coefficient  D,,  (via  the 
lime  constants  r,  and  the  constraint  (15)1  and  the  set  of  single¬ 
mode  pitch  angle  diffusion  amplitudes  c.  The  time  variations 
of  choice  are  thus  identified  in  the  data  In  the  next  step,  we 
assume  the  variations  of  the  residue  R,  in  (19)  to  be  in  /.  only; 
thus  proceeding  to  fit  the  set  of  residues  to  the  1.  eigenfunc¬ 
tions.  (211.  This  fit  of  the  residues  yields  the  eigenfunction 
parameters  (K  and  the  amplitude  parameters  tu,.  hj.  which 
completes  the  procedure  in  principle.  However,  these  parame¬ 
ters  are  not  entirely  free  because  the  final  constraint  l  I4i  must 
be  satisfied  I  he  natural  procedure  lo  satisfy  this  constraint, 
without  resorting  lo  the  matrix  diagonalization  method  dis¬ 
cussed  above,  is  lo  recycle  the  entire  procedure  by  using  |I4| 
as  the  adjustment  criterion  until  a  set  of  li,.  r,,  /(,  :l  that 
minimizes  the  variance  fiom  (14)  is  reached  This  can  be  done 
in  principle,  although  it  is  somewhat  tedious  l  or  this  test 
case,  however,  it  is  much  more  meaningful  to  show  the  con¬ 
vergence  of  our  procedure  by  displaying  the  v.iriance  Irom  i  Ni 
in  the  first  p.iss  of  the  procedure  I  o  demonstrate  the  van.ince, 
let  us  focus  on  (14)  as  an  equation  for  determining  ;.  .i  con¬ 
stant  independent  of  order  ii  In  anv  fit  tv'  a  constraint,  an 
equation  such  as  (14)  is  salishcd  lo  a  minimized  v.inaucc  all 
s>f  which  cun  be  ascribed  lo  .;  I  his  is  the  worse  case  fi'r  the 
variance  is  usually  shared  in  a  minimization  1  he  un- 
minimized  variance  of  .;  in  the  hrsi  pass  ('f  the  above  pri'- 
cedure.  determined  fiom  solving  ll4i  lor  all  the  orders  will 
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Fig  4  Comparison  between  observed  pilch-angle  distribution  for  the  830  keV  channel  at  L  =  5.42  (solid  curve)  and  its 
reconslituiion  after  decomposition  into  pitch  angle  eigenfunctions  (dash  curve).  The  fit  propagated  error  bars  are  shown. 


\ield  a  spread  of  the  values  of  4  as  function  of  order.  This 
spread  as  a  function  of  order  is  a  good  indicator  of  the  degree 
of  convergence  in  the  hrst  pass  of  the  fit.  Further,  since  {  is 
directly  proportional  to  the  radial  diffusion  coefficient,  the 
spread  of  the  values  determined  by  this  first  pass  of  the  pro¬ 
cedure  can  be  compared  vvitli  the  totality  of  values  and  vari¬ 
ances  of  the  radial  diffusion  cocfTicicnts  determined  in  the  past 
b\  methods  totally  unrelated  to  Ihts  analysis.  Thus  the  spread 
of  4  obtained  in  this  way  acts  as  an  independent  and  worse- 
case  indicator  of  the  goodness  of  fit  for  our  entire  analysis.  In 
this  regard,  it  is  crucial  to  note  that  4  does  not  appear  in  our 
solution  as  a  direct  fit  parameter,  rather  it  is  a  derived  quan¬ 
tity  due  to  the  constraint  (14). 

In  summary,  this  application  involves,  per  energy  channel.  6 
time  constants  (/„-4.  D,,).  14  amplitudes  (Uo-a-  ^o-a-  fi -a) 
and  5  eigenfunction  parameters  yielding  25  fit  parame¬ 

ters  per  energy  channel  or  a  total  of  175  fit  constants  for  the 
entire  data  set  These  are  to  be  compared  to  the  total  of  5040 
data  values  organized  into  168  distribution  functions.  It 
should  be  noted  that  no  attempt  at  organizing  the  fit  parame¬ 
ters  into  empirical  functions  of  energy  and/or  pilch  angle 
order  it  has  been  made  Further,  the  total  number  of  fit  pa¬ 
rameters  IS  exhausted  by  the  above  list;  thus  an  enlargement 
of  the  data  span  will  not  introduce  any  more  parameters  be¬ 
cause  the  above  exhausts  the  count  of  parameters  according  to 
theory  Our  ratio  of  fit  constants  to  data  values  gives  adequate 
constraints  on  the  lit  With  more  data  spans,  such  as  in  dy¬ 
namic  modeling,  the  energy  and  order  dependences  of  the  lit 
parameter  can  further  be  organized,  resulting  in  a  drastic  re¬ 
duction  in  the  number  of  parameters  to  data  values.  Flowever, 
as  It  IS.  the  results  of  the  fit  are  excellent. 

In  the  above  test  procedure,  special  attention  is  paid  to 
error  propagation  and  goodness  of  lit  at  each  intermediate 
stage  tach  individual  distribution  function  is  reconstituted 


from  the  fit  parameters  at  each  step  of  the  procedure  and 
compared  to  the  observed  distribution  function.  The  variances 
between  the  reconstituted  and  observed  distribution  functions 
are  calculated  and  minimized  before  going  to  the  next  step,  for 
which  the  cycle  is  again  repeated.  Figures  4  and  5  illustrate 
this  reconstitution  cycle  for  the  pitch-angle  distribution  of  the 
830-keV  channel.  In  Figure  4  the  comparison  between  ob¬ 
served  data  and  reconstitution  is  performed  for  the  pitch-angle 
decomposition  (16)  at  L  =  5.42  to  test  the  validity  of  trunca¬ 
tion  at  n  *  4.  Figure  5  shows  the  same  comparison  after  the 
entire  decomposition  and  fit  procedure  has  been  completed  for 
the  same  830-keV  channel  but  measured  at  /-  =  5.59  to  show 
the  nontrivial  variations  in  the  observed  data  and  the  fit. 
Compared  to  the  pitch  angle  distribution,  the  variance  be¬ 
tween  observed  data  and  fit  is  increased  in  the  complete  de¬ 
composition.  This  is  expected,  since  it  is  constrained  much 
more  by  the  solution  to  diffusion  theory:  however,  the  sub¬ 
stantial  changes  of  the  distribution  are  well  modeled  The 
errors  associated  with  the  fit  curves  in  these  two  figures  are 
errors  associated  with  the  data  which  have  been  propagated 
onto  the  fit  via  the  procedure.  The  original  errors  are  statis¬ 
tical  errors  of  the  instrument  count  rates  plus  3";.  data  com¬ 
pression  error  for  energy  channels  below  1.6  MeV.  In  addi¬ 
tion.  background  errors  are  included  for  channels  above  16 
MeV,  induced  by  background  counts  from  a  weak  radioactive 
calibration  source  in  the  instrument.  The  data  processing  soft¬ 
ware  does  make  a  correction  for  bremsstrahlung  contami¬ 
nation.  which  appears  primarily  at  448  keV,  caused  by  parti¬ 
cles  of  energy  above  1.6  MeV  partially  penetrating  the  detec¬ 
tor  shielding  Figure  6  shows  a  comparison  similar  to  Figure  5 
but  for  the  1616-keV  channel  at  L  =  5.42.  In  this  figure,  the 
original  errors  associated  with  the  observed  data  arc  shown 
for  comparison  with  the  fit  propagated  errors  in  the  previous 
figures.  Figure  7  shows  a  complete  decomposition  comparison 
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I'ifi  5  C'oinpiirison  helween  observed  pitch-angle  distribution  for  the  830-lieV  channel  al  L  =  S  59  (solid  curve)  and  ns 
rcconsiiiulion  afier  decomposiluin  inlo  Ihe  complete  solution  expansion  (dash  curve).  The  hi  propagaied  error  bars  are 
shown.  Note  the  differences  with  Figure  4. 


at  L  =  5.48  similar  to  Figure  5  again  but  performed  for  the 
lower  energy  channel  at  203  keV.  This  ht  is  worse  than  those 
shown  in  Figures  4-6.  We  speculate  that  convective  signatures 
may  be  present  in  these  lower  energy  channels  and  their  dy¬ 


namics  may  not  be  related  solely  to  diffusion,  which  is  w  hy  the 
energy  channels  below  256  keV  are  excluded  from  the  formal 
analysis.  Another  interpretation  is  that  storm  time  effects  were 
still  prevalent,  although  no  unusual  pitch  angle  distributions 


I  ig  A  .Same  as  Figure  5  estepi  ihc  comparison  is  for  Ihe  lAlh-keV  channel  al  t  ^  5  42  and  ihe  oiiginal  daia  erioi  bars 

are  shown. 
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[Siheik  el  ill .  1987]  were  found  in  our  data  base.  In  Figure  8. 
a  tally  of  the  fit  error  distributions  for  the  entire  set  of  data 
points  in  all  energy  channels  within  L  values  of  5.45  5.70  is 
shown  The  error  percentage.  100  x  |data  (it|  data,  for  each 
data  point  is  computed  for  each  data  value.  It  is  seen  that  the 
most  probable  error  is  about  5''.i  Wc  ludgc  this  to  be  ade¬ 
quate  accuracy  and  to  be  a  reasonable  basis  for  dynamic  mod¬ 
eling  of  the  radiation  belts.  With  the  method  of  analysis  devel¬ 
oped  above,  we  shall  be  able  to  distinguish  real  physical  vari¬ 
ations  of  the  radiation  belts  from  fit  noise  at  the  5'’''.  level. 

5.  Disc  ussioN  and  Summary 
In  the  foregoing,  we  have  demonstrated  that  a  solution  of 
the  bimodal  diffusion  equation  appropriate  to  a  finite  space¬ 
time  volume  describes  adequately  the  dynamics  of  outer-belt 
electrons  in  the  sp.icc-timc  span  of  a  selected  SCAFHA  data 
set  during  moder.ile  geomagnelic  activity  after  a  major  sub- 
storm  Although  simple-minded  assumptions  have  been  made 
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on  the  form  of  the  dilTusion  coefficients  in  order  to  reduce  the 
computational  cITori  and  although  the  computation  procedure 
IS  designed  to  test  and  illustrate  the  application  rather  than  to 
achieve  the  best  tit.  we  have  been  able  to  lit  the  entire  data  set 
to  a  general  solution  with  a  variance  level  of  about  5"n.  Our 
results  and  fit  accuracies  demonstrate  that  this  is  an  adequate 
and  potent  approach  to  dynamic  modeling  of  the  radiation 
bells  Obviously,  many  improvements  and  rehnements  can  be 
made  Some  of  these  will  be  dealt  w  ith  in  subsequent  work 
While  the  principal  goal  of  this  p.iper  is  to  demonstrate 
data  interpretation  using  the  general  solution  of  bimodal  dif¬ 
fusion  in  finite  space-lime  volumes  An  equally  important  goal 
IS  to  derive  the  physical  parameters,  i.c..  the  diffusion  coef- 
ticients.  and  to  assess  them  m  the  light  of  other  work  The 
pilch  angle  diffusion  cocfricicnl  />,,  is  derived  by  the  imposed 
fit  constraint  II.S|  independent  of  the  order  of  the  pitch  angle 
eigenfunction  Figure  9  shows  the  derived  together  with 
the  fit  errrirs  as  functions  of  energy.  In  the  solution.  I),,  is 
assumed  indepiendcnl  of  /.  in  order  to  achieve  separability 
f  rom  I  igurc  9.  it  is  seen  that  /2,,  is  remarkably  constant  with 
energy  also,  except  for  the  44X-kcA'  ch.innel,  which  has  large 
errors  possibly  introduced  by  brenisstrahlung  contamination, 
although  corrections  have  been  made  The  magnitude  of  /),, 
Ic.iding  to  time  sc.ilcs  via  1 1  5)  is  consonant  with  other  dctcrnii- 
nations  tsf  the  average  pitch  angle  diHusion  lilclimc  of  about  5 
days,  although  most  of  these  are  not  eigcnfunciion  specilic 
{c  c.  St  Iml:  iiiul  l.iinzcriil  1 1.  1974]  W  illtiiins  il  nl  [IhtiHjmc.i 
sored  the  rise  limes  to  half  m.ixim.i  of  electron  lliivcs  in  siinilai 
I.  shells  inimcdi.itclv  alter  several  magneiic  sli'ims  to  be  I  ' 
d.ivs  Magnetic  conditions  t>f  their  cases  arc  quite  dillciciil 
fri'in  ours  and  the  delinitions  of  time  scale  are  also  dillcieni 
from  ours  .Since  both  of  the  above  two  f.iclors  wr'iild  Iciul  to 
shorten  our  tinic  sc.ilcs  use  limes  ol  I  '  vl.ivs  m  It  i//i.iiiiv  it 
til  (  lUtvNj  c.iii  be  coin[s.ii,il>lc  lo  oiii  piltli  .ingle  ililliisioii  lime 
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Fig.  9.  Pilch  angle  diffusion  coefTicieni  as  function  of  energy. 

scales  The  L  dependence  in  their  storm  time  data  in  the  same 
L  shells  exhibited  a  variety  of  behavior,  showing  independence 
of  L  as  often  as  not;  therefore  it  is  dilTicult  to  sort  out  the 
energy  and  L  shell  dependence  of  0„  at  present. 

As  we  have  discussed  above,  the  derivation  of  the  radial 
dillusiun  coefficient  serves  the  dual  role  of  providing  for  physi¬ 
cal  data  in  the  outer  belt,  which  can  be  compared  with  other 
determinations  from  direct  particle  observations  [e  g..  Kava- 
iiau^h,  1968;  Newkirk  and  Wall,  I968<>;  Cornwall,  1968;  Lan- 
zeriiiii  el  a!.,  1970;  West  et  ai.  1981],  and  of  testing  indepen¬ 
dently  the  adequacy  of  our  application  to  the  data.  Figures  10 
and  1 1  show  scatterplots  of  0^.  for  the  energy  channels  indi¬ 
cated  and  for  L  values  of  5.3  and  6.1,  respectively.  As  has  been 
discus.sed  above,  scattered  values  of  result  from  our  delib¬ 
erately  not  enforcing  the  constraint  (14)  in  a  cycled  fit  pro¬ 
cedure  in  order  to  evaluate  the  worst  case  variance  of  the  fit 
by  solving  (14)  for  s  at  each  order.  This  deliberate  procedure 
produces  a  D,  ^  for  each  order  in  the  energy  channel,  account¬ 
ing  for  the  spread  of  values  on  these  figures  in  each  energy 
channel.  The  values  of  so  determined  have  a  spread  of  a 
factor  of  2-5.  with  the  largest  spreads  concentrated  in  the 
lower  energy  channels  where  we  have  noted  possible  convec- 


I  ig  It)  Radial  diffusion  ciieiricienl  as  funclion  of  energy  at 
s  to  The  spread  of  values  for  e.ish  energy  is  a  measure  of  the 
Horse-sase  variance  in  its  delermtnalion 


0  1  I  ,  ,  . _ L _ ■  ■  _ 

0  SOO  1000  1500  200i  25CO 
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Fig.  1 1.  Same  as  Figure  lOexcept  the  L  value  is  b  10 

lion  contamination.  Since  these  are  results  from  the  initial 
pass  of  the  procedure  without  enforcing  the  constraint  1 14).  we 
are  confident  that  a  reduction  of  the  spread  by  a  factor  of  at 
least  2  is  well  within  the  capabiiily  of  the  procedure  without 
jeopardizing  the  goodness  of  fit  demonstrated  in  the  last  sec¬ 
tion.  Still,  the  spreads  of  derived  in  this  initial  pass  arc 
already  considerably  smaller  than  most  previous  determi¬ 
nations  from  particle  data,  which  frequently  showed  spreads 
of  1-2  orders  of  magnitude  for  £)^^  within  the  same  data  set. 
Figure  12  shows  a  comparison  of  the  spread  of  values  of 
over  the  high-energy  range  determined  by  this  unconstrained 


Fig  12  Comparison  of  the  radial  diffusion  cocrTicieni  dciermmcd 
for  ihe  energy  range  of  0  8  2  0  MeV  (monied  area  sho» mg  decreasing 
shading  svilh  increasing  encigyl  »ilh  Ihe  tompilalion  of  1*  i 'i  er  .d 
(I9HI).  whose  individual  analysis  yields  curve  I  C  urves  2  8  arc  ev- 
perimenlal  cuve  2.  7 omu.v.vion  el  al  [1972).  curve  f  Fur/ci  [I9b')j 
curves  4  h.  Selkirk  anil  H  u/r  [1968hJ.  h allhamm,ir  [1970]  curve 
.\i-ivCirC  anil  Wall  [l968o].  curve  8.  Kaianauqh  [1968]  Curves  9  14 
are  Iheoreiical  or  semilheorelical  curve  9.  SakaJa  and  Mead  [  |96SI 
curve  10.  Treriknx  [1965].  curve  11.  Hirminaham  [1969],  turves  12, 
13,  Ciirnnall  [  1968] ,  curve  14  Hoi: north  and  Mo:er  (  1979  ) 
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Fig  13  Three-diitiensionjl  prcAenlalion  of  ihc  phase  space  densily  of  the  dala  base  and  ils  reconstiiulions  for  Ihe 
830-keV  channel  in  pilch  angle  and  space-lime  Figure  I3<i  shows  Ihe  dala  base.  Figure  I3f)  is  Ihe  reconslilulion  after 
Hlling  Ihe  pilch  angle  eigenfunclion  decomposiiion  (16)  Figure  I3f  is  Ihe  reconslilulion  afler  ihe  space-lime  unmixing 
procedure  and  radial  diffusion  analysis  (17)  have  been  applied  wiihout  ihe  benelils  of  a  filling  procedure. 


lesi  pass  plotted  as  function  of  L  and  compared  with  the 
up-to-date  compilation  of  all  determinations  given  by 
H'fs)  et  al.  [1981].  Note  that  the  determined  radial  dilTusion 
coefFicients  agree  well  with  the  storm  time  data  of  Timui.sMiiii 
el  al  [1972].  Since  Ihe  data  on  Figure  12  are  from  a  mixiure 
of  energies,  it  is  difFicult  to  test  the  energy  dependence  of  Ihc 
radial  dilTusion  coefFicienI,  which  decreases  with  increasing 
shading  (lower  energy)  Thus  this  independent  test  of  the  ade¬ 
quacy  of  the  general  solution  to  bimodal  diffusion  to  represent 
outer  bell  electron  data  is  not  only  satisfactory  but  actually 
indicates  an  improvement  over  previous  analyses,  which  were 
restricted  to  single-mode  diffusion 

Finally,  to  demonstrate  the  multidimensional  nature  of  our 
analysis  and  the  ability  of  Ihe  exact  bimodal  solution  to  repre- 
senl  dala.  we  show  on  Figure  1.3  Ihe  specification  of  the  entire 
dala  base  of  the  83()-keV  channel  m  pilch  angle  and  space¬ 
time  (Figure  I3ul.  ils  reconslilulion  after  filling  of  Ihc  pilch 
angle  eigenfunction  expansion  (Figure  13b)  and  its  reconslilii- 
lion  from  ihc  complete  eigenfunction  expansion  after  Ihc  ad 
hoc  unmixing  of  space  and  lime  in  Ihe  radial  diffusion  seg¬ 
ment  of  Ihe  .inalvsis  (F  igure  1 3< )  As  is  discussed  above,  a 
filling  procedure  has  not  been  applied  to  Figure  13/c  in  order 


to  determine  Ihe  maximal  spread  of  the  radial  diffusion  coef¬ 
ficients;  therefore.  Ihe  error  in  Figure  13c,  relative  to  Figure 
1 3o.  reflects  Ihc  error  distribution  in  Figure  8  and  manifests 
Itself  as  the  spread  of  radial  diffusion  coefficients  shown  in 
Figures  10  and  )  I  This  unfitted  reconstruction  explicitly  dem¬ 
onstrates  the  effects  of  the  liigh-pcrccntage  tail  of  Figure  8 
The  basic  trend  of  the  data  reconstruction  reflects  the  effects 
of  the  most  probable  part  (5°;,)  of  Ihc  error  distribution  in 
Figure  8.  It  is  clear  from  a  comparison  of  Figures  13a  and  l.3< 
that  their  basic  trends  agree.  It  is  also  likely  that  the  increased 
error  in  the  unfitted  second  (radial)  phase  of  the  analysis  orig¬ 
inates  from  Ihc  ad  hoc  unmixing  of  space  and  lime  variables 
necessary  in  any  multidimensional  analysis  of  dala  from  a 
moving  satellite.  Further  investigations  on  the  unmixing  pro¬ 
cedure  may  be  necessary,  if  a  filling  procedure  does  not  reduce 
the  unfitted  errors 

The  ultimate  goal  of  an  analysis  such  as  in  this  paper  is  to 
obtain  a  proper  theoretical  foundation  for  dynamic  modeling 
of  the  radiation  bells  We  have  not  achieved  this  goal  by  far, 
but  we  have  demonstrated  that  our  approach  has  vast  polcn- 
li.il  to  qiianlilativcly  inicrprcl  dala  and  achieve  a  higher  level 
of  ■Kcii'';icy  ■I  determination  of  physical  parameters  lo 
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achieve  dynamic  modeling  in  the  outer  belt,  an  important 
hurdle,  alluded  to  previously,  is  that  the  L  shells  are  time- 
dependent  even  during  moderate  geomagnetic  activity,  re¬ 
sulting  in  not  just  another  scale  to  be  di.sentangled  but  also  in 
dra.s:iw  cITccts  of  diift  loss  through  the  magnetopause  bound¬ 
ary.  Loss  boiiiiil.iries  can  he  iiccommodalcd  niilurully  in  the 
eigenfunction  solution  formulism;  however,  the  choice  of 
phase  space  as  representation  coordinates  may  not  be  as 
powerful  as  the  use  of  invariant  space  in  this  regard.  Our 
future  efforts  will  be  directed  toward  this  area. 
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HIGH-ENERGY  OUTER  RADIATION  BELT  DYNAMIC  MODELING 

Y.  T.  Chiu,  R.  VV.  Nightingale  and  M.  A.  RinaJdi 
Lockheed  Palo  Alto  Research  Laboratory 
3251  Hanover  St.,  Palo  Alto,  Ca.  94304 

ABSTRACT 

Specification  of  the  average  high-energy  radiation  belt  environment  in  terms 
of  phenomenological  montages  of  satellite  measurements  has  been  available  for 
some  time.  However,  for  many  reasons  both  scientific  and  applicational  (includ¬ 
ing  concerns  for  a  better  understanding  of  the  high-energy  radiation  background 
in  space),  it  is  desirable  to  model  the  dynamic  response  of  the  high-energy  ra¬ 
diation  belts  to  sources,  to  losses,  and  to  geomagnetic  activity.  Indeed,  in  the 
outer  electron  belt,  this  is  the  only  mode  of  modeling  that  can  handle  the  large 
intensity  fluctuations.  Anticipating  the  dynamic  modeling  objective  of  the  up¬ 
coming  Combined  Release  and  Radiation  Effects  Satellite  (CRRES)  program, 
we  have  undertaken  to  initiate  a  study  of  the  various  essential  elements  in  con¬ 
structing  a  dynamic  radiation  belt  model  based  on  interpretation  of  satellite 
data  according  to  simultaneous  radial  and  pitch-angle  diffusion  theory.  In  order 
to  prepare  for  the  dynamic  radiation  belt  modeling  based  on  a  large  data  set 
spanning  a  relatively  large  segment  of  L- values,  such  as  required  for  CRRES,  it 
is  important  to  study  a  number  of  test  cases  with  data  of  similar  characteristics 
but  more  restricted  in  space-time  coverage.  In  this  way,  models  of  increasing 
comprehensiveness  can  be  built  up  from  the  experience  of  elucidating  the  dy¬ 
namics  of  more  restrictive  data  sets.  The  principal  objectives  of  this  paper  are 
to  discuss  issues  concerning  dynamic  modeling  in  general  and  to  summarise  in 
particular  the  good  results  of  an  initial  attempt  at  constructing  the  dynamics  of 
the  outer  electron  radiation  belt  based  on  a  moderately  active  data  period  from 
Lockheed’s  SC-3  instrument  flown  on  board  the  SCATHA  (P78-2)  spacecraft. 
Further,  we  shall  discuss  the  issues  brought  out  and  lessons  learned  in  this  test 
case. 

FORMULATION 

Some  two  decades  of  satellite  observations  have  shown  that  the  electron  ra¬ 
diation  flux  (including  the  satellite-damaging  electrons  of  >500  keV  energy) 
in  the  outer  belt  fluctuates  by  orders  of  magnitude  in  response  to  geomag¬ 
netic  activity.  Some  typical  fluxes  and  spectra  measured  by  the  SCATHA  SC-3 
instrument'**'*  are  shown  in  Figure  1.  The  practical  need  to  characterise  and 
predict  these  large  radiation  fluctuations  for  satellite  survivability  considerations 
requires  an  understanding  of  their  dynamics  beyond  phenomenological  corre¬ 
lations;  thus,  the  need  for  dynamic  modeling  of  the  outer  electron  belts  was 
initiated. 

In  its  most  comprehensive  form  dynamic  modeling  involves  a  theoretical 
accounting  of  the  flux  responses  in  terms  of  fundamental  processes  of  electron 
transport  and  loss  in  the  magnetosphere.  Since  wc  are  dealing  with  energetic 
electrons,  a  Fokker-Planck  description  of  the  transport  processes  naturally  leads 
to  a  diffusion  picture  in  terms  of  several  independent  variables  characterizing 
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Diffusive  transport  of  energetic  particles  in  the  magnetosphere  is  an  excel¬ 
lent  theoretical  approach  to  outer-belt  dynamic  modeling,  but  the  latter  involves 
much  more.  The  constraint  of  satellite  data  and  conditions  of  the  magneto¬ 
sphere  play  roles  in  dynamic  modeling  which  are  as  important  as  the  diffusion 
theory  itself.  In  the  outer  belt  the  magnetic  field  is  transient.  Plasma  and 
field  conditions  encountered  by  a  drifting  electron  vary  drastically  and  parts  of 
the  energetic  particle  population  may  even  cross  outside  of  the  magnetosphere. 
Thus,  dynamic  modeling  involves  an  accounting  of  fields  and  magnetospheric 
boundary  configurations  as  much  as  diffusive  particle  transport.  All  of  these 
elements  are  integral  parts  of  a  dynamic  model  of  the  outer  electron  belt. 


Figure  1.  Sample  spectra  showing  the  variation  of  outer-belt  elec¬ 
trons  obtained  in  1979  and  1980  by  the  Lockheed  SC3 
instrument.  Phenomenological  models  based  on  6-month 
averages  of  the  data  can  be  in  error  by  an  order  of  mag¬ 
nitude  or  more. 
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To  illustrate  the  complexity  that  a  dynamic  model  must  deal  with,  Fig¬ 
ure  2  shows  a  survey  of  the  high-energy  electron  characteristics  encountered  by 
SCATHA  in  an  active  period  on  day  115  of  1979  in  which  a  moderate  magnetic 
storm  occurred,  with  a  sudden  commencement  just  prior  to  the  start  of  the  day 
and  a  maximum  Dst  of  -148  7.  Just  before  18:00  h  LT  there  was  a  brief  en¬ 
hancement  of  electron  fluxes  that  may  have  been  an  adiabatic  compression  event. 
“Butterfly”  pitch-angle  distributions,  indicating  loss  of  90°  pitch-angle  particles 
at  the  magnetospause,  were  seen  at  sJl  energies  for  almost  2  hours,  after  which 
the  lower-energy  electron  fluxes  recovered  to  their  previous  levels.  Shortly  after 


Magnetopause  Scatma  Orbit 


Figure  2.  A  montage  of  electron  flux  characteristics  encountered 
by  SCATHA  on  day  115  in  1979.  The  inset  figures  are 
schematics  of  pitch-angle  distribution  with  90°  pitch  an¬ 
gle  placed  at  the  center. 
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midnight,  there  began  a  succession  of  substorms,  with  deep  flux  dropouts  over 
the  portions  of  the  orbit  shown  as  dotted  lines  in  Figure  2.  The  dropouts  exhibit 
the  expected  correlation  with  the  geomagnetic  activity  indices,  particularly  AE. 
These  were  followed  by  rises  in  the  flux  levels  that  can  be  interpreted  as  a  fresh 
arrival  of  electrons  on  the  observed  flux  tubes  by  either  injection  or  compres- 
sional  changes  of  the  magnetic  field.  Such  active  periods  are  usually  recognized 
by  decreases  in  the  total  flux  and  the  occasional  appearance  of  nearly  isotropic 
pitch-angle  distributions  as  shown  on  the  Figure.  In  the  morning  hours,  typical 
“hat”  distributions,  indicating  enhanced  fluxes  of  90°  pitch-angle  particles,  are 
encountered. 

In  view  of  the  complexity  of  the  above  data,  it  is  sensible  to  separate  the 
elements  of  outer-belt  dynamic  modeling  into  modules  which  can  be  individually 
formulated.  Thus,  our  approach  to  dynamic  modeling  is  to  initially  treat  each 
of  the  following  elements  as  if  they  are  not  related,  but  bearing  in  mind  all  the 
time  that  they  are  to  be  integrated  into  a  whole.  These  elements  are: 

(i)  Transport,  injection  and  loss  of  trapped  energetic  electrons  in  terms  of  dif¬ 
fusion  in  phase  space; 

(ii)  Description  of  the  space-time  variations  of  the  geomagnetic  field; 

(iii)  Physics  of  transport  parameters. 

The  principal  topics  of  this  paper  deal  only  with  issues  arising  out  of  the 
diffusion  physics  (item  (i)].  The  diffusion  physics,  which  is  the  fundamental  ba¬ 
sis  of  a  dynamic  model,  appears  initially  to  be  mostly  straightforward,  although 
mathematically  cumbersome.  However,  it  turns  out  that  the  boundary  condi¬ 
tions  corresponding  to  satellite  data  can  become  a  major  issue  to  be  resolved 
irrespective  of  the  chosen  diffusion  formulation.  Since  the  boundary  data  issue 
is  central  to  sdl  schemes  of  dynamic  modeling,  we  discuss  it  at  the  outset. 

On  the  least  ambitious  level,  dynamic  modeling  can  be  viewed  as  a  means 
of  projecting  radiation  belt  data  obtained  along  a  satellite  track  by  physical 
means  to  specify  the  radiation  environment  throughout  some  reasonable  re¬ 
gions  not  on  the  satellite  track.  As  such,  the  selected  diffusion  treatment  is  the 
projection  mechanism  and  the  in-track  satellite  data  is  the  raw  material  to  be 
processed:  the  “boundary”  values  for  the  diffusion  physics.  However,  satellite 
measurements,  obtauned  on  a  space-time  trajectory,  are  not  standard  boundary 
conditions  stated  at  fixed  time  or  space  coordinates  to  define  either  initial-value 
problems  or  boundary- value  problems  of  diffusion  theory.  The  situation  for  dy¬ 
namic  modeling  is  different  from  these  simple  boundary-value  or  initial-value 
cases.  Figure  3  illustrates  the  satellite  data  spap  and  its  relation  to  the  appro¬ 
priate  solution  volume  in  (x,  L,t)  space,  where  i  is  the  cosine  of  the  pitch  angle 
and  L  is  the  McDwain  parameter.  The  fundamental  nature  of  nonsynchronous 
satellite  data  spans  is  that  they  lie  on  a  data  surface  (ABCD  on  the  figure) 
defined  by  the  satellite  trajectory  L{t).  Further,  the  time  span  defines  the  {£,  x) 
boundary  surfaces  of  the  solution  volume  labelled  by  the  Roman  letters  A  to  H 
on  the  figure.  In  such  a  source-free  volume,  particles  are  transported  into  and 
out  of  the  surfaces  ABGH  and  CDEF  by  radial  diffusion;  in  other  words,  these 
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Figure  3.  Schematic  illustration  of  the  boundary  volume  in  {x,L,t) 
space  in  which  the  general  bi-modaJ  diffusion  solution  is 
developed.  Satellite  data  to  be  matched  to  the  solution 
span  lie  on  the  data  surface. 

boundary  surfaces  are  also  source  surfaces.  Particles  are  lost  through  the  surface 
BGCF  through  pitch-angle  diffusion  into  the  loss  cone,  approximately  signified 
by  X  =  1  here.  If  the  data  time  span  is  less  than  or  comparable  to  the  injection 
and  loss  times  in  the  volume,  as  is  the  case  in  most  satellite  observations  such 
as  SCATHA,  the  appropriate  solution  to  be  developed  is  the  general  source-free 
solution  in  a  finite  space-time  interval,  and  not  an  initial-value  solution  in  the 
semi-infinite  time  domain  corresponding  to  response  to  a  transient  (i.e.,  delta 
function)  source.  In  the  finite  space-time  domain,  the  entire  spectrum  of  time 
constants  (rising  as  well  as  falling)  must  be  included  in  the  general  complete-set 
solution  expansion  because  unknown  sources  (and  sinks)  outside  of  the  volume 
are  diffusing  particles  into  (and  out  of)  the  finite  volume  causing  the  phase-space 
density  to  incre2ise  (decrease)  with  time  scales  that  can  be  comparable  to  the 
data  time  span.  In  contrast,  in  the  semi-infinite  time  ciomain  only  the  falling 
time  response  function  is  appropriate.  To  be  sure,  this  Green’s  function  can  be 
applied  to  obtain  the  finite  space-time  solution  by  applying  Green’s  theorem  onto 
the  boundary  surfaces*.  This  crucial  integration  imparts  the  source  time  scales 
onto  the  solution.  The  important  point  is  that  this  integration  convolves  the 
rising  and  falling  time  constants  of  the  surface  sources  with  the  natural  falling 
time  constants  of  diffusion  to  result  in  a  spectrum  of  both  rising  and  falling  time 
constants  for  the  spatial  eigenfunctions.  As  a  result  of  this  time  convolution,  the 
final  solution  for  finite  volume  need  not  look  anything  like  the  Green’s  function 
that  would  be  developed  for  the  initial-value  problem  over  semi-infinite  time.  To 


B-5 


be  sure,  all  excited  states  eventually  decay  by  diffusion  to  equilibrium,  but  with 
the  satellite  data  span,  over  a  finite  space-time  interval,  we  do  not  know  whether 
it  has  caught  the  rising  or  falling  phase  of  a  modal  disturbance  of  unknown  time 
scale.  For  the  SCATHA  orbit,  space  and  time  are  relatively  distinct  since  the 
range  of  L  covered  by  the  satellite  is  small.  However,  for  highly  elliptical  orbits 
(such  as  for  CRRES)  space  and  time  on  an  orbit  are  not  linearly  related.  We 
shall  see  that  under  such  circumstances  the  satellite  data  set  acts  as  a  highly 
nonlinear  boundary  for  the  dynamic  model. 

The  consequence  of  the  esoteric  points  discussed  above  is  that  it  would 
be  wrong  to  force  satellite  data  obtained  over  a  finite  space-time  volume  to 
fit  an  impulsive  transient  response  solution  such  as  the  Green’s  function.  For 
purposes  of  satellite  data  analysis  in  a  finite  space-time  interval,  the  appropriate 
solution  to  diffusion  theory  is  the  solution  corresponding  to  a  general  boundary- 
value  problem^,  and  not  corresponding  to  a  global  response  to  an  impulsive 
source.  Since  a  general  solution  can  be  obtained  by  eigenfunction  expansion 
with  unrestricted  separation  constants,  which  will  be  constrained  by  data,  this 
convenient  property  forms  the  basic  principle  of  data  charactrization  in  dynamic 
modeling.  It  will  be  shown  later,  based  on  the  success  of  our  data  analysis,  that 
this  principle  can  be  used  to  characterize  the  state  of  the  entire  radiation  belt 
in  a  piece-wise-continuous  fashion. 

Given  the  above  general  property  of  radiation  belt  dynamic  modeling,  let 
us  proceed  to  consider  our  test  case^  for  which  we  choose  to  specify  the  diffusion 
physics  in  a  static  magnetosphere  in  (z,L,t)  space  as  illustrated  in  Figure  3. 
[The  choice  of  spaces  to  characterize  radiation  belt  diffusion  is  not  arbitrary. 
This  complex  issue  is  treated  elsewhere^.)  Our  choice  of  diffusion  description 
is  that  of  simultaneous  multi-mode  diffusion  theory.  The  multimode  diffusion 
equation  can  be  simplified^’*'*  under  a  number  of  reasonable  assumptions.  One 
of  these  is  that  the  Z,-dependence  of  the  radial  diffusion  coefficient  Du  can 
be  parameterized  as  Du  =  ^  where  ^  and  fi  are  parameters  to  be  deter¬ 
mined  by  the  on-trajectory  satellite  data.  The  diffusion  equation  for  the  electron 
distribution  function  /  can  be  written 


dt 


■fz  + 

dz^  X  dz 


dL 


(1) 


where  i?,,  (the  pitch-angle  diffusion  coefficient),  (,  and  i/  are  considered  to 
be  constant  parameters. 

This  simplified  equation  is  to  be  solved  in  a  small  finite  domain  Lt  <  L  < 
I>u  expected  to  be  in  the  outer  belt  and  in  the  pitch-angle  domain  0  <  z  <  Ze, 
where  Zc  is  cosine  of  the  loss-cone  angle.  The  effect  of  pitch-angle  diffusion  is 
to  define  pitch-angle  behavior  of  /  inside  and  outside  of  the  loss-cone  boimdary 
z  =  Zc{L).  The  dependence  of  x^  on  L  is  yet  again  another  factor  which  mixes 
the  z  and  L  operators  in  (1);  however,  for  the  range  in  our  test  case  in  outer 
zone  region,  where  (L„  —  X/)  <  2,  z*  is  very  nearly  constant  and  nearly  unity. 
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We  assume  that  Xe  is  constant  in  the  data  span  so  that  separable  solutions  can 
be  obtained  for  (1)  with  the  present  test  case.  This  is  another  important  point 
when  we  deal  with  data  from  a  large  domain  of  L  (such  as  that  from  CRRES) 
which  involves  large  variations  of  the  loss  cone  both  from  the  large  variations  in 
L  and  in  distances  to  the  magnetic  mirror  points.  For  the  SCATHA  test  case 
the  loss-cone  dependence  on  L  is  entirely  negligible. 

Under  the  approximations  discussed  above,  a  general  solution  to  (1)  in  finite 
space-time  can  be  developed^  in  terms  of  complete-set  eigenfunction  expansions. 
In  general,  there  are  three  modal  expansions;  (1)  single-mode  pitch-angle  JT- 
fusion,  (2)  single-mode  radial  diffusion  and  (3)  simultaneous  bi-modal  diffusion. 
Using  the  device  of  the  Kronecker  delta  inoi  the  three  modal  expansions  can  be 
written  in  terms  of  a  single  eigenfunction  expansion 


S„,  +  (1  - 


JxiK) 


(2a) 


h{L,  t)  =  e‘/‘"  (a„y;.(L)  b„Z^{L)]  (2b) 

where  k„  is  the  zero  of  the  Bessel  function  Jo-  The  parameters  and  new 
functions  are  defined  as  follows: 


e  =  3(1  -  ^)/(m  -  2) 


y„(L)  =  I<'.[Ji+.(|^„|/L")]>-[/,+.(l/?nl/i")]<  (4) 

Z„(L)  =  ■  [yV,+.((  I  /L^)]^  .  /?„  (  fL^)]^  (5) 

J  is  the  regular  Bessel  function,  /  and  K  are  the  modified  Bessel  functions,  and 
N  is  the  Neumann  function.  In  the  above,  the  separation  constants  are  and 
0n.  The  usual  notation  [p]>  is  defined  in  this  case  by 


(ll]>  =9  if  >  0 
[5]>  =1  if  t„  <  0 


(6) 


and  analogously  for  [9]<.  The  separation  constants  and  parameters  of  the  solu¬ 
tion  are  related  by  the  algebraic  relations 
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tn'  =  (1  -  ^no)r-*  -  iPlX^ 


(7) 


T'^  =  -D,,kHxl  (8) 

which  will  be  referred  to  as  separation  equations.  It  is  seen  that  the  single 
pitch-angle  mode  time  constants  are  negative  as  expected,  since  single-mode 
pitch-angle  diffusion  always  leads  to  decay.  For  bi-modal  and  radial  Effusion, 
however,  the  time  constants  t„  can  be  either  positive  or  negative  over  a  finite 
time  interval.  This  is  due  to  the  fact  that  the  phase-space  density  in  the  finite 
L  domain  can  rise  or  fall  in  time  in  response  to  particle  transport  into  and 
out  of  the  adjacent  volumes,  as  discussed  above.  With  the  sign  of  each  tn  not 
predetermined,  the  separation  constant  0^  can  be  positive  or  negative.  This  is 
the  origin  of  the  {J,I)  and  (N,K)  options  of  (4)  and  (5).  The  solution  (2)  is 
not  the  most  general  solution,  for  a  separate  sum  of  pure  radial  modes  can  be 
added.  We  have  not  done  so  for  two  reasons:  (1)  we  expect  that  episodes  of 
radial  diffusion  must  accompany  pitch-angle  diffusion  because  of  the  association 
of  magnetospheric  fluctuations  of  all  scales  with  each  other  and  (2)  this  separate 
sum  creates  more  parameters.  As  it  presently  stands,  only  one  pure  radial  mode 
is  included,  i.  e.  for  n  =  0. 


Figuie  «. 

Illuctralion  of  the  nonlinear  nature  of  fitting  the 
representation  (2)  to  data.  In  all  panels,  the  func¬ 
tion  (2h),  it  plotted  with  the  same  set  of 

arbitrariljr  chosen  constant  parameters.  The  solid 
lines  are  contours  of  h[L,l)  and  the  dashed  lines 
are  W*)-  P»n«'  (»)  '*  independent  (t.t) 

grid.  Panel  (b)  it  for  L{t)  specified  by  SCATHA 
orbit  and  panel  (c)  is  for  the  CURES  orbit.  Note 
the  drastic  change  of  the  time  scale  required  in  or¬ 
der  for  the  CRRES  case  to  have  similar  contours  for 
the  h{L,t)  function  at  the  nearly  identical  panels 
(a)  and  (b). 
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Having  now  specified  a  general  solution  of  the  chosen  diffusion  theory,  we 
can  quantitively  illustrate  how  the  satellite  trajectory  L{t)  impacts  the  charac¬ 
terization  of  radiation  fluxes  away  from  the  trajectory.  The  idea  is  of  course 
central  to  dynamic  modeling  and  is  very  simple:  determine  the  parameters  of 
the  solution  with  on-trajectory  satellite  data  then  the  solution  yields  fluxes  at 
aJl  other  (x,L,<)  points.  But  as  pointed  out  above,  the  problem  is  that  the 
satellite  trajectory  mixes  L  and  1,  i.e.,  L(t),  which  complicates  the  fit.  For 
near-synchronous  satellites  such  as  SCATHA  for  which  only  a  small  range  of  L 
is  covered  in  an  orbit,  the  effect  of  space-time  mixing  is  quite  minimal.  This 
allows  us  to  complete  the  test  case  shown  in  the  next  section.  For  the  more 
comprehensive  situation  of  highly-elliptic  orbits  such  as  that  for  CRRES,  which 
attempts  to  specify  the  radiation  belts  over  a  wide  range  of  T,  the  trajectory 
L{t)  is  a  non-linear  function  of  t,  requiring  a  non-linear  fit  for  the  determination 
of  the  diffusion  parameters.  This  significent  complication  is  illustrated  in  Figure 
4  which  shows  plots  of  the  function  h{L,i)  in  (2b)  for  three  situations;  (a)  L 
and  t  independent,  (b)  SCATHA  orbit  on  which  L  is  nearly  independent  of  t 
and  (c)  CRRES-like  orbit  on  which  L  is  strongly  dependent  on  i.  The  strong 
distortions  on  panel  fc)  for  h(L{i),t)  relative  to  the  same  function  h{L,t)  on 
panel  (a)  signifies  that  dynamic  modeling  in  conjunction  with  a  satellite  data 
base  involves  much  more  than  obtaining  solutions  to  diffusion  theory. 
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Figure  5.  Ten-minute  averaged  electron  fluxes  integrated  over  all 
pitch  angles  for  7  energy  channels  measured  by  the  Lock¬ 
heed  SC3  instrument  on  board  SCATHA  during  day  165, 
1980.  This  is  the  data  set  to  be  analyzed  in  this  paper. 
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TEST  CASE  VERIFICATION 

Notwithstanding  the  difficulties  of  applying  radiation  belt  diffusion  theory 
to  dynamic  modeling  in  the  comprehensive  situation  of  the  expected  CRRES 
data  base,  it  is  of  interest  to  test  whether  multi-mode  diffusion  theory  ade¬ 
quately  describes  the  simpler  on-trajectory  SCATHA  data  base.  This  section 
summarizes  the  results  of  such  an  investigation*. 

The  Lockheed  SCATHA  SC3  data  base^**’*  consists  of  ~  950  days  of  dig¬ 
itized  data  on  magnetic  tape  from  launch  in  1979  to  the  present  in  several 
intervals.  The  pitch  angle  for  each  accumulation  period  is  determined  using 
data  from  the  on-board  vector  magnetometer  (SCll).  From  the  data  tapes  a 
condensed  data  base  of  10-min  averaged,  pitch-angle  binned  electron  fluxes  for 
the  24  energy  channels  over  the  full  energy  range  of  47  to  4970  keV,  has  been 
generated  for  on-line  storage  on  a  VAX  11/780  computer.  For  each  10-min  time 
interval  the  data  are  grouped  by  pitch-angle  values  into  3*  bins  from  0®  to  90®, 
making  30  pitch-angle  bins  for  each  energy  channel  in  the  interval. 

The  SCATHA  SC3  temporal  data  set  for  this  verification  was  from  the  latter 
portion  of  day  165,  June  13,  1980,  a  moderately  active  period  with  Kp  =  -(-4 
and  Dst  =  -50.  Two  days  earlier  a  geomagnetic  storm  had  occurred,  but  the 
magnetosphere  had  since  relaxed  enough  that  the  observed  electron  fluxes  on 
day  165  were  smoothly  varying  with  no  pitch-angle  signatures  of  particle  loss  in 
the  magnetopause.  Figure  5  shows  the  10-min  averaged  data  integrated  over  all 
pitch  angles  for  7  channels  of  electron  flux.  During  this  period  the  spacecraft 


Figure  6.  Typical  pitch-angle  distributions  for  the  7  energy  channels 
of  the  data  set  referred  to  in  Figure  5. 


B-10 


was  drifting  inward  from  L  =  7.2  to  L  =  5.3.  The  L-shell  range  included  in  the 
analysis  has  been  restricted  to  L  <  6.48  to  eliminate  the  large  variability  of  the 
electron  fluxes  beyond  the  geosynchronous  region^  and  the  low  count  rates  due 
to  the  approximately  exponential  fall  off  with  increasing  L-shell' .  In  addition  we 
have  selected  for  the  formal  analysis  only  those  data  channels  above  256  keV  so 
as  to  eliminate  channels  which  are  most  likely  to  be  affected  by  convective  forces 
resulting  from  consequences  of  the  major  storm  preceding  the  data  period.  Data 
with  electron  energies  above  2208  keV  have  low  statistics  and  are  also  excluded 
from  the  data  set.  Figure  6  illustrates  typical  pitch-angle  distributions  during 
this  time  for  the  selected  energies.  Although  the  period  was  moderately  active, 
the  magnetic  field  was  smoothly  varpng,  with  a  field  direction  that  remained 
perpendicular  to  the  spacecraft  spin  vector.  This  allowed  SC3  to  look  into  the 
loss  cone  for  the  whole  period,  providing  data  over  the  full  pitch-angle  range. 
In  addition,  the  pitch-angle  distributions  were  smoothly  varying  and  without 
signatures  that  indicate  loss  at  the  magnetopause. 

The  application  of  the  diffusion  theory  summarized  above  to  the  SC3  elec¬ 
tron  data  for  day  165  has  been  discussed  in  detail  in  an  earlier  paper^.  The 
verification  that  the  theory  is  a  good  description  of  the  on-trajectory  SCATHA 
data  is  indicated  in  Figure  7,  in  which  the  entire  data  set  is  shown  in  Panel(a) 
as  a  three-dimensional  plot  in  terms  of  pitch  angle  and  L{t).  Panel(b)  shows  the 
reconstructed  data  ?*t  after  it  has  been  decomposed  into  a  sum  of  pitch-angle 
eigenfunctions  truncated  at  tne  fourth  order.  Comparison  of  Panels  (a)  and  (b) 
shows  that  the  pitch-angle  eigenfunction  series  is  a  very  efficient  representation 
of  the  data  set.  In  Panel  (c)  we  also  show  the  reconstruction  of  the  data  set 
after  decomposition  in  space,  using  the  property  that  L  and  t  are  ap¬ 

proximately  decoupled  for  the  SCATHA  orbit  as  illustrated  by  Panels  (a)  and 
(b)  of  Figure  4.  In  order  to  show  the  worst  case  scenario  we  show  in  Panel  (c) 
of  Figure  7  the  reconstruction  after  the  first  iteration  of  the  decomposition;  we 
have  not  attempted  to  obtain  the  best  solution  at  this  time.  As  we  can  see  by 
a  comparison  of  the  panels  in  Figure  7,  we  conclude  that  multi-mode  diffusion 
theory  in  a  fairly  simplified  form  is  a  good  representation  of  the  electron  outer 
belt  fluxes.  This  implies  that  multi-mode  diffusion  theory  can  be  used  as  a 
fundamental  basis  of  radiation  belt  dynamic  modeling. 

Verification  that  multi-mode  diffusion  theory  is  a  good  basis  for  radiation 
belt  dynamic  modeling  involves  more  than  just  data  representation  and  the 
implied  power  of  flux  mapping  throughout  the  entire  belt,  although  the  appli- 
cational  needs  are  focused  onto  this  area.  From  the  science  view  point,  an  im¬ 
portant  part  of  dynamic  modeling  is  the  determination  of  transport  parameters 
which  are  related  to  the  microscopic  physics  of  wave-particle  interactions  in  the 
magnetospheric  plasma.  The  above  data  representation  exercise  yields  physical 
diffusion  parzuneters  with  maximum  error  spread  factors  of  2-5  that  can  be  com¬ 
pared  to  other  determinations.  The  radial  diffusion  coefficient  Du  determined 
in  our  exercise  and  shown  on  Figure  8  as  the  shaded  trapezoid  region  compares 
well  with  other  determinations.  Indeed,  our  determination  is  able  to  obtain  the 
energy  dependence  from  0.8  MeV  to  2.0  MeV.  Our  pitch-angle  diffusion  coeffi¬ 
cients  are  also  in  good  agreement  with  other  studies*.  These  give  confidence  that 
our  dynamic  modeling  exercise  is  not  in  conflict  with  independent  expectations. 
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Figure  8.  Comparison  of  the  radial  diffusion  coefficient  determined 
for  the  energy  range  0. 8-2.0  MeV  (mottled  area  show¬ 
ing  decreasing  shading  with  increasing  energy)  with  a 
previous  compilation^  which  yields  curve  1.  Curves  2- 
8  are  experimental*  Curves  9-14  are  theoretical  or 
semitheoretical'*"'*. 

FURTHER  ISSUES  IN  DYNAMIC  MODELING 

The  test  case  exercise*  summarized  above  may  have  resolved  an  important 
issue  in  radiation  belt  dynamic  modeling,  namely  that  multi-mode  diffusion 
theory  is  a  viable  basis  for  dynamic  modeling.  However,  there  are  several  other 
issues  that  must  be  resolved  before  CRRES-like  applications  can  be  made.  Here 
we  shall  discuss  what  we  preceive  to  be  the  two  most  important  issues. 
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(A)  Flux  mapping;  The  representational  use  of  dynamic  modeling  must 
clear  another  hurdle  before  it  can  be  made  operational.  The  issue  arises  be¬ 
cause  of  the  multi-dimensional  nature  of  phase  space.  Instead  of  basing  the 
discussion  on  phase  space  mapping  in  the  abstract,  we  can  illustrate  the  issue 
with  our  familiar  (£,  t)-space.  Armed  with  a  diffusion  theory  representation 
tested  on-trajectory,  as  in  the  exercise  above,  we  can  obtain  fluxes  at  any  given 
off-trajectory  point  by  “mapping”  the  on-trajectory  flux  utilizing  the  representa¬ 
tion.  Without  a  general  representation  such  as  (2),  the  mapping  is  not  unique  if 
it  is  not  based  on  a  solution  of  the  assumed  diffusion  physics.  Figure  9  illustrates 
this  mapping  based  on  the  assumed  diffusion  theory  of  the  SCATHA  test  case 
above.  The  fluxes  at  Point  2,  outside  of  the  SCATHA  trajectory,  are  uniquely 
specified  by  the  diffusion  parameters  based  on  the  fitted  fluxes  (dashed  curves) 
on  all  points  along  the  SCATHA  trajectory  between  Points  1  and  3.  In  this  re¬ 
gard,  uniqueness  is  somewhat  relative  because  it  depends  on  the  errors  involved 
in  specifying  the  on-trajectory  diffusion  peuameters.  A  case  of  extreme  constrast 
to  Figure  9  is  to  imagine  that  no  solution  of  diffusion  theory  was  available  as  the 
basis  for  dynamic  modeling.  In  such  a  case,  Point  2  can  still  be  mapped  in  (L,i) 


Figure  9.  Illustration  of  flux  mapping  using  a  nominal  dynamic 
model:  our  test  case  verified  on  the  SCATHA  trajectory. 
To  the  extent  that  the  on-trajectories  fits  are  good,  the 
mapping  bj  the  diffusion  representation  (2)  is  “unique”. 
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from  Point  1  or  Point  3  along  t  and  L  axes,  respectively,  by  some  phenomeno¬ 
logical  construct.  Indeed,  a  large  multitude  of  such  mappings  can  be  made  to 
Point  2  from  data  on  the  trajectory.  Since  these  individual  mappings  need  not 
be  consistent  with  each  other,  the  uniqueness  (within  observational  errors)  of 
the  predicted  flux  at  Point  2  will  be  lost  without  the  diffusion  representation  as 
a  consistent  basis  for  the  mapping.  Having  a  diffusion  representation,  however, 
does  not  free  us  from  the  fetters  of  non-uniqueness.  Indeed,  the  errors  incurred 
in  finding  the  on-trajectory  diffusion  parameters  operate  in  the  same  way  as 
the  inconsistent  multiple  mappings.  Similar  to  these,  the  errors  tend  to  vitiate 
the  validity  of  a  given  representation  for  dynamical  modeling.  There  are  two 
necessary  elements  in  ensuring  “uniqueness'^  m  daia  mapping.  The  first  is  to 
minimize  errors  in  diffusion  parameters  for  a  given  test  diffusion  representation. 
The  second  is  to  determine  the  most  “unique”  representation  among  severed 
theoretical  specifications  of  the  diffusion  physics.  The  first  element  is  much  eas¬ 
ier  to  carry  out  than  the  second  since  the  only  way  to  insure  that  a  diffusion 
representation  is  “unique”  is  to  show  that  it  works  in  practice.  The  theory  of 
transport  coefficients  in  the  magnetosphere  is  presently  not  advanced  enough  to 
allow  us  to  choose  the  “unique”  physics  on  an  a  priori  bzisis. 

(B)  Invariant  Representations;  An  area  of  dynamic  modeling  not  touched 
upon  above  is  the  validity  of  the  implicit  assumption  in  the  test  case  that  fluxes 
on  the  same  T-shell  are  the  same  at  the  same  time  and  energy.  We  know  from 
Figure  2  that  on  the  time  scale  of  a  day  this  is  not  true  as  effects  of  magnetic 
and  electric  shell-splitting,  local  injections  and  losses  introduce  a  local  time 
(azimuthal)  dependence  upon  the  data.  This  crucial  aspect  of  radiation  data  is 
under  current  investigation  on  a  case  by  case  basis*®;  however,  there  is  currently 
little  work  done  on  how  these  effects  can  be  incorporated  into  a  dynamic  response 
model.  A  possible  solution  to  this  important  problem  may  be  to  return  to  the 
most  fundanittii.au  representation  of  radiation  belt  diffusion  physics,  i.e.,  that 
of  diffusion  in  adiabatic  invariant  space  rather  than  in  its  various  reductions  to 
configuration  space^,  such  as  in  the  above.  The  abstraction  required  to  treat 
data  at  such  a  level  is  somewhat  of  a  drawback;  however,  a  test  case  attempt  is 
needed  to  identify  the  detauled  issues  involved. 

CONCLUSIONS 

Various  elements  required  in  radiation  belt  dynamic  modeling  have  been 
discussed.  It  is  seen  that  dynamic  modeling  differs  fundamentally  from  radiation 
belt  data  analysis  in  that  a  theoretical  representation  of  the  on-trajectory  data 
set  is  required  to  predict  fluxes  outside  of  the  satellite  orbit.  The  faithfulness  of 
the  projection  depends  on  whether  the  representation  is  a  good  description  of  the 
physical  diffusive  processes  seen  in  the  on-trajectory  data.  W’e  have  summarized 
the  results  of  a  test  case  study^  on  such  a  tractable  representation  which  yields 
results  consistent  with  basically  all  other  measurements  in  the  moderately  active 
outer  electron  radiation  belt.  The  resolution  of  a  number  of  important  issues 
brought  out  by  the  examination  of  the  test  case  awaits  further  work  before 
operational  dynamic  models  can  be  constructed. 
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Abilract.  The  interactions  of  waves  and  trapped  elec¬ 
trons  are  described  by  a  closed  system  of  three  nonlin¬ 
ear  coupled  equations.  The  model  has  applications  to 
many  aspects  of  wave  particle  interactions  in  the  magne¬ 
tosphere.  Nonlinear  numerical  solutions  pertinent  to  au¬ 
roral  pulsations  have  been  computed  for  realistic  ranges 
of  physical  parameters.  The  results  confirm  that  precipi¬ 
tation  pulsations  can  be  initiated  by  either  an  injection  of 
energetic  trapped  particles  or  an  increase  in  the  plasma 
ionization  density.  Relaxation  of  the  system  results  in 
cyclic  phase  trajectories  about  a  (quasi  equilibrium)  point 
that  is  not  located  at  the  initial  equilibrium  point.  That 
observation  may  explain  the  tendency  of  auroral  pulsa¬ 
tions  to  appear  superimposed  upon  an  increase  in  the 
general  rate  of  electron  precipitation.  The  results  also 
imply  that  precipitation  pulsations  are  most  readily  initi¬ 
ated  by  an  encounter  of  drifting  energetic  electrons  with 
a  region  of  enhanced  cold  plasma  density. 

1.  Introduction 


it  provides,  for  the  first  time,  explanations  to  a  number 
of  nonlinear  features  of  the  morningsidr  aurora. 

Nonlinear  processes  are  expected  to  govern  the  inter 
actions  of  waves  and  particles  in  the  outer  regions  of 
the  magnetosphere.  In  the  approach  adopted  here,  the 
growth  of  VLF  waves  and  the  precipitation  of  trapped 
electrons  is  described  by  a  system  of  nonlinear  differential 
equations  [Davidson  and  Chiu,  1986].  The  system  com 
prises  three  equations:  an  equalirm  descrilung  the  evolu 
lion  of  the  total  electron  flux  /  [Schulx,  1971], 

dt  1  ■+  Dt,  ^  1  -f  /for,  ’ 

an  equation  describing  the  evolution  of  the  wave  intensity, 
as  represented  by  the  diffusion  rate  D  j.Schiilz,  197T|, 


The  interaction  between  waves  and  trapped  particles 
in  the  magnetosphere  is  a  key  problem  in  magnetospheric 
physics  [Kennel  and  Petschek,  1966]  whose  basic  elements 
are  considered  to  be  well  understood.  However,  the  com¬ 
plexity  of  a  full  treatment,  based  on  first  principles,  of  the 
diffusion  of  particles  into  the  loss  cone  and  the  response 
of  waves  to  particle  diffusion  has  restricted  quanlllalive 
studies  to  a  handful  of  artificial  situations  [e  g.,  Iiian  el  al., 
1985].  Other  modelling  applications  for  which  a  full  treat¬ 
ment  is  needed  are:  the  dynamic  response  of  the  radiation 
belts  to  injections  of  cold  plasma,  kinetic  modelling  of  the 
generator  which  supplies  a  steady  stream  of  electrons  to 
the  loss  cone  opened  up  by  a  parallel  potential  drop  in  the 
eveningside  aurora,  and  modelling  of  the  morningside  au¬ 
rora.  Among  th^se,  the  morningside  aiirc-ra  problem  is  the 
best  natural  laboratory  because  it  is  generally  believed  to 
be  a  ca.se  of  pristine  wave  particle  interactions,  relatively 
unaffected  by  the  complications  of  parallel  electric  fields. 
Given  the  difhcullies  of  a  first  principles  treatment,  an 
heuristic  model,  in  which  detailed  processes  are  repre¬ 
sented  by  simple  empirical  relations  je.g.  Schulz,  197-1|, 
is  appropriate  because,  once  tested,  it  ran  then  be  pack¬ 
aged  as  a  subsidiary  unit  to  be  incorporated  into  the  more 
complex  magnetospheric  infidels  mentioned  above.  This 
paper  demonstrates  that  such  an  approach  is  feasible,  for 
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And  an  fqualion  describing  the  evi)hiliott  of  the  pil(l>- 
angle  anisotropy,  as  represented  by  tbe  wnvr  gr«>wth  rate 
7  (Davidson  and  Chiu,  1986|, 


d-)  T,  I)q  -  I)  1  (1/ 

d(  ^  r,  1  -t  Dt,  ^  7  »!/ 


(3) 


T  he  iubscripl  o  rlrn<»trs  the  vnliies  of  the  vnriablr  param¬ 
eters  al  the  cquilil)rinin  point,  whrrr  the  liin<'  drrivallvrs 
vanish.  'I  he  physical  p.irari>elerB  are  the  strmig  diffusion 
lifetime,  t,,  the  wave  loss  or  escape  li?nr,  t,,.,  the  wave 
source,  ir,  and  a  resjinnse  lime,  r, ,  that  is  a  measure 
of  the  lime  taken  ff>r  pilch  angle  diffnsinn  f*»  liil  ihr  loss 
cone.  Since  we  are  mainly  interested  in  purely  lrm|)«> 
ral  behavior,  the  relative  transverse  mf>llon»-  f*f  waves  and 
particles  are  ignorefl. 

'I  he  closure  of  the  system  through  I'ccpialion  (3)  is  of 
primary  importance  For  example,  in  the  rase  r)f  morn 
ingside  aurora,  the  first  Iwfi  equations  yield  long  prriful 
oscillatory  solutions  ff>r  large  values  r>f  but  llie  third 
equation  is  critical  to  the  formation  of  prerqdt alion  pulsa 
lions  at  the  relatively  short  periods  observed  in  pulsating 
aurorae  (Davidson  and  ('hiii,  While  (he  exlslenre 

of  short  period  pul.salif>ns  was  drmonsl rat r  tl  by  I'avnI 
son  |i986|  in  the  context  of  a  two  {liriiensitinal  system, 
the  t wo-iiimensional  morlel  offers  no  clues  to  the  pnl 

salions  might  be  initialed  *1  hr  rhararlerisl  u  periods  re 
suiting  fr<»m  a  linear  analvsiA  have  been  di‘7(nssrd  in  a 
prcvitMis  paper  (Davldsoii  and  ('Inn,  I9H(»|,  the  work  rle 
scribed  here  is  ran  exlen'non  of  lliat  stinl)  to  filjlain  fully 
nonlinear  sfilulions 
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Fig.  1.  Solutions  of  the  nonlinear  equa^ns  for  equi 
l^rium  diffusion  rates  (dimensionless)  o{  D  =  0.06  and 
D  =  0.1 ,  started  with  a  displacement  io  D  =  2Do  :  7,  D 
projections. 

In  Section  2  we  describe  the  results  of  numerical  com¬ 
putations  showing  how  the  system  responds  to  various 
assumed  initial  perturbations,  and  the  assumptions  upon 
which  the  computations  were  based.  Our  conclusions  are 
discussed  in  Section  3. 

2.  Results  Of  Numerical  Computations 
It  has  been  found  useful  to  reduce  the  number  of  pa¬ 
rameters  by  scaling  the  solutions  to  the  strong  diffusion 
lifetime  (Davidson  and  Chiu,  1986|.  The  physical  param¬ 
eters  are  thus  replaced  by  their  dimensionless  forms: 


V  H  L)tT,/(l  f  Dot,),  (-la) 

W  =  Wt,  =  D,r’(l/r„-27o)  >  0,  (46) 

7  =  7oT, ,  (4c) 

f,  =  r,/T,  .  (4d) 


The  behavior  of  a  system  of  differential  equations  de¬ 
pends  strongly  on  the  location  of  the  phase-spare  attrac¬ 
tors  or  critical  points  (Heinler  and  Orszag,  1978|.  The 
system  of  three  first  order  equations  (1 — 3)  has  three  crit¬ 
ical  points.  The  nonlinear  equations  are  similar  to  a  well 
known  set,  the  Lorenz  system  (Lorenz,  I963(  that  has  a 
strange  attractor.  It  is  therefore  worthwhile  to  inquire 
whether  any  of  the  rriliral  points  in  the  present  system 
constitute  a  strange  attractor. 

I'here  is  a  simple  criliral  point  at  the  equilibrium  val¬ 
ues  of  =  Io,  D*'*  =  Do,  and  7'''  =  70-  All  the  phase 
trajectories  eventually  spiral  inward  toward  that  point. 
The  most  interesting  solutions  in  that  region  are  those  for 
which  the  linear  analysis  yields  a  characteristic  frequency, 
II,  having  a  real  part.  Solutions  for  which  1 11 1  <t:  1  /t,  can 
be  ignored  because  they  will  not  lead  to  perceptible  ef¬ 
fects  within  realistic  time  scales.  The  linearized  solutions 
for  II  at  small  values  of  D  have  two  branches  with  purely 
imaginary  11.  1  he  upper  branch  gives  for  11  ■  Ur, 


II  =  Itrillj  ~  (5) 

the  lower  branch  gives 

n  =  ImH,  ~  I  Wf,.  (6) 

The  numerical  designation  of  the  roots  of  11  is  the  same  as 
in  Davidson  and  Chiu  (I986|,  in  order  of  increasing  |  fl  |. 

Cyclic  solutions  are  are  found  above  some  critical  dif¬ 
fusion  rate,  D  =  D* .  Above  D*  complex  roots  appear 
in  the  characteristic  equation,  whose  real  and  imaginary 
parts  are 

Dj.j  ~  ±  ^2yD/fr  +  tW/2D.  (7) 

The  traiisition  between  the  two  types  of  solutions  occurs 
where  the  four  roots  meet,  or,  more  precisely,  at 

D*  ^  TF’^’fy’/27'^V  (8) 

We  concentrate  on  parameter  regimes  appropriate  to 
the  outer  trapping  regions  and  the  morningside  auroral 
regions.  In  the  vicinity  of  D*,  the  behavior  becomes  more 
regular  toward  either  higher  or  lower  values  of  V.  For  the 
values  of  the  parameters  adopted  by  Davidson  and  Chiu 
(1986),  r,  =  240  s,  f,  =s  IO’’,  7  «  10^  F  =:  1  -  100, 
the  resulting  pulsation  period  at  D*  is  about  1.5  —  7  b, 
which  agrees  well  with  the  observed  periods  of  pulsating 
aurorae  (Johnstone,  1979,  1983]. 

One  of  the  other  critical  points  i.s  ;>l)y sically  uninter¬ 
esting,  because  it  corresponds  to  a  vanishing  total  flux, 
1  =  0  The  third  critical  point,  at  ■  Io,  D*’*  = 

Vo\2yD  -t  1],  7***  =  0  is  not  a  simple  attractor.  Al¬ 
though  there  is  a  |iossibility  that  an  intrriti  tion  between 
the  first  and  third  critical  points  could  result  in  a  strange 
attractor,  the  phase  space  trajectories  for  realistic  sim- 
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Fig.  n.  Solution  of  the  nonlinear  equalioni  for  an  equilibrium  diffusion  rale  (dimensionless)  of  D  =  0.5,  started  with 
a  (lisplaceinenl  to  7  =  I.OI70  (panels  li  and  c).  Panel  a  shows  a  case  with  the  same  physical  parameters,  but  starting 
conditions  similar  to  Figure  2.  'i  he  innermost  loops  of  the  spirals  are  omitted  for  clarity;  the  trajectories  all  eventual¬ 
ly  decay  l]a(  k  to  the  origin. 


Illations  in  the  parameter  regimes  appropriate  to  morn- 
ingside  aurorae  do  not  approach  the  third  critical  point 
closely  enough  that  they  could  he  significantly  influenced 
by  it, 

F.xamples  of  the  numerical  results  are  presented  in  F'ig- 
ures  1  through  3  I  hey  show  projections  of  the  phase  tra¬ 
jectories  on  the  (  /,  f>]  and  (  /),  7  )  planes.  ('Fhe  figures  are 
plotterl  to  dimensionless  scales  to  emphasize  the  univer¬ 
sality  of  the  solutions.)  'I  he  figures  were  constructed  for 
the  parameters  --  10"^,  7  -  lO’,  IF  10.  Those  val¬ 
ues  are  plivsirally  realistic,  and  were  chosen  because  they 
provide  the  iimst  easily  visualized  trajectories.  Figures  t, 
and  2  depict  rases  that  start  with  a  displacement  along 
the  I  axis.  An  initial  displacement  of  /  is  equivalent  to 
a  siirlden  inrre.ise  in  the  total  energetic  trapped  parti¬ 
cle  intensity.  I  he  (/,  I))  projections  are  not  presented  in 
Figure  1  because  they  are  characterized  by  a  rapid  rise 
followed  by  a  ne.irly  monotonic  decay,  similar  to  the  tra¬ 
jectory  in  Figure  2  Three  cases  are  shown  in  order  of 
increasing  vulucs  of  I),  The  trajectories  are  qualitatively 
similar  to  Figure  2  for  larger  values  of  7b  At  the  approach 
to  strong  diffusion,  /)  ~  1,  the  spirals  become  loo  tightly 
wrapped  to  rr()rodiice  on  the  scale  of  the  figures. 

f  igure  3,  p.iiK-U  b  and  c,  begins  with  an  initial  dis¬ 
placement  coiilined  to  the  (/•,  7)  plane,  f  igure  3  is  in- 
tenilerl  to  repri  sent  a  sudden  Increase  in  the  plasma  den¬ 
sity.  fhe  wave  growth  rate  Is  proportional  to  the  total 
plasma  ficiisily  IKenriel  and  I’elsrhek,  I‘)(i6|,  so  such  a 
perliirbalioii  icurrsponds  to  a  sudden  change  in  7,  from 
an  equilibriniii  slate  having  a  relatively  low  value  of  7  to 
a  new  state  having  an  Increased  value  of  7.  llecause  the 
equilibrium  solutions  of  (2)  Involve  both  D  and  7,  the 
initial  displacement  must  be  related  to  an  initial  value  of 
1)  consistent  with  the  earlier  equilibrium  state.  Thus,  for 
a  constant  It  ,  the  value  of  Pg  in  eai  h  slate  must  be  in¬ 
versely  prop* .r  1 1-  ui .il  to  (270  1  / c,,, )  for  comparison  with 


the  results  of  Figures  1  and  2,  panel  3a  shows  the  phase 
trajectory  that  would  result  from  an  initial  displacement 
in  /  for  the  same  physical  parameters  as  panels  b  and 
c.  Combinations  of  the  two  types  of  initial  perturbation 
yield  trajectories  that  resemble  the  (I,  D)  projection  of 
Figure  2  and  the  (/),  7)  projection  of  Figure  3. 

3.  Conclusions 

The  results  of  numerical  solutions  of  the  nonlinear  dif¬ 
ferential  equations,  (1  —  3),  indicate  general  classes  of 
behavior  of  waves  and  trapped  electrons  expected  in  the 
auroral  regions  and  outer  trapping  regions.  They  confirm 
that  cyclic  solutions  exist,  and  that  cyclic  solutions  may 
follow  either  an  increase  in  the  trapped  electron  density 
or  an  increase  in  the  cold  plasma  density.  Such  solutions 
reproduce  the  gross  features  of  auroral  pulsations  inso¬ 
far  as  the  electron  precipitation  rate  is  a  reflection  of  the 
diffusion  rale,  D. 

Several  interesting  aspects  of  the  detailed  solutions 
could  not  have  been  predicted  from  the  linear  analysis. 
For  diffusion  rales  higher  than  D*  an  initial  perturbation 
in  the  total  trapped  flux  is  followed  by  spiral  orbits  about 
a  point  that  is  displaced  from  (/o,  /lo,  7o)-  The  periods 
of  those  orbits  are  similar  In  the  periods  derived  from  the 
linear  analysis;  the  periods  do  not  strongly  depend  on  the 
amplitude  of  the  oscillations.  The  examples  of  Figure  1, 
with  few  or  no  complete  pulsation  cycles,  are  probably  ap¬ 
propriate  to  the  radiation  bells,  where  isolated  precipita¬ 
tion  "luirsts"  are  occasionally  seen  |  rrcfall  and  Williams, 
19791.  As  the  diffusion  rate  is  further  Increased,  in  Fig¬ 
ure  2,  the  damping  rate  is  diminished  and  the  number 
of  loops  incresises.  The  spiral  phase  trajectories  do  not 
enclose  the  stable  critical  point,  but  lend  to  approach  a 
point  that  is  located  at  larger  values  of  D  and  7.  The  ef 
feci  of  the  displacement  of  the  “center"  of  the  spirals  is  to 
shift  the  average  diffusion  rale  so  that  the  instantaneous 
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diffusion  rate  remains  greater  than  the  equilibrium  value. 
For  morningside  aurora  that  result  is  consistent  with  the 
observation  by  Davidson  and  Sears  |19B0|  that  auroral 
pulsations  are  superimposed  upon  a  general  elevation  of 
the  emission  intensity  resulting  from  electron  precipita¬ 
tion;  the  “background”  intensity  tends  to  rise  before  the 
pulsations  begin.  Such  a  finding  is  important,  because  it 
confirms  the  need  for  a  three-parameter  system  (Equa¬ 
tions  1  —  3),  explicitly  including  the  total  trapped  flux 
as  a  variable. 

A  comparison  of  Figures  3a  and  3c  reveals  that  the 
amplitude  of  the  pulsation  amplitude,  as  measured  by 
the  precipitation  rate  or  D,  is  most  sensitve  to  the  to¬ 
tal  plasma  density.  While  it  is  difficult  to  speculate  on 
what  constitutes  a  “typical”  initial  perturbation,  the  sim¬ 
ulations  in  which  the  initial  total  flux  is  increased  by  a 
factor  of  2  seem  reasonable.  The  initial  perturbation  of 
Figure  3,  on  the  other  hand,  corresponds  to  a  mere  1% 
increase  in  the  plasma  density  or  7.  The  implication,  ap¬ 
plicable  to  a  broad  range  of  physical  parameters,  is  that 
an  abrupt  increase  in  the  plasma  density  experienced  by 
drifting  electrons  is  more  effective  in  initiating  precipita¬ 
tion  pulsations  than  an  increase  in  the  total  intensity  of 
energetic  electrons. 

The  sensitivity  of  the  pulsation  amplitude  to  a  sud¬ 
den  increase  in  70  raises  an  interesting  question:  whether 
a  compression  of  the  magnetic  field  might  initiate  pulsa¬ 
tions,  as  suggested  by  Coroniti  and  Kennel  |1970|.  While 
a  detailed  treatment  of  that  issue  is  beyond  the  scope 
of  the  present  work,  the  results  discussed  above  indicate 
that  such  a  mechanism  might  be  very  effective.  Any  per- 
tubation  in  the  initial  value  of  70  would  serve  to  initiate 
pulsations. 

The  picture  implied  by  the  above  results  for  pulsat¬ 
ing  aurorae,  is  that  they  follow  an  injection  of  energetic 
trapped  electrons  in  the  outer  trapping  regions.  Repet¬ 
itive  pulsations  turn  on  when  the  drifting  trapped  elec¬ 
trons  encounter  a  region  of  enhanced  cold  plasma  density. 
These  conclusions  follow  from  a  closed  nonlinear  model 
that  dues  not  depend  on  the  specification  of  wave  modes; 
predictions  from  the  model  are  fully  general  and  can  be 
applied  to  the  various  cases  mentioned  in  the  introduc¬ 
tion.  The  nature  of  the  detailed  wave  modes  enters  only 
in  the  determination  of  the  physical  parameters  (4). 

A  strange  attractor  was  not  found  for  simulations  in 
the  time  domain  alone.  We  leave  open  the  question 
whether  a  strange  attractor  may  occur  for  conditions  ap¬ 
propriate  to  spatial-temporal  structures.  Although  our 
model  explains  the  semi-regular  temporal  structure  of 
pulsating  aurorae,  the  apparently  random  spatial  struc¬ 
ture  remains  a  major  unsolved  problem. 
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Towards  Dynamic  Modeling  of  the  Outer  Electron  Radiation  Belt 
Y.  T.  Chiu,  M.  A.  Rinaldi,  W.  E.  Francis  and  R.  W.  Nightingale 
Lockheed  Palo  Alto  Research  Laboratory.  California 

Abstract 

We  show  that  the  general  solution  of  a  simple  form  of  the  simultaneous  bimodal  (ra¬ 
dial  and  pitch-angle)  diffusion  equation  can  represent  the  SCATHA  outer  belt  electron 
distribution  to  a  high  degree  of  accuracy  for  periods  up  to  more  than  10  hours.  The 
data  representation  requires  diffusion  parameters  that  are  in  agreement  with  all  previous 
experimental  and  theoretical  determinations  of  such  diffusion  coefficients.  Such  a  represen¬ 
tation  satisfies  the  basic  requirements  towards  the  dynamic  modeling  of  the  outer  electron 
radiation  belts  for  Z  <  7  at  quiet  and  moderate  geomagnetic  conditions.  For  Z  ~  7 
the  representation  has  difficulties  with  “butterfly”  distributions  that  may  signify  energetic 
electron  encounters  with  the  magnetopause.  It  is  demonstrated  with  a  typical  example  of 
computer  simulation  of  electron  encounters  with  the  magnetopause  that  the  effects  of  the 
encounter  may  be  highly  anisotropic  in  configuration  or  in  invariant  space. 

1.  Introduction 

The  variability  of  the  electron  fluxes  in  the  outer  radiation  belt  (Z  >  5)  is  not  only 
a  matter  of  concern  for  space  activities  in  the  region  but  is  also  a  major  challenge  to 
understanding  the  macroscopic  consequences  of  wave-particle  interactions  in  the  magne¬ 
tosphere.  Theories  on  the  generation  of  plasma  wave  disturbances  and  on  the  scattering 
of  individual  electrons  by  such  waves  exist;  however,  a  systematic  test  of  the  integration 
of  such  microscopic  effects  into  a  diffusion  description  of  the  globail  dynamics  of  the  outer 
electron  radiation  belt  has  yet  to  be  attempted.  To  be  sure,  episodic  w’ork  je.g.,  Reagan  et 
al..  1981;  West  et  al..  1981;  Sibeck  et  al.,  1987]  have  established  the  plausibility  of  various 
aspects  of  the  diffusion  description  in  some  regions  of  space.  On  the  global  scale,  even  a  de¬ 
scription  of  the  observationsd  data  organized  according  lo  the  generic  concepts  of  diffusion 
theory  does  not  presently  exist.  In  anticipation  to  the  upcoming  .A-ir  Force-NASA  Com¬ 
bined  Release  and  Radiation  Effects  Satellite  (CRRES),  we  have  initialed  an  attempt  lo 
determine  if  diffusion  theory  can  provide  an  adequately  quantitative  and  computationally 
efficient  description  of  the  behavior  of  the  outer  electron  belt  as  observed  by  the  SCATH.A 
satellite  Chiu  et  al..  1988;  hereafter  referred  lo  as  Paper  L.  This  attempt  was  not  directed 
towards  the  microphysics  of  wave-particle  interaction  theory;  rather,  it  addressed  a  much 
more  primitive  question;  Were  the  observational  data  organizable  in  terms  of  simultaneous 
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pitch-angle  and  radial  diffusion  theory,  even  in  the  restrictive  context  of  piece-wise  slices 
of  space-time  as  in  satellite  observations? 

This  work  and  the  preceeding  effort  in  Paper  I  are  focused  upon  this  central  issue 
of  whether  a  computationally  efficient  dynamic  model  of  the  outer  electron  belt  can  be 
constructed  in  principle,  or  not.  The  resvilts  of  Paper  I  were  encouraging  in  that  solutions 
of  simult2meous  bimodaJ  (pitch-angle  and  radial)  diffusion  theory  were  shown  to  have  the 
potential  to  provide  an  efficient  description  of  a  SCATHA  outer-belt  electron  data  set 
over  the  phase  space  volume:  (5.3  <  L  <  6.7;  At  =  5  hrs.;  200  keV  <  E  <  2000  keV; 
Qc  <  a  <  90°),  where  a  is  the  pitch-angle,  is  the  edge  of  the  loss  cone  and  E  is  the 
energy.  The  description  was  achieved  in  essentially  the  initial  pass  without  a  systematic 
fitting  procedure.  The  diffusion  parameters  required  happened  to  cluster  around  values  in 
agreement  with  those  obtained  by  other  workers.  These  initial  results  established  plausibity 
for  the  approach.  However,  to  go  beyond  plausibility  towards  dynamic  modeling  of  the 
outer  electron  belt,  a  number  of  deficiencies  of  Paper  1  need  to  be  addressed.  Among  these, 
at  least  three  major  issues  require  attention  if  the  bimodal  diffusion  approach  to  dynamic 
modeling  is  to  be  pursued. 

First,  we  must  establish  a  systematic  iterative  procedure  to  achieve  a  fit  of  the  data 
consistent  with  all  the  constraints  of  the  solution.  This  requirement  involves  not  only 
specification  of  the  accuracy  of  the  fit  in  four-dimensional  phase  space  but  also  satisfaction 
of  diffusion  parameter  constraints  imposed  by  bimodal  theory.  Second,  the  form  of  the 
bimodsd  diffusion  theory  used  was  computationally  efficient  but  its  validity  was  recognized 
to  be  restricted  to  a  sufficiently  thin  slice  of  i-space  [Walt,  1970j.  Such  a  property  is  not 
incompatible  with  dynamic  modeling  with  a  satellite  data  base  since  satellite  data  are  nec¬ 
essarily  restricted  in  {L,t)  space.  Nevertheless,  it  is  crucial  to  determine  the  observational 
extent  of  L-space  over  which  the  version  of  bimodaJ  diffusion  theory  is  valid.  Third  and 
most  important,  the  degree  of  uniqueness  of  the  data  representation  must  be  established 
since  dynamic  modeling  implies  the  projection  of  electron  fluxes  from  measurements  on 
a  satellite  trajectory  to  other  regions  of  phase  space  where  direct  measurements  are  not 
available.  Note  that  the  solution  representation  of  the  theory  is  necessarily  unique  every¬ 
where  if  the  on-trajectory  specification  is  unique.  However,  the  data  fit  to  the  solution  is 
over  a  trajectory  on  which  L  and  i  are  mixed;  therefore,  the  on-trajectory  specification 
is  not  necessarily  unique  insofar  as  the  separation  of  L  and  i  effects  are  concerned.  In 
order  to  insure  that  our  iterative  procedure  for  untangling  L  and  t  effects  does  not  alias 
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the  projection  to  other  regions  of  space-lime,  we  must  at  least  demonstrate  that  our  data 
representation  yields  flux  projections  independent  of  the  projection  path. 

In  this  paper,  we  shall  address  the  above  three  issues.  VVe  shall  demonstrate  that 
simultaneous  bimodal  diffusion  theory  is  an  efficient  representation  of  SCATHA  outer-belt 
electron  data  in  the  configuration  space  range  (5  <  i  <  7;  At  <  10  hrs.)  for  Day  165, 
1980  and  Day  110,  1979.  Iterative  procedures  have  been  devised  to  improve  the  fit  and  to 
satisfy  diffusion  parameter  constraints  that  were  not  enforced  in  Paper  I.  The  result  is  a 
reduction  of  parameters  (by  ~  15%)  to  achieve  better  fits  over  a  much  larger  data  base  (by 
~  100%).  We  are  also  able  to  demonstrate  that  fit  errors  do  not  impact  the  uniqueness  of 
off-trajectory  flux  projections  based  on  the  solution  representation. 

2.  Constraints  on  Bimodal  Diffusion  Representation 

In  order  to  organize  the  discussion  in  this  paper,  we  briefly  summarize  the  salient 
features  of  simultaneous  bimodcd  diffusion  theory  [Walt,  1970;  Falthammar  and  Walt,  1969; 
Roederer  and  Schulz,  1969  and  its  specific  form  developed  for  efficient  data  representation 
in  Paper  I. 


By  introducing  an  approximate  invariant  valid  over  a  sufficiently  thin  slice  of  L-space, 
Walt  T970  showed  that  the  bimodal  diffusion  equation  can  be  formulated  in  configuration 
and  pitch-angle  coordinates  where  x  =  cos  a.  Such  a  formulation  is  particularly 

convenient  for  data  applications  because  the  coordinates  are  simple,  although  it  needs  to 
be  emphasized  that  its  validity  is  restricted  to  slices  in  T-space,  over  which  energy  geuned 
or  lost  from  radial  diffusion  does  not  impact  Walt’s  invariant.  In  Paper  I,  we  showed  that  if 
the  pitch-angle  and  radial  diffusion  coefficients,  Du  and  Dni  respectively,  lake  on  certain 
generic  forms 


=  DUE) 

D^t  =  i{E)L- 


(n 


where  p  is  a  given  constant,  then  the  general  solution  to  the  bimoded  diffusion  equation  can 
be  expanded  in  terms  of  eigenfunctions  of  various  Bessel  forms.  This  convenient  property 
forms  our  basic  principle  to  characterize  the  state  of  the  outer  electron  belt  in  a  piecewise- 
continuous  fashion  in  terms  of  a  data  base  measured  along  a  satellite  trajectory  L{i). 


The  general  solution  consists  of  the  linear  sum  of  three  modal  expansions:  (1)  single¬ 
mode  pilch-angle  diffusion.  (2)  single-mode  radial  diffusion  and  (3)  simultaneous  bi-modaJ 
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diffusion.  Using  the  device  of  the  Kronecker  delta  6„oi  the  three  modal  expansions  can  be 
written  in  terms  of  a  single  eigenfunction  expansion 


/  =  E 


^nO  +  (1  —  Sno) 


Jo(knX/Xc) 

JlM 


(cn(l-SnOy^^” 


where  k„  is  the  zero  of  the  Bessel  function  Jo  and  Xc  = 
functions  are  defined  as  follows: 


+  e‘/‘"  la„y„(l)  ^  l>nZ„(l):) 

(2) 

cosoc-  The  parameters  and 


e  =  3(1  -  ^)/(m  -  2) 


Vn(l)  =  X"  •  [Jl^.(|  I  /!")]>  •  [Jl..(i  ;  /X")]< 

Zn(l)  =  X"  ■  /?„  i  /X^)]^  •  [XTi+.d  /3„  i  /L^)]^ 


(3) 

(4) 

(5) 


J  is  the  regvdar  Bessel  function,  I  and  K  are  the  modified  Bessel  functions,  and  N  is  the 
Neumann  function.  The  usual  notation  is  defined  in  this  case  by 


[5)>  =  g  if  <n  >  0 


(6) 


[5]>  =1  if  fn  <  0 

and  analogously  for  [5j<.  The  separation  constants  {in,0n)  and  parameters  of  the  solution 
are  related  by  the  algebraic  relations 


={1  -  <5„o)r-^  -  (7) 

=  -D..kl/xl  (8) 

The  parameter  value  A  =  4  (i.e.,  /i  =  10)  is  chosen  to  conform  with  the  consensual  X- 
dependence  ol  the  radial  ditlusion  coefficient  [Schulz  and  Lanzerotti,  1974  .  From  Paper  1 
the  parameter  i/  is  set  to  2  (i.e.,  t  =  0)  for  ease  of  computing  the  Bessel  functions.  Other 
values  of  i/  are  easily  accommodated  at  the  expense  of  computational  speed.  Our  work  so 
far  has  demonstrated  the  economy  of  using  =  2,  but  further  refinements  can  be  readily 
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accommodated  in  our  procedure.  Equation  (7)  will  be  referred  to  as  the  separation  equation 
because  the  entire  separable  diffusion  differential  equation  is  reduced  by  the  eigenfunction 
expansion  into  this  constraint  on  the  diffusion  parameters.  It  is  seen  from  (8)  that  the 
single  pitch-angle  mode  time  constants  are  negative  as  expected,  since  single-mode  pitch- 
angle  diffusion  aJways  leads  to  decay.  For  radial  diffusion,  however,  the  time  constants 
in  can  be  either  positive  or  negative  over  a  finite  time  intervaJ.  This  is  due  to  the  fact 
that  the  phase  space  density  in  the  finite  L  domain  can  rise  or  fall  in  time  in  response 
to  particle  transport  into  and  out  of  the  adjacent  volumes.  With  the  sign  of  each  in  not 
predetermined,  the  separation  constant  0^  can  be  positive  or  negative.  This  is  the  origin 
of  the  {J.I)  ajid  {N.K)  options  of  (4)  and  (5). 

When  a  data  set  is  fitted  to  the  generaJ  solution  expansion  exclusive  of  the  separation 
equation  (7),  the  electron  distribution  function  is  characterized  by  the  set  of  parameters 
(an-bn-  Cn-in-'^'ni^n)-  in  Paper  1  it  was  found  that  five  orders  o^  circnfunctions  (n  =0-4) 
per  energy  channel,  i.e.,  30  parameters  per  energy  channel,  are  sufficient  to  characterize 
a  data  base  of  5040  data  points.  For  this  paper,  we  include  the  effects  of  the  separation 
equation  (7)  in  constraining  the  diffusion  parameters  Note  that  the  diffusion 

parameter  (  is  independent  of  order  n,  whereas  the  others  are  order  dependent.  In  Paper  1, 
r„  and  3n  were  independently  derived  from  the  fit  and  (7)  was  not  imposed.  No  iterations 
were  made  and  a  value  of  (  was  determined  for  each  order  from  (7),  i.e.  a  set  of  It  was 
found  that  these  order-dependent  values  of  ^  clustered  closely  around  those  determined 
by  other  workers,  an  indication  that  the  representation  is  reasonable.  Here  we  use  the 
constraint  (7)  as  the  focus  for  an  iterative  fitting  scheme  in  determining  the  best  order- 
independent  value  of  Since  (7)  is  imposed,  and  are  no  longer  independent, 

thus  the  number  of  required  parameters  per  energy  channel  is  reduced  by  4  in  this  work. 
The  iteration  also  improves  the  fit  to  the  data  set. 

3.  Iterative  Fit  of  Data  Sets 

The  process  of  fitting  a  data  set  to  the  representation  begins  with  decomposing  the 
electron  distribution  functions  into  pitch-angle  eigenfunctions  and  their  amplitudes.  We 
have  demonstrated  in  Paper  I  that  this  procedure  is  highly  accurate.  We  shall  focus 
our  discussion  upon  the  decomposition  of  the  pitch-angle  amplitudes  into  eigenfunction 
expansions  in  L. 

We  have  devised  a  nonlinear  fitting  procedure  keyed  to  the  iterative  determination  of 
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the  order-independent  radiaJ  diffusion  parameter  (  in  the  separation  equation  (7).  In  the 
first  pass,  the  procedure  described  in  Paper  I  is  applied  and  the  separation  equation  (7) 
is  solved  for  each  order,  resulting  in  a  set  of  order-dependent  parameters  Iterations 
begin  with  the  choice  of  an  interim  mean  value  (  =  (^n)  where  the  averaging  is  weighed  by 
the  eigenfunction  amplitudes  {an,  bn).  The  interim  value  (  is  enforced  in  (7)  to  determine 
new  values  of  the  parameters  (T„,y9„)  in  a  new  fit.  With  the  new  parameters,  a  new  set  of 
^n  is  determined  by  solving  (7)  exactly  and  ajtother  iteration  cycle  begins.  The  iterative 
procedure  continues  until  the  residual  error  in  the  equality  between  the  left  and  right  hand 
side  of  (7)  can  no  longer  be  improved  upon.  The  finaJ  value  of  (  is  taken  to  be  the  fit  value 
of  the  diffusion  parameter.  This  is  a  nonlinear  iterative  procedure;  therefore,  there  is  no 
a  priori  guarantee  that  it  is  convergent.  In  practice,  it  converges  because  the  set  of  in 
the  first  pass  is  closely  clustered  to  begin  with.  Figure  1  shows  an  example  of  the  degree 
to  which  the  constraint  (7)  is  obeyed  by  the  initial  and  final  iterations  of  the  procedure. 

The  iterative  scheme  has  been  applied  to  ~  4  hours  of  SCATHA  electron  belt  data 
for  Day  165,  1980  and  to  ~  9  hours  of  data  for  Day  110,  1979.  Figures  2  and  3  show  three- 
dimensional  plots  of  the  data  and  the  fit  reconstructions  for  these  two  days,  respectively.  It 
should  be  noted  that  reconstruction  errors  for  both  days  are  primarily  concentrated  in  the 
high-£  and  a  ~  90°  region.  We  found  that  this  region  is  rich  in  “butterfly”  distributions 
[Sibeck  et  al.,  1987].  These  wiU  be  discussed  in  more  detail  in  Section  5  below.  Figures 
4  and  5  show  the  distributions  of  the  fit  errors  for  each  data  point  in  the  initial  and  final 
passes  of  the  iterations  for  the  two  days,  respectively.  The  data  points  associated  with  the 
large-error  tail  of  the  error  distributions,  less  than  1000  data  points  in  each  data  set,  are 
again  related  to  “butterfly”  distributions,  although  a  substantial  portion  of  the  large  fit 
errors  are  also  contributed  to  by  some  noisy  data  in  the  data  base.  Unsmoothed  raw  data 
is  used  to  tBilly  the  percentage  error  in  these  figures.  Figure  6  shows  the  comparison  of 
radial  diffusion  coefficients  determined  in  this  work  with  all  other  known  observational  and 
theoretical  determinations.  The  mottled  band  of  Dll  values  in  the  figure  corresponds  to 
this  work.  The  spread  of  the  band  is  due  to  energy  dependence  of  the  diffusion  coefficient. 
The  larger  cross-hatched  band  of  Dll  values  was  determined  in  Paper  I  and  is  due  not 
only  to  the  spread  of  energy  channels  but  also  primarily  to  the  dependence  of  ^  on  order. 
It  is  seen  that  our  enforcement  of  (7)  allows  for  a  much  more  precise  determination  of  the 
radial  diffusion  coefficient.  Energy  dependence  is  not  always  addressed  in  other  work. 

4.  OfT-Trajectory  Projection  of  Electron  Distributions 
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For  a  data  repesentation  scheme  to  qualify  as  a  dynamic  modeling  scheme,  it  must 
demonstrate  the  ability  to  uniquely  predict  (or  obtain  by  projection),  the  off-trajectory 
electron  distributions,  from  a  set  of  on-trajectory  electron  distributions.  The  use  of  the 
general  solution  to  bimodal  diffusion  theory  yields  a  naturally  unique  off-trajectory  projec¬ 
tion  wiihin  the  limitations  of  the  assumed  physics.  Whether  the  assumed  physics  correctly 
describe  the  actuaJ  processes  prevailing  in  the  particular  region  of  the  outer  electron  belt 
is  a  separate  issue  to  be  discussed  in  Section  5.  Within  the  scope  of  diffusion  theory,  if 
the  on-trajectory  representation  is  exact,  the  projection  to  any  off-trajectory  point  is  also 
independent  of  path  because  the  solution  is  uniquely  specified  everywhere  that  it  is  valid. 
Since  the  on-trajectory  representation  is  not  given  a  priori  (but  given  as  a  fit  based  on 
a  data  set  that  mixes  L  and  f),  an  important  question  arises:  Can  the  fit  errors  due  to 
aliaising  L  and  t  destroy  the  uniqueness  of  the  projection  by  making  it  path  dependent? 

We  have  investigated  this  problem  of  projection  consistency  (or  path  dependence) 
for  our  representation  by  projecting  to  an  off-trajectory  point  from  two  well-separated 
on-trajectory  points:  one  having  the  same  L  and  the  other  having  the  same  t  as  the  off- 
trajectory  point.  The  projection  paths  are  thus  cdong  the  t  and  the  L  axes,  respectively. 
Figure  7  shows  the  results  of  such  a  imiqueness  test  for  Day  110,  1979.  The  projection 
in  our  representation  is  indeed  path  independent.  This  test  clears  a  hurdle  for  using  the 
bimodal  diffusion  solution  for  outer-belt  dynamic  modeling,  since  it  demonstrates  that 
uniqueness  is  not  destroyed  by  the  fit  errors.  It  must  be  noted  that  this  is  only  an  internal 
consistency  test,  not  a  validity  test  that  diffusion  theory  correctly  describes  the  prevailing 
physics  of  the  outer  belt.  A  validity  test  requires  data  sets  from  more  than  one  satellite. 
A  vadid  dynamic  model  should  be  able  to  “predict”  one  satellite  data  set  from  the  solution 
fit  to  the  other  satellite’s  data  set,  and  vice  versa. 

5,  Effects  of  Drastic  Violations  of  the  Third  Invariant 

We  have  repeatedly  stated  from  the  outset  that  the  validity  of  the  forn.  of  simultaneous 
bimodal  diffusion  theory  chosen  for  this  study  is  very  severely  limited,  particularly  in 
regions  of  the  outer  electron  belt  at  T  >  7.  The  first  limitation  is  that  Walt’s  invariant 
fl970]  remains  so  only  within  a  limited  slice  of  L.  This  property  of  piecewise  invariance  is 
exploited  as  much  as  possible  in  our  application  in  the  region  5  <  L  <  7,  but  there  must  be 
constraints  to  the  exploitation.  In  the  stable  regions  of  the  outer  belt  at  T  <  5.  we  do  not 
see  any  reason  why  the  application  of  the  representation  will  be  different  from  that  shown 
here,  although  it  is  necessary  to  test  it  with  another  satellite  data  set  because  SCATHA 
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does  not  cover  that  region.  For  X  >  6,  we  sometimes  encounter  difficulties  with  “butterfly” 
pitch-angle  distributions.  We  recognized  at  the  outset  that  an  important  factor  limiting 
the  validity  of  the  representation  is  the  effects  of  third  invariant  violations  that  may  or 
may  not  be  describable  in  terms  of  diffusion.  Such  effects  can  be  due  to  severe  magnetic 
or  electric  shell  splitting  or  to  encounters  with  the  magnetopause.  In  such  encounters, 
the  electron  distribution  on  a  given  X-shell  can  be  changed  catastrophically.  Since  such 
third-invariant  violating  effects  are  deemed  to  be  most  effective  at  the  equatori2J  region,  it 
is  usually  assumed  that  their  signature  is  a  loss  to  the  drift  shell  of  equatorially  mirroring 
particles  (those  with  ~  90°  pitch  angle)  leading  to  the  so-called  “butterfly”  distributions 
[West,  1979:  Sibeck  et  al.,  1987].  Such  an  expectation  seems  reasonable  as  far  as  the 
pitch-angle  distribution  is  concerned;  however,  the  transport  in  the  configuration  part  of 
phase  space  is  a  matter  that  has  not  been  investigated  beyond  the  expectation  that  the 
third  invariant  is  drasticaUy  violated. 

Since  our  representation  is  not  expected  to  be  valid  for  describing  such  drastic  viola¬ 
tions  of  the  third  invariant,  the  inclusion  of  such  distribution  functions  in  our  data  base 
incurs  large  errors  in  the  representation.  These  constitute  a  part  of  the  large-error  tail  of 
the  error  distributions  shown  in  Figures  4  and  5.  The  successful  treatment  of  the  “butter¬ 
fly”  events  constitutes  the  remaining  hurdle  towards  dynamic  outer  electron  belt  modeling. 
The  solution  to  this  important  problem  is  beyond  the  scope  of  this  brief  report  but  we 
shall  address  some  salient  points  here.  An  initial  reaction  to  the  failure  of  the  simple  form 
of  radial  diffuion  used  is  to  assume  that  another  form  of  radial  diffusion,  perhaps  in  third 
invariant  coordinates,  may  be  a  better  description  of  the  physics  involved  in  “butterflies”. 
Such  may  indeed  be  the  case,  but  a  more  primitive  question  needs  to  be  considered  first: 
Are  encounters  involving  drastic  violations  of  the  third  invariant  necessarily  diffusive  in 
the  configuration  part  of  phase  space? 

In  order  to  answer  this  and  related  questions,  we  have  initiated  computer  simulations 
of  energetic  electron  encounters  with  magnetopauses  in  realistic  geomagnetic  field  models. 
Figure  8  shows  the  encounter  of  a  10  MeV  electron  with  a  non-turbulent  paraboloidal 
magnetopause  [Stern,  1985]  in  the  configuration  of  a  Tsyganenko  1987  field  model.  The 
magnetic  field  magnitude  drops  by  30  nT  across  the  magnetopause  at  the  sub-solar  point. 
The  field  configuration  outside  merges  smoothly  into  a  uniform  southward  interplanetary 
configuration.  Figure  8  shows  the  projection  of  the  three-dimensional  electron  trajectory 
onto  the  equatorial  plane.  The  important  feature  of  this  magnetopause  encounter  is  that 
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the  electron  executes  trapped  bounce  motion  occasionally  on  inward  field  lines  and,  in 
between  bounces,  it  skips  along  the  equatorial  regions  sometimes  inward  and  sometimes 
outward  of  the  magnetopause  boundary;  i.e.,  there  is  little  evidence  of  scatter  perpendicu¬ 
lar  to  the  boundary.  The  electron  trajectory  shown  on  Figure  8  is  one  example  in  a  Monte 
Carlo  simulation  of  encounters  with  a  smooth  magnetopause  in  which  all  scatterings  seem 
to  be  directed  along  the  boundary,  even  though  individual  trajectories  differ  in  the  amount 
of  residence  time  inside  or  outside  of  the  boundary  surface,  depending  on  their  pitch  an¬ 
gle  and  gyration  phases  as  they  encounter  the  boundary.  Such  anisotropic  scatterings  in 
encounters  with  the  magnetopause  can  hardly  be  described  in  simple  terms  of  isotropic 
diffusion,  in  invariant  coordinates  or  otherwise.  The  encounter  effects  are  somewhat  more 
isotropic  in  configuration  space  if  the  magnetopause  is  not  smooth,  but  the  anisotropy  will 
remain  since  the  average  magnetic  field  is  non-zero  both  inside  and  outside.  Clearly  the 
equatorial  skippings  are  likely  to  lead  to  enhanced  pitch-angle  diffusion;  however,  pitch- 
angle  diffusion  cannot  produce  a  “butterfly”  population.  On  the  other  hand,  skippings 
along  the  equatorial  magnetopause,  which  cuts  through  many  L-sheils,  clearly  constitute 
transport  in  L  and,  consequently,  violation  of  the  third  invariant,  as  expected.  The  loss 
to  a  given  L-shell  of  electrons,  which  skip  to  a  different  L-sheU  along  the  equatorial  mag¬ 
netopause,  will  clearly  produce  a  “butterfly”  distribution.  The  surprising  feature  of  the 
simulation  is  that  this  transport  in  L  is  anisotropic  in  configuration  space  and  can  hardly 
be  considered  diffusive.  Based  on  computer  simulations  of  magnetopause  encounters,  the 
“butterfly”  distribution  is  but  a  superficicd  signature  peculiar  to  the  process  in  which  vi¬ 
olations  of  invariants  occur  with  drastic  anisotropy  in  configuration  space.  Further,  the 
channeling  of  electron  motion  provided  by  the  sharp  magnetic  field  gradient  at  the  magne¬ 
topause  boundary  surface  indicates  that  it  acts  macroscopically  as  a  scattering  boundary 
to  outer-belt  electrons,  having  simultaneously  the  properties  of  a  sink  to  electrons  ap¬ 
proaching  the  boundary  on  given  initial  T-sheUs  intercepting  the  magnetopause  and  of  a 
highly  anisotropic  source  of  electrons  to  L-shells  intercepting  the  magnetopause  inward  of 
the  initial  L-shell.  In  the  non-turbulent  magnetopause  case,  virtually  no  electrons  were 
lost  to  interplanetary  space.  It  appears  that  diffusion  in  invariant  space  is  not  the  panacea 
that  cursory  considerations  make  it  out  to  be.  Detailed  results  of  the  computer  simulation 
study  will  be  published  elsewhere. 

6.  Conclusions 

By  testing  against  a  SCATHA  electron  data  base,  we  have  found  a  computationally 


practicable  representation  for  the  outer-belt  electron  distribution  from  the  general  solution 
of  a  simple  version  of  the  simultaneous  bimodal  diffusion  theory  that  satisfies  the  basic  re¬ 
quirements  for  outer-belt  dynamic  modehng  inwards  of  Z  ~  7.  The  representation  fails  to 
characterize  the  small  portion  of  outer-belt  electron  populations  known  as  “butterfly”  dis¬ 
tributions.  The  nature  of  outer-radiation-belt  electron  encounters  with  the  magnetopause 
needs  to  be  investigated  beyond  the  superficial  formation  of  “butterfly”  pitch-angle  dis¬ 
tributions.  In  particular,  the  distribution  of  “butterflies”  in  configuration  space  can  shed 
light  on  the  scattering  processes  m  encounters  with  the  magnetopause. 
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Figure  Captions 

Figure  1.  Values  of  the  right  hand  side  (RHS)  and  left  hand  side  (LHS)  of  the  separation 
equation  (7)  are  plotted  as  functions  of  order  n  for  (a)  the  initial  iteration  and 
(b)  the  final  iteration.  The  equality  is  enforced  cis  a  fit  constraint  in  this  paper. 
Note  the  large  divergence  at  order  2  in  the  initial  iteration. 

Figure  2.  Three-dimensional  presentation  of  the  phase  space  density  of  (a)  the  data  base 
and  (b)  its  reconstitutions  for  the  288.5  keV  channel  of  Day  165,  1980  in  pitch 
angle  and  L.  Large  errors  are  found  primarily  in  the  comer  of  (b)  pertaining  to 
large  L  and  large  pitch  angles.  Noisy  data  near  L  ~  5.7  also  resulted  in  some 
error. 

Figure  3.  Three-dimensional  presentation  of  the  phase  space  density  of  (a)  the  data  bast 
and  (b)  its  reconstitutions  for  the  448  keV  channel  of  Day  110.  1979  in  pitch 
angle  and  L.  As  in  Figure  2,  the  large  errors  are  mainly  confined  to  the  corner 
of  (b)  pertaining  to  large  L  and  large  pitch  angles. 

Figure  4.  The  distribution  of  fit  percentage  errors  in  the  reconstitution  of  all  the  data 
points  in  the  Day  165,  1980  data  base.  The  error  distributions  of  the  initial 
and  final  iterations  are  shown  as  dots  and  asterisks,  respectively. 

Figure  5.  Same  as  Figure  4  except  the  distribution  is  for  the  Day  110,  1979  data  base. 

Figure  G.  Comparison  of  the  radial  diffusion  coefficients  determined  in  this  paper  for  th' 
energy  range  0. 2-2.0  ke\'  (mottled  area)  with  Paper  1  (cross-hatched  area). and 
with  the  compilation  of  West  et  al.  1981  .  whose  individual  analysis  yields 
curve  1.  Curves  2-8  are  experimental:  2.Tomassian  et  al.  1972';  3,  Farley 
I9G9  :  4.  Newkirk  and  Walt  lytiSb  ,  5.  Lan?ero;ti  et  al.  1970,;  6.  Newkirk  and 
Walt  1968b  .  Faithammar  j970  ;  7.  Newkirk  and  Walt  1968a  ;  8.  Kavanaugh 
1968  .  Curves  9  14  are  theoretical  or  serni-theoretical:  9.  Nakada  and  .Mead 
1965  ;  10.  1  verskoy  1965  ;  11.  Birmingham  1969  ;  12.  13.  Cornwall  1968::  14. 
Hoi7W"rth  and  Mozei'  1979  . 

Figure  7.  Ltemonst ration  of  path  independence  in  projecting  the  electron  distributions 
from  on-trajectory  to  off-trajectory  points.  Panels  k  and  B  are  measured 
and  fitted  phase  spare  di.st ritiutions  at  the  trajectory  locations  shown  on  the 
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SCATHA  trajectory  panel.  Using  the  eigenfunction  representation,  the  distri¬ 
butions  are  projected  along  the  t  and  L  axes  respectively  to  the  off-trajectory 
point  C.  The  comparison  of  the  two  projected  distributions  is  shown  on  Panel 

C. 

Figure  8.  Projection  on  the  equatorial  plane  of  the  trajectory  of  a  10  MeV  electron  as 
it  makes  encounters  with  the  non-turbulent  magnetopause  indicated  by  the 
parabola.  The  electron  is  initially  bouncing  and  drifting  towards  the  magne¬ 
topause  on  an  T-shell  inside  the  magnetosphere.  Upon  encounter,  the  electron 
leaves  the  initial  T-shell  skips  along  the  equatorial  magnetopause  partially  in¬ 
side  and  partially  outside  of  the  magnetosphere  until  it  finds  itself  on  a  new 
magnetospheric  L-sheU  on  which  it  executes  its  regular  bounce-drift  motion. 
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